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2001-2005  Final  Report 


This  report  presents  the  specific  aims  and  accomplishments  of  our  breast  cancer  research  project  during 
the  years  of  funding  sponsored  by  the  US  Department  of  the  Army.  It  covers  our  activities  from  May  2001 
to  June  2005. 

1.  Introduction 

Tumor  oxygenation  is  believed  to  be  crucial  in  clinical  prognosis  as  evidenced  by  needle  electrodes.  It 
has  been  suggested  that  the  Blood  Oxygenation  Level  dependent  (BOLD)  method  can  provide  a  non- 
invaisve  means  for  tumor  oxygenation  monitoring.  However,  interpretation  of  BOLD  measurements  is 
confounded  by  variation  in  hemoglobin  concentration,  oxygen  saturation,  and  blood  flow.  Near  infrared 
(NIR)  spectroscopy  (NIRS)  can  differentiate  between  changes  in  hemoglobin  saturation  and 
concentration  in  breast  tumors  caused  by  respiratory  intervention. 

The  objective  of  this  study  is  to  investigate  vascular  oxygen  dynamics  in  breast  tumors  by  correlating 
the  results  of  NIRS  and  of  the  BOLD  method.  In  the  research,  we  wish  to  utilize  the  NIR  technique  for  a 
better  understanding  of  the  BOLD  method  in  order  to  improve  its  clinical  use  as  a  diagnostic  and/or 
prognostic  tool  for  breast  cancer  research  and  clinical  practice. 

Specifically,  we  have  proposed  two  hypotheses  and  three  aims,  as  follows: 

Hypothesis  7:  Three-channel  NIRS  combined  with  a  fiber  optical  needle  can  provide  spatial 
information  of  hemoglobin  saturation  (S02)  and  concentration  (Hb)  of  breast  tumors  during  intervention. 

Hypothesis  2.  Contrast  agents  in  breast  tumors  based  on  the  BOLD  effect  can  be  used  to  monitor 
variation  in  tumor  oxygenation. 

Aim  1\  To  investigate  heterogeneity  of  S02  in  the  tumor  vascular  bed  of  breast  tumors,  using  a  three- 
channel  NIR  system  against  the  fiber  optic  needle  measurements. 

Aim  2.  To  compare  and  correlate  the  measurement  results  of  the  breast  tumors  under  100%  02 
intervention  taken  simultaneously  from  the  three-channel  NIR  oximeter  and  from  the  BOLD  method. 

Aim  3:  To  study  the  influence  of  several  respiratory  interventions  on  S02  and  BOLD  signals  of  breast 
tumors  using  both  the  three-channel  NIR  oximeter  and  the  MRI  BOLD  method. 

2.  Body  of  the  report 

The  PI  along  with  her  research  collaborators  as  well  as  her  graduate  students  have  made  significant 
efforts  during  the  funding  period  to  accomplish  the  proposed  aims,  resulting  in  multiple  published 
journal  papers  and  manuscripts.  This  report  is  organized  in  such  a  way  that  we  will  use  our  published 
papers  and  manuscripts  to  provide  detailed  information  on  specific  aspects  for  each  research  aim, 
while  the  papers  and  manuscripts  are  organized  as  Appendices  1  to  4  for  referral,  as  listed  below: 

Appendix  1: 

Yueqing  Gu,  Vincent  Bourke,  Jae  Kim,  Anca  Constantinescu,  Ralph  P.  Mason,  Hanli  Liu,  “Dynamic 
Response  of  Breast  Tumor  Oxygenation  to  Hyperoxic  Respiratory  Challenge  Monitored  with  Three 
Oxygen-Sensitive  Parameters,”  Applied  Optics,  42(16),  2960-2967  (2003). 

Appendix  2: 

Hanli  Liu,  Yueqing  Gu,  Jae  G.  Kim,  and  Ralph  P.  Mason,  “Near  Infrared  Spectroscopy  and  Imaging  of 
Tumor  Vascular  Oxygenation,”  Methods  in  Enzymology:  Imaging,  ed.  by  Michael  Conn,  vol.  385-386,  pp. 
349-378  (2004).  '  “  ~ 
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Yueqing  Gu,  Ralph  Mason,  and  Hanli  Liu,  “Estimated  fraction  of  tumor  vascular  blood  contents  sampled 
by  near  infrared  spectroscopy  and  19F  magnetic  resonance  spectroscopy,”  Optics  Express  13(5),  1724- 
1733  (2005).  http://www.opticsexpress.org/abstract.cfm7URI-OPEX-13-5-1724. 
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Mengna  Xia,  Yueqing  Gu,  Hanli  Liu,  Vikram  Kodibagkar,  Anca  Constantinescu,  Ralph  Mason,  “Tumor 
oxygen  dynamics  measured  simultaneously  by  near-infrared  spectroscopy  and  MR  EPI  imaging,”  in 
Biomedical  Topical  Meetings  on  CD-ROM  (The  Optical  Society  of  America,  Washington,  DC,  2004), 
ThF33. 
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status,”  in  Biomedical  Topical  Meetings  on  CD-ROM  (The  Optical  Society  of  America,  Washington,  DC, 
2004),  FB6. 
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Aim  1:  To  investigate  heterogeneity  of  S02  in  the  tumor  vascular  bed  of  breast  tumors,  using  a  three- 
channel  NIR  system  against  the  fiber  optic  needle  measurements. 

For  Aim  1,  we  have  designed,  implemented,  and  refined  a  multi-distance,  optical  spectroscopic  system, 
which  allows  us  to  perform  simultaneous  measurements  of  three  oxygen-sensitive  parameters:  arterial 
hemoglobin  oxygen  saturation, Sa02,  oxygenated  hemoglobin  concentration,  Hb02  in  tumor 
vasculature,  and  tumor  oxygen  tension,  p02  in  response  to  hyperoxic  respiratory  challenge  in  rat 
breast  tumors.  The  effects  of  two  hyperoxic  gases,  oxygen  and  carbogen  (5%  C02  and  95%  02)  were 
compared,  by  use  of  two  groups  of  Fisher  rats  with  subcutaneous  13762NF  breast  tumors  implanted  in 
pedicles  on  the  foreback.  Also,  two  different  gas-inhalation  sequences  were  compared,  i.e.,  air— 
carbogen-air-oxygen-air  and  air-oxygen-air-carbogen-air.  The  results  demonstrate  that  both  of  the 
inhaled,  hyperoxic  gases  significantly  improved  the  tumor  oxygen  status.  Detailed  development, 
experimental  descriptions,  and  data  analysis  for  Aim  1  can  be  found  in  Appendices  1  and  2. 

Aim  2:  To  compare  and  correlate  the  measurement  results  of  the  breast  tumors  under  100%  02 
intervention  taken  simultaneously  from  the  three-channel  NIR  oximeter  and  from  the  BOLD  method. 

For  Aim  2,  we  have  performed  animal  measurements  on  Fisher  rats  with  subcutaneous  13762NF 
breast  tumors  implanted  in  pedicles  on  the  foreback  using  both  the  NIR  oximeter  and  19F  magnetic 
resonance  spectroscopy.  This  study  introduces  an  experimental  approach  to  estimate  percentage  of 
hemoglobin  content  and  volume  sampled  by  near  infrared  spectroscopy  (NIRS).  Carbogen  (5%  C02, 
95%  02)  respiratory  intervention  was  used  to  induce  physiological  changes  in  a  group  of  six  Fisher  rat 
breast  tumors.  Changes  in  total  hemoglobin  concentration,  [Hb]totai,  and  in  total  blood  volume,  AVT_b|00d, 
of  the  tumors  were  measured  by  NIRS  and  19F  magnetic  resonance  spectroscopy  of  perflubron, 
respectively.  The  ratio  of  [Hb]totai/  AVT_b|00d  was  used  to  calculate  the  fraction  of  hemoglobin  contents 
sampled  by  NIRS.  The  results  showed  that  the  mean  value  of  estimated  fractions  is  within  a  range  of 
15-30%  of  total  hemoglobin  content  in  the  tumor  tissues.  Based  on  the  results,  we  suggest  that  NIRS 
does  not  sample  the  entire  hemoglobin  volume  of  the  tissue  vasculature,  but  is  more  sensitive  to 
microvasculature.  This  study  helps  to  understand  the  blood  vascular  volume  sampled  by  NIRS,  and 
demonstrates  that  the  low  cost,  portable  NIRS  system  may  be  a  reliable,  non-invasive,  real-time, 
monitoring  tool  for  changes  in  tumor  blood  contents.  Refer  to  Appendix  3  for  more  details. 

Also,  we  have  demonstrated  the  ability  to  simultaneously  investigate  tumor  oxygen  dynamics  by  Near- 
infrared  spectroscopy  and  Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen 
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Mapping  in  breast  tumors,  with  close  correlations  between  the  two  techniques.  More  details  are  given 

in  Appendices  4  and  6. 


Aim  3:  To  study  the  influence  of  several  respiratory  interventions  on  S02  and  BOLD  signals  of  breast 
tumors  using  both  the  three-channel  NIR  oximeter  and  the  MR I  BOLD  method. 

For  Aim  3,  we  have  studied  dynamic  changes  of  oxygenated  and  deoxygenated  hemoglobin 
concentrations  in  response  to  hyperoxic  gas  interventions  on  rat  breast  tumors  simultaneously 
measured  by  near  infrared  spectroscopy  and  BOLD  (blood  oxygenation  level  dependent)  contrast  MR 
imaging.  We  have  demonstrated  the  ability  to  correlate  NIR  and  ^  MRI  investigations  simultaneously. 
This  preliminary  study  demonstrates  that  the  global  signals  of  breast  tumors  obtained  by  NIRS  and 
BOLD  MRI  in  response  to  respiratory  intervention  are  in  good  agreement.  These  two  techniques  may 
be  complementary  to  one  another  for  in  vivo  measurements  of  tumor  oxygenation  dynamics.  MRI  can 
be  applied  to  deep  tumors  with  special  pulse  sequences,  whereas  NIRS  can  be  used  to  monitor  acute 
responses  of  tumors  to  different  therapeutic  interventions.  The  correlation  of  NIRS  and  BOLD  MRI 
techniques  could  provide  a  new  prognostic  tool  for  assessing  efficacy  of  tumor  therapies  and  could  help 
us  better  interpret  BOLD  MRI  signals.  Furthermore,  the  heterogeneity  of  tumor  vascular  oxygenation 
dynamics  to  external  interventions  can  be  investigated  directly  by  BOLD  MRI,  while  the  development  of 
NIR  imaging  for  tumor  oximetry  has  also  great  potential  to  reveal  tumor  vascular  heterogeneities. 
Please  refer  to  Appendix  5  for  details. 


Furthermore,  during  the  course  of  conducting  this  proposed  research  project,  we  have  explored  new 
extra  findings,  as  listed  below: 

We  have  investigated  the  effect  of  photothermal  laser  irradiation  on  rat  breast  tumor  (DMBA-4)  vascular 
contents.  An  805-nm  diode  laser  was  used  in  our  experiment  with  a  power  density  ranging  from  0.32  to 
1.27  W/cm2.  The  dynamic  changes  of  oxygenated  hemoglobin  and  total  hemoglobin  concentrations, 
A[Hb02]  and  A[Hb]totai,  in  rat  tumors  during  photothermal  irradiation  were  noninvasively  monitored  by  a 
near-infrared  spectroscopy  system.  A  multichannel  thermal  detection  system  was  also  used 
simultaneously  to  record  temperatures  at  different  locations  within  the  tumors.  Our  experimental  results 
showed  that:  (1)  photoirradiation  did  have  the  ability  to  induce  hyperthermic  effects  inside  the  rat  breast 
tumors  in  a  single  exponential  trend;  (2)  the  significant  changes  (p<0.005)  of  A[Hb02]  and  A[Hb]totai  in 
response  to  a  low  dosage  of  laser  irradiation  (0.32  W/cm2)  have  a  single  exponential  increasing  trend, 
similar  to  that  seen  in  the  tumor  interior  temperature;  and  (3)  the  increase  in  magnitude  of  A[Hb02]  is 
nearly  two  times  greater  than  that  of  A[Hb]totai,  suggesting  that  photoirradiation  may  enhance  tumor 
vascular  oxygenation.  The  last  observation  may  be  important  to  reveal  the  hidden  mechanism  of 
photoirradiation  on  tumors,  leading  to  improvement  of  tumor  treatment  efficiency.  The  detailed  study  is 
reported  in  Appendix  7:  Yueqing  Gu,  Wei  R.  Chen,  Mengna  Xia,  Sang  W.  Jeong,  and  Hanli  Liu,  “Effect  of 
Photothermal  therapy  on  breast  tumor  vascular  contents:  non-invasive  monitoring  by  near  infrared 
spectroscopy,”  Photochemistry  and  Photobiology  81  (4),  1002-1009  (2005). 
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4.  Conclusions 


The  objective  of  this  study  was  to  investigate  vascular  oxygen  dynamics  in  breast  tumors  by  correlating 
the  results  of  NIRS  and  of  the  BOLD  method.  In  the  research,  we  have  utilized  the  NIR  technique  for  a 
better  understanding  of  the  BOLD  method  in  order  to  improve  its  clinical  use  as  a  diagnostic  and/or 
prognostic  tool  for  breast  cancer  research  and  clinical  practice. 

Throughout  the  study,  we  have  demonstrated  that  both  of  the  inhaled,  hyperoxic  gases  significantly 
improved  the  tumor  oxygen  status.  Arterial  hemoglobin  oxygen  saturation, Sa02,  oxygenated 
emoglobin  concentration,  Hb02  in  tumor  vasculature,  and  tumor  oxygen  tension,  p02  displayed  similar 
dynamic  response  in  rat  breast  tumors  to  two  hyperoxic  gases,  oxygen  and  carbogen  (5%  C02  and 
95%  02).  But  the  response  times  were  different:  the  fastest  for  arterial  Sa02,  followed  by  biphasic 
changes  in  tumor  vascular  Hb02,  and  then  delayed  responses  for  p02.  Both  of  the  gases  induced 
similar  changes  in  vascular  oxygenation  and  regional  tissue  p02  in  the  rat  tumors,  and  changes  in 
Hb02  and  mean  p02  showed  a  linear  correlation  with  large  standard  deviations,  which  presumably 
results  from  global  versus  local  measurements.  Indeed,  the  p02  data  revealed  heterogeneous  regional 
response  to  hyperoxic  interventions.  Although  preliminary  near-infrared  measurements  had  been 
demonstrated  previously  in  this  model,  the  addition  of  the  p02  optical  fiber  probes  provides  a  link 
between  the  noninvasive  relative  measurements  of  vascular  phenomena  based  on  endogenous 
reporter  molecules,  with  the  quantitative,  albeit,  invasive  p02  determinations. 

Also,  our  preliminary  study  investigated  dynamic  changes  of  oxygenated  and  deoxygenated 
hemoglobin  concentrations  in  response  to  hyperoxic  gas  interventions  on  rat  breast  tumors 
simultaneously  measured  by  near  infrared  spectroscopy  and  BOLD  (blood  oxygenation  level 
dependent)  contrast  MR  imaging.  We  have  demonstrated  the  ability  to  correlate  NIR  and  MRI 
investigations  simultaneously.  This  preliminary  study  demonstrates  that  the  global  signals  of  breast 
tumors  obtained  by  NIRS  and  BOLD  MRI  in  response  to  respiratory  intervention  are  in  good 
agreement.  These  two  techniques  may  be  complementary  to  one  another  for  in  vivo  measurements  of 
tumor  oxygenation  dynamics.  MRI  can  be  applied  to  deep  tumors  with  special  pulse  sequences, 
whereas  NIRS  can  be  used  to  monitor  acute  responses  of  tumors  to  different  therapeutic  interventions. 
The  correlation  of  NIRS  and  BOLD  MRI  techniques  could  provide  a  new  prognostic  tool  for  assessing 
efficacy  of  tumor  therapies  and  could  help  us  better  interpret  BOLD  MRI  signals.  Furthermore,  the 
heterogeneity  of  tumor  vascular  oxygenation  dynamics  to  external  interventions  can  be  investigated 
directly  by  BOLD  MRI,  while  the  development  of  NIR  imaging  for  tumor  oximetry  has  also  great 
potential  to  reveal  tumor  vascular  heterogeneities. 
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The  simultaneous  measurement  of  three  oxygen-sensitive  parameters  [arterial  hemoglobin  oxygen  sat¬ 
uration  (Sa02),  tumor  vascular-oxygenated  hemoglobin  concentration  ([Hb02]),  and  tumor  oxygen  ten¬ 
sion  (p02)]  in  response  to  hyperoxic  respiratory  challenge  is  demonstrated  in  rat  breast  tumors.  The 
effects  of  two  hyperoxic  gases  [oxygen  and  carbogen  (5%  C02  and  95%  02)]  were  compared,  by  use  of  two 
groups  of  Fisher  rats  with  subcutaneous  13762NF  breast  tumors  implanted  in  pedicles  on  the  foreback. 
Two  different  gas-inhalation  sequences  were  compared,  i.e.,  air- carbogen-air- oxygen-air  and  air- 
oxygen-air- carbogen-air.  The  results  demonstrate  that  both  of  the  inhaled,  hyperoxic  gases  signifi¬ 
cantly  improved  the  tumor  oxygen  status.  All  three  parameters  displayed  similar  dynamic  response  to 
hyperoxic  gas  interventions,  but  with  different  response  times:  the  fastest  for  arterial  Sa02,  followed  by 
biphasic  changes  in  tumor  vascular  [Hb02],  and  then  delayed  responses  for  p02.  Both  of  the  gases 
induced  similar  changes  in  vascular  oxygenation  and  regional  tissue  p02  in  the  rat  tumors,  and  changes 
in  [Hb02]  and  mean  p02  showed  a  linear  correlation  with  large  standard  deviations,  which  presumably 
results  from  global  versus  local  measurements.  Indeed,  the  p02  data  revealed  heterogeneous  regional 
response  to  hyperoxic  interventions.  Although  preliminary  near-infrared  measurements  had  been  dem¬ 
onstrated  previously  in  this  model,  the  addition  of  the  p02  optical  fiber  probes  provides  a  link  between 
the  noninvasive  relative  measurements  of  vascular  phenomena  based  on  endogenous  reporter  molecules, 
with  the  quantitative,  albeit,  invasive  p02  determinations.  ©  2003  Optical  Society  of  America 
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1.  Introduction 

It  is  widely  recognized  that  hypoxic  regions  in  solid 
tumors  may  limit  the  efficacy  of  nonsurgical  therapy, 
including  radiotherapy,  photodynamic  therapy,  and 
chemotherapy.1-4  Many  adjuvant  interventions 
have  been  tested,  including  simple  strategies  such  as 
breathing  hyperoxic  gases.5-7  However,  a  meta¬ 
analysis  of  some  10,000  patients  showed  only  a  mod¬ 
est  benefit,  and  this  benefit  was  restricted  to  specific 
tumor  types.8  It  is  thought  that  the  failure  of  such 
interventions  was  largely  due  to  the  inability  to  iden- 
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tify  those  patients  who  would  benefit.  Indeed,  there 
is  growing  emphasis  on  tailoring  therapy  to  the  indi¬ 
vidual  characteristics  of  each  patient’s  tumor.  Fur¬ 
thermore,  carbogen  (5%  C02  and  95%  02)  and  oxygen 
have  been  used  on  experimental  tumors  in  animals  as 
well  as  on  clinical  trials  in  patients  for  many 
years.9-10  But  the  therapeutic  benefits  of  the  two 
kinds  of  respiratory  hyperoxic  gases  are  diverse,  de¬ 
pending  on  the  tumor  types  and  individuals.11-13 
Accordingly,  accurate  assessment  of  tumor  oxygen¬ 
ation  at  various  stages  of  tumor  growth  and  in  re¬ 
sponse  to  interventions  may  provide  a  better 
understanding  of  tumor  development  and  may  serve 
as  a  prognostic  indicator  for  treatment  outcome,  po¬ 
tentially  allowing  therapy  to  be  tailored  to  individual 
characteristics. 

Various  techniques  have  been  developed  to  mea¬ 
sure  oxygen  tension  (p02)  or  vascular  oxygenation  of 
tumors.14  Many  methods  are  invasive,  and  those 
requiring  biopsy  preclude  dynamic  investigations. 
Optical  techniques  based  on  light  absorption  of  en¬ 
dogenous  chromophores,  e.g.,  near-infrared  (NIR) 
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spectroscopy  of  oxygenated  and  deoxygenated  hemo¬ 
globin,  are  entirely  noninvasive  and  allow  real-time 
monitoring  of  tumor  vascular  oxygenation.15-17 
However,  NIR  has  limited  spatial  resolution,  and  it 
remains  to  be  determined  whether  vascular  oxygen¬ 
ation  is  related  to  therapeutic  outcome.  Hitherto, 
quantitative  p02  has  been  shown  to  have  prognostic 
value,18-21  but  p02  represents  a  balance  between  ox¬ 
ygen  delivery  and  consumption.  Thus,  we  seek  to 
explore  the  interplay  of  vascular  and  tissue  oxygen¬ 
ation.  Electrodes  have  been  used  widely  to  study 
tumor  oxygen  dynamics  with  respect  to 
interventions,22-24  but  they  are  generally  limited  to  a 
single  location  and  small  probes  can  be  fragile.  We 
have  ourselves  recently  shown  a  correlation  between 
p02  and  AHb02  in  some  tumors,  but  we  noted  dis¬ 
tinct  heterogeneity,  and  thus,  the  global  NIR  mea¬ 
surements  were  not  always  related  to  local  p02.25 
Multiple  fiber-optic  probes  may  be  inserted  into  a 
tumor,26-28  and  we  have  now  investigated  correlation 
between  NIR  measurements  and  multiple  (three)  si¬ 
multaneous  p02  measurements. 

We  now  report  simultaneous  measurements  of 
three  oxygen-related  parameters,  i.e.,  arterial  hemo¬ 
globin  oxygen  saturation,  Sa02;  tumor  oxygenated 
hemoglobin  concentration,  [Hb02];  and  tumor  oxygen 
tension,  p02,  to  assess  dynamic  responses  of  rat 
breast  tumors  to  hyperoxic  gases.  Changes  in  tumor 
vascular  [Hb02]  were  measured  by  NIR  spectroscopy 
(NIRS)  using  a  photon-migration,  frequency-domain 
device;  changes  in  regional  p02  were  monitored  by  a 
fluorescence-quenched,  oxygen-sensing,  fiber-optic 
system  (FOXY);  the  arterial  Sa02  values  were  re¬ 
corded  with  a  fiber-based,  pulse  oximeter. 


FOXY  System 


Fig.  1.  Experimental  setup  for  simultaneous  oximetry.  The 
3-mm-diameter  fiber  bundles  of  the  NIRS  system  deliver  and  de¬ 
tect  the  laser  light  through  the  tumor  in  transmittance  geometry. 
PMT  represents  a  photomultiplier  tube.  I/Q  is  an  in-phase  and 
quadrature  phase  demodulator  for  retrieving  amplitude  and  phase 
information.  The  FOXY  system  comprises  three  fiber-optic 
oxygen-sensing  probes  that  are  inserted  into  different  regions  of 
the  tumor.  The  pulse  oximeter  probe  is  placed  on  the  hind  foot  of 
the  rat. 


pression)  in  the  middle  parts  of  the  tumors,  providing 
optimal  geometry  to  interrogate  deep  tumor  tissue. 

Based  on  modified  Beer-Lambert’s  law,29  changes 
in  oxygenated  and  deoxygenated  hemoglobin  concen¬ 
trations,  A[Hb02]  and  A[Hb],  due  to  respiratory  in¬ 
tervention  were  derived  from  the  measured 
amplitudes  at  the  two  wavelengths  and  calculated 
with  the  following  equations25: 


A[Hb02] 


—  10.63  log(^)758  +  14.97  log(^)785 

cT 


2.  Materials  and  Methods 


(1) 


A.  Near-Infrared  Spectroscopy  System  for  Measurement 
of  Changes  in  [Hb02] 

NIR  light  (700  to  900  nm)  has  considerable  tissue 
penetration  depth  (several  centimeters)  and  permits 
in  vivo  sampling  of  large  tissue  volumes  (e.g.,  human 
breast,  brain,  skeletal  muscle,  or  tumors),  since  pho¬ 
ton  transport  in  tissue  is  dominated  by  scattering 
rather  than  by  absorption.  Absorption  of  NIR  light 
by  the  oxygenated  and  the  deoxygenated  hemoglobin 
chromophores  may  be  used  to  determine  hemoglobin 
oxygenation  and  blood  concentration  changes.  As 
described  in  detail  previously,16’25  a  homodyne 
frequency-domain  system  (NIM,  Philadelphia,  Penn¬ 
sylvania)  was  used  to  monitor  the  global  changes  in 
oxygenated  and  deoxygenated  hemoglobin  concentra¬ 
tions,  A[Hb02]  and  A[Hb],  respectively,  in  rat  breast 
tumors  in  response  to  variations  in  inhaled  gas. 
Briefly,  the  light  from  two  NIR  laser  diodes  (758  nm 
and  785  nm)  was  coupled  into  a  bifurcated  fiber  bun¬ 
dle  and  illuminated  on  the  tumor,  and  the  transmit¬ 
ted  light  was  collected  and  propagated  to  a 
photomultiplier  tube  (Fig.  1).  The  fiber  bundles 
were  placed  on  the  surface  of  the  tumors  in  a  trans¬ 
mittance  mode  parallel  to  the  body  of  the  rat.  The 
fiber  tips  touched  firmly  on  the  skin  (without  com- 
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where  AB  and  AT  are  the  baseline  and  transient  am¬ 
plitudes  measured  from  the  NIR  system,  respective¬ 
ly;  d  is  the  source- detector  separation;  the  unit  for 
both  A[Hb02]  and  A[Hb]  is  millimolar  per  differential 
path-length  factor  (DPF);  and  the  DPF  is  for  tumor 
tissues.  As  demonstrated  in  our  previous  study, 
normalization  of  A[Hb02]  and  A[Hb]  to  their  maximal 
values  can  eliminate  the  effects  of  d  and  DPF  on  the 
results.25 


B.  Fiber-Optic  Oxygen-Sensing  System  for  Measurement 
of  Changes  in  p02 

Regional  p02  in  tumors  was  monitored  with  a  mul¬ 
tichannel,  fiber-optic,  oxygen-sensing  system  (FOXY, 
Ocean  Optics,  Inc.,  Dunedin,  Florida).30  Three 
fluorescence-quenched,  optical  fiber  probes  (AL300, 
tip  diameter  410  p,m)  were  inserted  into  different 
regions  of  the  tumors  (Fig.  1).  Probes  were  posi¬ 
tioned  so  that  at  least  one  was  in  a  poorly  oxygenated 
region  (low  baseline  p02)  and  at  least  one  in  a  well- 
oxygenated  region  (high  baseline  p02).  If  necessary, 
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the  probes  were  gently  moved  through  the  tumor 
until  such  diverse  regions  were  located.  In  some 
cases,  the  mean  p02  derived  from  the  three  individ¬ 
ual  measurements  is  presented.  Although  this  is  a 
commercial  system,  few  details  have  been  published 
previously,31  and  no  applications  to  in  vivo  tumor 
oximetry  have  been  published  to  our  knowledge. 
Light  from  a  pulsed  blue  LED  (475  nm)  was  coupled 
into  one  branch  of  a  bifurcated  optical  fiber  bundle 
and  propagated  to  the  probe  tip.  The  distal  end  of 
the  probe  is  coated  with  a  thin  layer  of  a  hydrophobic 
solgel  material,  where  an  oxygen-sensing  ruthenium 
complex  is  effectively  trapped.  Illumination  of  the 
ruthenium  complex  causes  fluorescence  at  —600  nm. 
If  the  excited  ruthenium  complex  encounters  an  ox¬ 
ygen  molecule,  the  excess  energy  is  transferred  to  the 
oxygen  molecule  in  a  nonradiative  transition,  de¬ 
creasing  or  quenching  the  fluorescence  signal.  The 
degree  of  quenching  correlates  with  the  oxygen  con¬ 
centration,  and  hence,  p02. 

The  fluorescence  response  of  the  ruthenium  crystal 
complex  is  highly  temperature  dependent,  so  to  ac¬ 
complish  probe  calibration  it  was  necessary  to  stream 
gases  of  known  oxygen  concentrations  (100%,  20.9%, 
10%,  2%,  and  0%)  through  a  cylindrical  water  jacket 
heated  to  37  °C.  Calibration  curves  were  automati¬ 
cally  calculated  by  means  of  the  vendor-supplied  soft¬ 
ware,  with  the  second-order,  polynomial  calibration: 

j  =  1  +  KIO]  +  K2{0f  (3) 

where,  I0  is  the  fluorescence  intensity  at  zero  oxygen 
concentration  (nitrogen),  I  is  the  measured  intensity 
of  fluorescence  at  a  pressure  of  oxygen,  [O]  represents 
the  oxygen  concentration  (related  to  p02),  K1  and  K2 
are  the  first-  and  the  second-order  coefficients  and  are 
automatically  supplied  by  the  curve-fitting  routine 
from  the  calibration  measurements. 

C.  Pulse  Oximeter  for  Measurement  of  Arterial  Sa02 

Arterial  Sa02  of  the  breast-tumor-bearing  rats  was 
also  monitored  with  a  fiber-optic  pulse  oximeter 
(Nonin  Medical,  Inc.,  Plymouth,  Minnesota)  placed 
on  the  hind  foot  of  the  rats.  The  system  consisted  of 
two  optical  fibers  used  for  delivering  and  receiving 
the  light.  The  tips  were  placed  on  either  side  of  the 
foot  in  transmission  mode. 

D.  Animal  Model 

Mammary  adenocarcinomas  13762NF  (originally  ob¬ 
tained  from  the  Division  of  Cancer  Therapeutics, 
NIH,  Bethesda,  Maryland)  were  implanted  in  skin 
pedicles32  on  the  foreback  of  adult  female  Fisher  344 
rats  (—150  g).  Once  the  tumors  reached  1-2  cm  di¬ 
ameter,  rats  were  anesthetized  with  150-p.l  ketamine 
hydrochloride  (100  mg/ml,  i.p.)  and  maintained  un¬ 
der  general  gaseous  anesthesia  with  1.3%  isoflurane 
in  air  (1  dm3/min).  Body  temperature  was  main¬ 
tained  at  37  °C  by  a  warm  water  blanket.  Tumors 
were  shaved  to  improve  optical  contact  for  transmit¬ 
ting  NIR  light.  The  tumor  diameters  along  the 


air  carbogen  air 
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Fig.  2.  Time  profile  of  the  three  oxygen-sensitive  parameters,  i.e., 
the  normalized  changes  of  tumor  A[Hb02],  the  mean  changes  of 
tumor  Ap02,  and  the  arterial  Sa02  with  respect  to  carbogen 
breathing  in  a  representative  13762NF  rat  breast  tumor  (No.  1,  3.2 
cm3). 


three  major  orthogonal  axes  (a,  b,  c)  were  measured 
with  calipers  and  volume  estimated  with  an  ellipsoid 
approximation  with  the  formula:  V  =  ('n/6)(abc). 

Two  groups  of  rats  (n  =  7  in  each  group)  were  used 
to  compare  the  effects  of  carbogen  and  oxygen  on 
vascular  oxygenation  of  breast  tumors.  Group  1  ex¬ 
perienced  the  gas-inhalation  sequence  of  air- 
carbogen-air- oxygen-air.  Group  2  was  exposed  to 
the  reverse  sequence  of  air- oxygen-air- carbogen- 
air.  Each  gas  was  maintained  for  20  min.  In  addi¬ 
tion,  the  FOXY  p02  probes  were  applied  to  five  rats 
from  Group  1,  and  the  dynamics  of  the  three  oxygen- 
related  parameters  were  measured  simultaneously. 

3.  Results 

A.  Dynamic  Responses  of  Three  Oxygen-Related 
Parameters  to  Carbogen  Intervention 
Typical  time  profiles  of  the  normalized  A[Hb02], 
mean  Ap02,  and  Sa02  in  response  to  carbogen  inter¬ 
vention  are  shown  for  a  representative  13762NF 
breast  tumor  (No.  1,  3.2  cm3)  in  Fig.  2.  When  the 
inspired  gas  was  switched  from  air  to  carbogen,  the 
Sa02  readings  increased  rapidly  and  significantly 
from  the  baseline  value  of  85%  to  the  maximum  of 
100%  within  2.5  minutes  ( p  <  0.0001).  The  normal¬ 
ized  A[Hb02]  showed  a  sharp  initial  rise  in  the  first 
minute  (p  <  0.0001)  followed  by  a  slower,  gradual, 
but  further  significant  increase  over  the  next  19  min 
(p  <  0.001).  Mean  Ap02  increased  rapidly  by  ap¬ 
proximately  50  Torr  within  8  min  (p  <  0.0005)  and 
also  continued  a  slower  and  gradual  increase  over  the 
next  12  min  {p  <  0.005).  Return  to  breathing  air 
produced  a  significant  decline  for  all  three  signals 

(p  <  0.0001). 

Sa02  and  p02  displayed  a  single-phase  dynamic 
behavior  in  response  to  carbogen  intervention, 
whereas  A[Hb02]  showed  an  apparent  biphasic  re¬ 
sponse.  These  dynamics  may  be  characterized  by 
time  constants  of  a  single-exponential  response.  For 
the  tumor  in  Fig.  3,  Sa02  had  the  fastest  response, 
with  a  time  constant  of  T(Sa02)  =  1.1  ±  0.2  min  (R  = 
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Fig.  3.  Dynamic  responses  of  the  three  oxygen-sensitive  param¬ 
eters  to  carbogen  intervention  in  a  rat  breast  tumor  (No.  1,  3.2 
cm3).  Single-exponential  curve  fitting  yielded  Sa02  = 
0.204{1  -  exp[-(t  -  20.02)/1.1]}  +  0.85  ( R  =  0.93),  A[Hb02]  = 
0.655{1  -  exp[-(f  -  20.36)/2.59]}  +  0.125  ( R  =  0.89),  and 
Ap02  =  42.68{1  -  exp[-(£  -  21.01)/4.56]}  +  16.66  ( R  =  0.98); 
biexponential  fitting  resulted  in  A[Hb02]  =  0.373{1  -  exp[— (t  — 
20.36)/0.61]}  +  0.648{1  -  exp[-(f  -  20.36)/21]}  ( R  =  0.97). 


0.93),  followed  by  [Hb02]  with  T(A[Hb02])  =  2.59  ± 
0.06  min  ( R  =  0.89),  whereas  Ap02  yielded  the  slow¬ 
est  response  T(Ap02)  =  4.56  ±  0.06  min  ( R  =  0.98). 
Time  constants  for  Group  1  are  listed  in  Table  1.  In 
every  case  T(Sa02)  <  T(A[Hb02])  T(Ap02),  based  on 

the  single-exponential  fitting.  No  apparent  relation 
between  the  time  constant  and  the  tumor  volume  was 
observed. 

It  is  clear  that  the  response  of  AHb02  is  not  well 
represented  by  a  single  exponential,  and  thus,  a 
double-exponential  expression  with  two  time  con¬ 
stants,  tx  and  t2,  was  also  used  (Fig.  3).  Comparison 
between  the  biexponential  fitting  for  A[Hb02]  and  the 
single-exponential  results  for  both  Sa02  and  Ap02  in 
the  first  five  rat  tumors  (Table  1)  shows  that  the  time 
constants  of  Sa02  (—1.2  ±  0.4  min)  are  similar  to 
those  of  the  first  phase  of  A[Hb02]  (—0.5  ±  0.2  min), 
whereas  the  second  phase  is  longer  and  highly  vari¬ 


able  (—14  ±  11  min).  No  significant  correlations 
were  found  between  any  of  the  time  constants  in 
Table  1. 

Time  delay,  td,  between  the  time  when  the  gas 
intervention  was  initiated  and  the  time  when  the 
changes  in  signals  were  detected,  reveals  another 
difference  among  the  three  oxygen-sensitive  param¬ 
eters.  For  tumor  1  (Fig.  2),  the  Sa02  signal  was  the 
first  to  respond  to  the  intervention.  Change  in 
A[HbOa]  was  observed  30  s  later  with  td  =  30  s, 
followed  by  changes  in  Ap02  another  30  s  later  (tr,  = 
60  s).  Similarly,  when  the  gas  was  returned  from 
carbogen  to  air,  the  Sa02  signal  decreased  immedi¬ 
ately,  followed  by  declines  in  A[Hb02]  and  in  Ap02 
with  td  of  30  and  120  s  later,  respectively.  As  ex¬ 
pected,  changes  in  Sa02  always  preceded  AHb02,  and 
Ap02  occurred  last  for  all  tumors. 

B.  Comparison  of  the  Effects  of  Carbogen  and  Oxygen 
Intervention  on  Tumor  Oxygenation 

Switching  from  air  breathing  to  carbogen  or  oxygen 
produced  similar  changes  in  AHb02  [Fig.  4(a)]. 
However,  the  time  course  was  substantially  different, 
requiring  a  biphasic  exponential  fit  for  carbogen,  but 
a  single  exponential  for  oxygen  [Fig.  4(b)].  For  the 
seven  tumors  in  Group  1,  there  was  no  significant 
difference  ( p  >  0.3)  in  the  maximum  magnitude  of 
A[Hb02]  caused  by  carbogen  or  oxygen  interventions 
[Fig.  4(c)]. 

To  examine  the  possible  effect  of  preconditioning 
required  that  Group  2  experience  a  reversed  gas  in¬ 
tervention,  with  exposure  to  oxygen  prior  to  carbogen 
[Fig.  5(a)].  In  this  case,  the  time  constants  of  the 
normalized  tumor  vascular  A[Hb02]  were  now  simi¬ 
lar  for  both  gas  challenges:  indeed,  for  six  of  seven 
tumors,  carbogen  no  longer  induced  the  biphasic  be¬ 
havior.  Figure  5(b)  shows  that  changes  in 
(AHb02)max  were  similar  to  those  in  Group  1,  and 
again,  the  two  gases  did  not  produce  significantly 
different  response  in  (AHb02)max.  This  is  empha¬ 
sized  for  both  Groups  1  and  2  by  a  strong  linear 
correlation  (slope  =  1.16)  between  the  A[Hb02]max 
values  observed  in  response  to  each  of  the  two  con- 


Table  1.  Time  Constants  of  Sa02,  A[Hb02],  and  Ap02  Response  to  Carbogen  and  Oxygen  Intervention  in  the  Breast  Tumors" 


Breast  Tumors 
Volume  (cm3) 

Single-Exponential  Fitting  of  Sa02,  A[Hb02]  and  Ap02 
(Carbogen  Intervention) 

Double-Exponential  Fitting  for 
A[Hb02]  (Carbogen 
Intervention) 

Single-Exponential 
Fitting  of  A[Hb02] 
(02  Intervention) 

Sa02 

A[Hb02] 

Ap02 

t  (min) 

R 

t  (min) 

R 

t  (min) 

R 

t1  (min) 

t2  (min) 

R 

t  (min) 

R 

No.  1  (3.2) 

1.1  ±  0.2 

0.93 

2.59  ±  0.06 

0.89 

4.56  ±  0.04 

0.98 

0.61  ±  0.03 

21  ±  3 

0.97 

0.35  ±  0.01 

0.92 

No.  2  (3.0) 

1.6  ±  0.2 

0.98 

3.40  ±  0.07 

0.91 

4.6  ±  0.1 

0.82 

0.62  ±  0.06 

11  ±  1 

0.96 

0.51  ±  0.01 

0.91 

No.  3  (4.6) 

1.2  ±  0.2 

0.97 

2.12  ±  0.06 

0.76 

2.26  ±  0.02 

0.98 

0.6  ±  0.1 

37  ±  3 

0.96 

1.52  ±  0.02 

0.89 

No.  4  (2.6) 

1.9  ±  0.3 

0.98 

2.68  ±  0.05 

0.93 

3.5  ±  0.1 

0.86 

0.12  ±  0.02 

5.2  ±  0.1 

0.98 

1.71  ±  0.03 

0.94 

No.  5  (5.6) 

0.8  ±  0.2 

0.91 

2.68  ±  0.05 

0.74 

4.51  ±  0.02 

0.99 

0.17  ±  0.03 

12.5  ±  0.6 

0.99 

5.49  ±  0.03 

0.98 

No.  6  (1.9) 

0.9  ±  0.2 

0.81 

1.62  ±  0.01 

0.95 

nd 

/ 

0.63  ±  0.08 

2.3  ±  0.1 

0.96 

5.16  ±  0.06 

0.93 

No.  7  (0.72) 

1.0  ±  0.5 

0.95 

3.60  ±  0.03 

0.93 

nd 

/ 

0.61  ±  0.02 

10.5  ±  0.3 

0.98 

3.54  ±  0.03 

0.95 

Mean 

1.2  ±  0.4 

2.7  ±  0.6 

4  ±  1 

0.5  ±  0.2 

14  ±  11 

2.5  ±  2 

“Under  the  inhalation  sequence  of  air-  carbogen-air-  oxygen-air. 
6nd,  not  determined. 
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Fig.  4.  (a)  Time  course  of  changes  in  tumor  vascular  A[Hb02]  for 
a  representative  13762NF  breast  tumor  from  Group  1  (No.  2,  3.0 
cm3)  with  respect  to  altering  inhaled  gas.  (b)  Respective  curve  fits 
for  the  carbogen  and  oxygen  interventions,  (c)  Average  maximum 
values  of  normalized  A[Hb02]  for  the  seven  breast  tumors  in  Group 
1. 


secutive  interventions  [Fig.  5(c)].  In  this  case,  non- 
normalized  data  are  shown  for  specific  comparison  of 
the  absolute  A[Hb02]max  produced  by  oxygen  and  car¬ 
bogen  for  each  of  the  tumors. 

C.  Tumor  p02 

The  FOXY  p02  probes  generally  indicated  distinct 
heterogeneity  in  p02.  Moreover,  response  to  the  hy- 
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Fig.  5.  (a)  Dynamic  changes  in  tumor  vascular  A[Hb02]  for  a 

representative  13762NF  breast  tumor  from  Group  2  (No.  9,  2.6 
cm3)  with  gas-inhalation  sequence  reversed  compared  with  Group 

1.  (b)  Average  maximum  values  of  normalized  A[Hb02]  for  Group 

2.  Gas-inhalation  sequence  reversed  compared  with  Group  1. 

(c)  Correlation  between  maximum  A[Hb02]  achieved  with  carbo¬ 
gen  inhalation  versus  that  with  oxygen  (R  =  0.97):  ♦,  carbogen 

first;  A,  oxygen  first. 


peroxic  gas  was  diverse:  those  probes  that  indicated 
apparently  well-oxygenated  regions  usually  showed  a 
large  and  rapid  response,  whereas  those  with  lower 
baseline  p02  often  showed  little  change  [Fig.  6(a)]. 
The  p02  responses  to  the  two  interventions  showed  a 
highly  consistent  behavior  at  each  individual  location 
[Fig.  6(b)].  There  was  also  a  distinct  correlation  be¬ 
tween  the  global  NIR  measurements  and  the  mean 
Ap02  (Fig.  7).  Because  of  heterogeneity  in  regional 
p02,  the  standard  deviations  of  the  mean  p02  values 
were  large. 
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Fig.  6.  (a)  Time  profiles  of  tumor  Ap02,  measured  with  the  three 
channels  of  the  FOXY  fiber-optic,  oxygen-sensing  system  with  re¬ 
spect  to  different  gas  inhalations  for  breast  tumor  No.  3  (4.6  cm3). 
The  mean  signal  for  the  three  channels  was  calculated  and  is 
plotted  by  the  thicker  trace,  (b)  Correlation  between  Ap02  at 
individual  locations  in  the  tumors  in  response  to  carbogen  or  oxy¬ 
gen  for  the  five  tumors  in  Group  1  ( R  >  0.8). 


4.  Discussion 

In  this  study,  we  have  simultaneously  measured  the 
arterial  Sa02,  the  global  changes  in  the  A[Hb02]  of 
tumor  vasculature,  and  the  regional  changes  in  the 
Ap02  of  tumor  tissue,  in  response  to  hyperoxic  (i.e., 
carbogen  and  oxygen)  gas  interventions  with  a  pulse 
oximeter,  an  NIRS  system  and  a  multichannel,  fiber¬ 
optic,  oxygen-sensing  system,  respectively.  All 
three  oxygen-sensitive  indicators  displayed  similar 


Fig.  7.  Correlation  between  mean  Ap02  and  A[Hb02]  for  the  five 
breast  tumors  (R  >  0.86):  ♦,  transition  from  air  to  carbogen;  A, 
transition  from  air  to  oxygen. 


dynamic  tendency  in  response  to  carbogen  interven¬ 
tion  (Fig.  2). 

The  simultaneous  measurements  demonstrate  the 
compatibility  of  the  NIRS  system  with  the  FOXY 
fiber-optic  oxygen-sensing  system,  without  interfer¬ 
ence.  Both  systems  are  relatively  inexpensive  and 
provide  real-time  measurements,  but  the  multichan¬ 
nel  FOXY  fiber-optic  system  monitors  Ap02  in  spe¬ 
cific  locations,  whereas  the  NIRS  system  provides 
global  measurements.  Whether  AHb02  determined 
with  this  methodology  will  be  a  clinically  useful  pre¬ 
dictor  for  tumor  response  to  oxygen-dependent  inter¬ 
ventions  and  therapies  remains  to  be  determined. 
However,  it  is  established  that  measurements  of  p02 
have  prognostic  value  in  the  clinic18-20  so  that  corre¬ 
lations  between  p02  and  NIR  measurements  would 
be  very  important. 

We  have  previously  applied  a  polarographic  oxygen 
electrode  simultaneously  with  NIR.25  However, 
that  study  provided  only  a  single  local  p02  value,  and 
in  some  cases  correlations  with  global  NIR  measure¬ 
ments  were  very  poor.  The  optical  fiber  system  used 
here  allows  multiple  locations  to  be  interrogated  si¬ 
multaneously.  The  device  can  be  expanded  to  many 
channels,  but  our  system  uses  four  channels.  Un¬ 
fortunately,  probes  are  fragile,  and  the  oxygen- 
sensitive  coating  on  the  tips  is  readily  damaged. 
Thus,  we  only  had  three  probes  available  for  this 
study.  Indeed,  fiber-optic  probe  fragility  is  a  well- 
recognized  problem,  and  our  previous  experience 
with  the  more  expensive  OxyLite  system  was  also 
restricted  to  three  channels  owing  to  probe  damage.26 

The  FOXY  system  (~$13k)  is  much  less  expensive 
than  the  OxyLite  (~$48k),  and  its  mode  of  action  is 
also  simpler,  detecting  fluorescent  signal  intensity 
rather  than  fluorescence  lifetime.  It  seems  capable 
of  measuring  p02  across  the  whole  range  of  atmo¬ 
spheric  oxygen  tensions  (0-760  Torr),  whereas  the 
OxyLite  is  restricted  and  becomes  very  insensitive 
above  approximately  100  Torr.  However,  our  expe¬ 
rience  shows  that  although  the  FOXY  system  pro¬ 
vides  precise  measurements  of  Ap02,  absolute  values 
of  p02  may  not  be  reliable.  We  continue  to  perform 
validation  experiments.  By  contrast,  the  OxyLite 
system  seems  to  give  very  accurate  p02  values. 

Our  experience  shows  that  the  FOXY  probes  are 
much  easier  to  use  than  electrodes,  particularly,  in 
terms  of  calibration  and  stability.  Since  the  probes 
are  fragile,  we  insert  them  into  tumors  through  a  fine 
needle  (25  gauge),  which  readily  punctures  the  sur¬ 
rounding  skin  and  penetrates  tough  fibrous  tissues. 
The  needle  is  then  backed  up  from  the  tip  to  facilitate 
measurements.  The  probes  often  require  a  few  min¬ 
utes  to  settle  at  a  stable  baseline  value,  but  then  show 
good  baseline  stability  and  rapid  response  to  inter¬ 
ventions  [Figs.  2  and  6(a)].  They  are  easily  moved 
within  the  tumor  to  locate  regions,  presenting  a  par¬ 
ticular  p02  of  interest,  e.g.  hypoxic  or  well  oxygen¬ 
ated.  In  the  search  for  appropriate  locations,  probes 
are  moved  forward  to  interrogate  fresh  tissues  rather 
than  in  reverse,  since  blood  may  pool  in  the  tracks 
owing  to  vascular  damage.  However,  we  observe 
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minimal  bleeding  on  removal  of  the  probes  from  the 
tumors. 

We  have  found  no  interference  between  the  NIR 
and  FOXY  instruments,  although  any  tumor  motion 
associated  with  moving  the  fiber  probes  can  alter  the 
optical  contact  of  the  NIR  optrodes,  and  thus,  alter 
apparent  AHb02.  Thus,  baseline  AHb02  is  deter¬ 
mined  once  the  fiber  probes  are  situated.  New  fiber¬ 
optic  probes  of  the  FOXY  system  have  a  thick  coating 
of  fluorescent  gel  and  a  black  covering,  but  this  wears 
with  use  and  gradually  allows  reception  of  the  NIR 
light.  Since  the  LEDs  of  the  two  systems  operate  at 
very  different  wavelengths,  viz.  475  versus  760  nm, 
there  is  no  interference  for  detection.  The  detection 
of  local  NIR  light  by  the  FOXY  probe  opens  the  ex¬ 
citing  possibility  of  detecting  regional  hemoglobin  ox¬ 
ygen  saturation.  We  believe  the  FOXY  optical 
probes  could  be  moved  within  the  tumor  to  map  the 
distribution  and  path  of  the  transmitted  NIR  light, 
helping  to  explore  and  validate  the  optical  character¬ 
istics  of  the  tumor.  This  can  simultaneously  provide 
a  correlation  between  local  AHb02  and  Ap02. 

In  this  study,  we  have  examined  a  much  larger 
group  of  rats  than  previously.16’25  We  have  now 
shown  rigorously  that  the  two  hyperoxic  gases  induce 
similar  changes  in  vascular  oxygenation  (NIR)  and 
regional  tissue  p02  (FOXY)  in  this  type  of  rat  breast 
tumor.  These  data  are  consistent  with  our  previous 
observations  using  19F  NMR  imaging  (FREDOM)33  in 
this  tumor  type  and  also  in  rat  prostate  tumors.34’35 
If  the  two  gases  are  indeed  equivalent  in  terms  of 
manipulation  of  tumor  oxygenation,  it  could  have 
great  therapeutic  benefit  since  the  popular  carbogen, 
which  is  in  use  in  clinical  trials,36  can  cause  respira¬ 
tory  discomfort. 

The  current  data  show  that  AHb02  and  Ap02  are 
correlated  (Fig.  7),  and  thus,  such  noninvasive  obser¬ 
vations  could  have  value  in  the  clinic.  The  major 
deficiency  in  our  current  NIR  approach  is  lack  of 
spatial  discrimination,  and  thus  efforts  to  implement 
NIR  imaging  will  be  of  great  value.  It  will  also  be 
interesting  to  correlate  other  measurements,  such  as 
blood-oxygen-level-dependent  (BOLD)  proton  mag¬ 
netic  resonance  imaging,  which  provide  high  spatial 
resolution,  but  which  are  sensitive  to  vascular  flow 
and  volume  as  well  as  oxygenation.37 

The  biphasic  response  of  AHb02  to  carbogen  is  in¬ 
triguing,  and  we  believe  it  represents  the  distinct 
vascular  compartments  of  arterioles  (high  flow)  and 
capillaries.  However,  the  change  to  monophasic  be¬ 
havior,  when  carbogen  is  administered  second,  re¬ 
quires  further  exploration;  in  the  future,  we  propose 
to  test  various  concentrations  of  oxygen  and  carbon 
dioxide  and  air  to  separate  the  components  of  the 
response.  The  carbogen  dioxide  component  of  car¬ 
bogen  is  known  to  be  vasoactive;  however,  the  specific 
effects  may  depend  on  tumor  type,  site  of  growth,  and 
other  factors.9’38 

In  terms  of  vascular  oxygen  delivery,  the  data  in 
Table  1  reveal  the  progressive  movement  of  oxygen: 
td  (Sa02)  <  td  (A[Hb02J)  <  td  (Ap02).  As  expected, 
switching  to  hyperoxic  gas  caused  the  systemic  arte¬ 


rial  Sa02  to  increase,  as  a  result  of  the  immediate 
combination  of  deoxyhemoglobin  with  oxygen.  The 
highly  oxygenated  blood  circulated  in  the  systemic 
vasculature  of  the  rats  (including  the  capillary  bed  of 
the  tumor  tissue),  resulting  in  a  delayed  increase  in 
[Hb02]  in  the  tumor  vasculature,  and  led  to  an  un¬ 
loading  of  oxygen  to  the  tumor  tissue.  For  the  biex¬ 
ponential  model  of  A[Hb02],  the  fast  component  has 
a  similar  time  constant  to  the  Sa02  measured  with 
the  pulse  oximeter  on  the  hind  leg,  strongly  suggest¬ 
ing  that  it  represents  arteriolar  oxygenation  in  the 
tumor.  In  this  study,  tumor  volumes  do  not  show 
any  direct  relation  with  time  constants  or  changes  of 
amplitude  in  response  to  hyperoxic  gas  interventions. 

It  is  increasingly  evident  that  oxygen  and  hypoxia 
play  important  roles  in  tumor  development  and  re¬ 
sponse  to  therapy.18  NIR  offers  an  attractive  non¬ 
invasive  means  of  investigating  tumor  oxygenation, 
particularly  in  terms  of  dynamic  response  to  inter¬ 
ventions,  but  we  had  previously  shown  a  potential 
mismatch  between  global  AHb02  and  local  Ap02.25 
The  data  presented  here  indicate  a  correlation  be¬ 
tween  the  global  NIR  measurements  and  mean  p02 
values  with  even  as  few  as  three  representative  loca¬ 
tions  per  tumor.  This  does  suggest  that  it  will  be 
important  to  develop  regional  NIR  measurements 
and  that  even  relatively  crude  mapping  could  reveal 
heterogeneity.  In  the  meantime,  we  believe  these 
studies  provide  further  evidence  for  the  value  of 
NIRS  to  explore  tumor  physiology. 

This  study  was  supported  in  part  by  the  Depart¬ 
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ICP.  Such  processes  include  hydrocephalus,  pseudotumor  cerebri,  intracra¬ 
nial  mass  lesions,  and  toxic-metabolic  encephalopathy  (where  a  depressed 
level  of  consciousness  may  or  may  not  correspond  to  an  increased  ICP). 
MR-ICP  measurement  may  prevent  unnecessary  invasive  monitoring  in 
an  increased  risk  setting.  Single  time  measurement  of  ICP  may  also  be 
helpful  in  the  evaluation  of  patients  with  possible  ventriculoperitoneal 
shunt  malfunction — particularly  young  children  who  present  with  nonspe¬ 
cific  complaints  and/or  who  are  unable  to  communicate  their  symptoms. 
The  MRI-based  technique  would  limit  exposure  to  harmful  radiation 
from  repeated  CT  scans.  Also,  a  finding  of  normal  ICP  by  the  MRI-based 
method  might  prevent  unnecessary  shunt  exploration. 

If  proven  accurate  and  reliable,  noninvasive,  MRI-based  measurement 
of  intracranial  compliance  and  pressure  and  TCBF  may  provide  a  valuable 
adjunct  to  the  evaluation  of  patients  with  various  neurologic  problems.  In 
the  setting  of  TBI,  a  finding  of  increased  ICP  may  provide  an  early  warning 
signal  for  a  physician  to  adjust  therapy.  This  tool  also  may  prove  valuable 
in  situations  where  external  factors  confound  clinical  evaluation  (i.e.,  drug 
intoxication),  where  “coma”  does  not  necessarily  reflect  an  increased  pres¬ 
sure  state  (i.e.,  diffuse  axonal  injury),  or  where  coexisting  medical  prob¬ 
lems  make  placement  of  an  invasive  monitor  too  dangerous  (i.e., 
bleeding  disorders).  In  each  case,  noninvasive  MR-ICP  measurement 
would  enhance  the  array  of  currently  available  evaluation  modalities,  pro¬ 
viding  a  novel  adjunct  to  patient  diagnosis  and  serial  monitoring  with 
limited  expense  and  morbidity. 


[17]  Near-Infrared  Spectroscopy  and  Imaging  of 
Tumor  Vascular  Oxygenation 

By  Hanli  Liu,  Yueqing  Gu,  Jae  G.  Kim,  and  Ralph  P.  Mason 
Introduction 

Current  Development  in  Optical  Methods  for  Vascular 
Hemoglobin  Oxygenation 

In  recent  years,  a  large  number  of  investigations  have  been  conducted 
in  both  laboratory  and  clinical  settings  to  noninvasively  monitor  tissue  vas¬ 
cular  oxygenation  using  near-infrared  (NIR)  spectroscopy  (NIRS)  and  im¬ 
aging.  Although  the  NIR  imaging  techniques  are  limited  by  their  spatial 
resolutions,  they  have  a  great  potential  to  be  developed  as  a  new  imaging 
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modality  because  of  their  capabilities  to  provide  functional  images.  NIR 
spectroscopy  and  imaging  research  has  been  mainly  focused  on  two  organs: 
(1)  the  brain  and  (2)  the  breast. 

The  NIR  studies  of  the  brain  include  detection  of  brain  injury  or 
trauma,1 2 3  determination  of  cerebrovascular  hemodynamics  and  oxygen¬ 
ation,2,3  and  functional  brain  imaging  in  response  to  a  variety  of  neurologic 
activations.4 5 6'7  NIR  functional  brain  imaging  increasingly  has  become  of 
great  interest  in  studying  hemodynamic  response  to  brain  activation.8  This 
is  mainly  because  the  optical  signals  of  the  NIR  techniques  are  able  to  non- 
invasively  penetrate  through  the  scalp  and  skull  of  an  adult  human  and  are 
sensitive  to  changes  in  the  concentration  of  oxygenated  (HbO)  and  deoxy- 
genated  hemoglobin  (Hb).  Although  it  is  difficult  to  obtain  very  accurate 
quantification  of  cerebral  HbO  and  Hb  concentrations  from  the  NIR  im¬ 
aging  techniques  due  to  rigorous  requirements  of  theory  and  boundary 
conditions,9 10 11  the  techniques  can  offer  relatively  accurate  measurements  of 
changes  of  HbO  and  Hb,  thus  providing  quantitative  changes  in  total  cere¬ 
bral  blood  volume  (CBV,  which  is  assumed  to  be  proportional  to  total 
hemoglobin  concentration,  HbT,  where  HbT  =  HbO  +  Hb). 

The  objective  of  NIR  breast  imaging  is  to  develop  a  novel  functional 
imaging  modality  for  early  breast  cancer  detection  and  diagnosis  beyond 
currently  available  techniques.  Various  efforts  by  several  research 
groups10-15  have  been  made  in  either  laboratory  or  clinical  studies.  For 
example,  the  research  groups  of  Paulsen  and  colleagues12  and  of  Jiang15 


1  S.  Gopinath,  C.  S.  Robertson,  R.  G.  Grossman,  and  B.  Chance,  J.  Neurosurg.  79, 43  (1993). 

2  C.  Cheung,  J.  P.  Culver,  K.  Takahashi,  J.  H.  Greenberg,  and  A.  G.  Yodh,  Phys.  Med.  Biol. 
46,  2053  (2001). 

3  G.  Zhang,  A.  Katz,  R.  R.  Alfano,  A.  D.  Kofinas,  P.  G.  Stubblefield,  W.  Rosenfeld, 
D.  Beyer,  D.  Maulik,  and  M.  R.  Stankovic,  Phys.  Med.  Biol  45,  3143  (2000). 

4  M.  Fabiani,  G.  Gratton,  and  P.  M.  Corballis,  J.  Biomedi.  Opti.  1,  387  (1996). 

5  R.  Wenzel,  H.  Obrig,  J.  Ruben,  K.  Villringer,  A.  Thiel,  J.  Bemarding,  U.  Dimagl,  and 
A.  Villringer,  J.  Biomed.  Opt.  1,  399  (1996). 

6  B.  Chance,  E.  Anday,  S.  Nioka,  S.  Zhou,  L.  Hong,  K.  Worden,  C.  Li,  T.  Murray, 
Y.  Ovetsky,  D.  Pidikiti,  and  R.  Thomas,  Opt.  Express.  2,  411  (1998). 

7  D.  A.  Boas,  T.  Gaudette,  G.  Strangman,  X.  Cheng,  J.  J.  A.  Marota,  and  J.  B.  Mandeville, 
Neuroimage  13,  76  (2001). 

8  D.  A.  Boas,  G.  Jasdzewski,  G.  Strangman,  J.  P.  Culver,  and  R.  Poldrack,  “OSA  Biomedical 
Topical  Meetings,  Technical  Digest,”  p.  307.  Miami,  FL,  April  7-10,  2002. 

9  J.  C.  Hebden,  E.  M.  C.  Hillman,  A.  Gibson,  N.  Everdell,  R.  M.  D.  Yusof,  D.  T.  Delpy,  S.  R. 
Arridge,  T.  Austin,  and  J.  H.  Meek,  “OSA  Biomedical  Topical  Meetings,  Technical 
Digest,”  p.  587.  Miami,  FL,  April  7-10,  2002. 

10  B.  Tromberg,  N.  Shah,  R.  Lanning,  A.  Cerussi,  J.  Espinoza,  T.  Pham,  L.  Svaasand,  and 
J.  Butler,  Neoplasia  2,  26  (2000). 

11 S.  Fantini,  S.  A.  Walker,  M.  A.  Franceschini,  M.  Kaschke,  P.  M.  Schlag,  and  K.  T.  Moesta, 
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have  developed  a  frequency-domain  (FD)  16-source,  16-detector  breast 
imager  and  have  reported  their  in  vivo  results  of  optical  properties  of  ab¬ 
normalities  from  female  volunteers  and  patients.  The  research  group  led 
by  Chance  has  employed  a  time-domain  (TD)  32-channel  imaging  system 
in  conjunction  with  magnetic  resonance  imaging  (MRI)  to  increase  specifi¬ 
city  and  sensitivity  of  breast  cancer  detection.13  Because  of  simplicity  and 
low  cost  in  comparison  to  the  FD  and  TD  imaging  systems,  continuous- 
wave  (CW)  NIR  breast  imaging  systems  also  have  been  developed  by  the 
groups  of  Barbour  et  al.14  and  Chance.16  These  systems  are  currently  under 
clinical  tests  for  better  breast  cancer  detection  and  diagnosis. 

Although  Hb,  HbO,  and  HbT  concentrations  and  light-scattering  prop¬ 
erties  of  tumors  may  be  different  from  those  of  surrounding  tissues,  the  op¬ 
tical  contrast  between  tumor  and  surrounding  tissue  is  about  2-3-fold  at 
most  in  absorption,  and  much  less  in  light  scattering.17  Thus  much  effort 
on  increasing  the  optical  contrast  between  tumor  and  healthy  surrounding 
tissues  also  has  been  made  using  fluorescence  imaging18  or  molecular 
beacons16  to  detect  and  diagnose  cancer  or  tumor  with  improved  sensitivity 
and  specificity. 

Flowever,  efforts  in  using  NIR  techniques  to  monitor  tumor  responses 
to  therapeutic  interventions19-21  and  therapy,  such  as  to  chemotherapy22 
or  photodynamic  therapy,23  are  very  limited  and  preliminary.  Moreover, 


12  T.  O.  McBride,  B.  W.  Pogue,  S.  Jiang,  U.  L.  Osterberg,  and  K.  D.  Paulsen,  Opt.  Lett.  26, 822 
(2001). 

13  V.  Ntziachristos  and  B.  Chance,  “SPIE  Proceedings  of  Optical  Tomography  and 
Spectroscopy  of  Tissue  III,”  Vol.  3597,  p.  565,  1999. 
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Technical  Digest,”  p.  468.  Miami,  FL,  April  7-10,  2002. 

16  B.  Chance,  “OSA  Biomedical  Topical  Meetings,  Technical  Digest”  p.  450.  Miami,  FL, 
April  7-10,  2002. 

17  J.  B.  Fishkin,  O.  Coquoz,  E.  R.  Anderson,  M.  Brenner,  and  B.  J.  Tromberg,  Appl.  Opt.  36, 
10  (1997). 
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19  H.  Liu,  Y.  Song,  K.  L.  Worden,  X.  Jiang,  A.  Constantinescu,  and  R.  P.  Mason,  Appl.  Opt. 
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8,  53  (2003). 
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few  reports  have  been  found  on  using  NIR  spectroscopy  and  imaging  as  a 
prognostic  tool  for  therapy  planning  and  optimization,  or  for  tumor  prog¬ 
nosis.  Based  on  the  existing  knowledge  and  development  of  NIR  spectros¬ 
copy  and  imaging  for  the  brain  and  breast,  we  were  motivated  a  few  years 
ago  by  the  possibility  of  using  NIRS  as  an  efficient,  real-time,  noninvasive 
means  to  monitor  tumor  vascular  oxygenation  during  respiratory  interven¬ 
tions.  The  recent  results  in  our  study  were  obtained  from  animal  breast  and 
prostate  tumors  in  vivo  with  a  one-channel  NIR  spectrometer.  The  data 
have  demonstrated  that  the  NIR  techniques  could  have  applications  as  a 
prognostic  means  accompanying  cancer  therapy.19,20,24 

Tumor  Oxygenation  as  an  Indicator  for  Tumor  Response 
to  Interventions 

The  presence  and  significance  of  tumor  hypoxia  has  been  recognized 
since  the  1950s.  Hypoxic  cells  in  vitro  and  in  animal  tumors  in  vivo  are 
documented  to  be  three  times  more  resistant  to  radiation-induced  killing 
compared  with  aerobic  cells.25  Recent  studies  have  shown  that  tumor  hyp¬ 
oxia  is  a  possible  prognostic  indicator  and  is  related  to  the  aggressiveness  of 
a  tumor,  the  clinical  stage,  and  poor  clinical  outcome.26-30  To  improve  the 
efficacy  of  oxygen-dependent  treatment  modalities,  one  possible  strategy  is 
the  reduction  of  tumor  hypoxia  by  raising  the  arterial  oxygen  partial  pres¬ 
sure,  Pa02  to  overcome  diffusion-limited  hypoxia.  Raising  the  arterial  p02 
by  breathing  hyperoxic  or  hyperbaric  gas  mixtures  could  be  effective.31,32 
However,  attempts  to  apply  increased  oxygen  breathing  in  the  clinic  have 
not  always  been  successful,  and  this  may  be  attributed  to  the  inability  to 
identify  those  patients  with  hypoxic  tumors.33  More  specifically,  it  may  be 


24  Y.  Gu,  V.  Bourke,  J.  G.  Kim,  A.  Constantinescu,  R.  P.  Mason,  and  H.  Liu,  Appl.  Opt.  42, 
2960  (2003). 

23  S.  M.  Evans,  W.  T.  Jenkins,  M.  Shapiro,  and  C.  J.  Koch,  Adv.  Exp.  Med.  Biol.  411,  215 
(1997). 

26  M.  Hockel,  K.  Schlenger,  B.  Aral,  M.  Mitze,  U.  Schaffer,  and  P.  Vaupel,  Cancer  Res.  56, 
4509  (1996). 

27  E.  K.  Rofstad  and  T.  Danielsen,  Br.  J.  Cancer  80,  1697  (1999). 

28  A.  W.  Fyles,  M.  Milosevic,  R.  Wong,  M.  C.  Kavanagh,  M.  Pintile,  A.  Sun,  W.  Chapman, 
W.  Levin,  L.  Manchul,  T.  J.  Keane,  and  R.  P.  Hill,  Radiother.  Oncol.  48,  149  (1998). 

29  D.  M.  Brizel,  S.  P.  Scully,  J.  M.  Harrelson,  L.  J.  Layfield,  J.  M.  Bean,  L.  R.  Prosnitz,  and 
M.  W.  Dewhirst,  Cancer  Res.  56,  941  (1996). 

30  J.  H.  A.  M.  Kaanders,  K.  I.  E.  M.  Wijffels,  H.  A.  M.  Marres,  A.  S.  E.  Ljungkvist,  L.  A.  M. 
Pop,  F.  J.  A.  van  den  Hoogen,  P.  C.  M.  de  Wilde,  J.  Bussink,  J.  A.  Raleigh,  and  A.  J.  van  der 
Kogel,  Cancer  Res.  62,  7066  (2002). 

31  S.  Dische,  Br.  J.  Radiol.  51,  888  (1979). 

32  O.  Thews,  D.  K.  Kelleher,  and  P.  Vaupel,  Radiother.  Oncol.  62, 77  (2002). 
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attributed  to  the  inability  to  identify  those  patients  who  would  benefit  from 
such  interventions.  Therefore  evaluation  of  tumor  oxygenation  distribu¬ 
tions  and  their  changes  during  various  stages  of  tumor  growth,  during 
therapeutic  interventions,  and  during  therapy  is  needed.  Such  evaluation 
can  provide  a  better  understanding  of  tumor  development  and  tumor  re¬ 
sponse  to  therapy,  potentially  allowing  therapy  to  be  planned  to  individual 
characteristics. 

Numerous  studies  on  tumor  oxygen  tension  (p02)  measurements  have 
been  conducted  in  recent  years  using  a  variety  of  methods,34"36  such  as  mi¬ 
croelectrodes,37  optical  reflectance,38  electron  spin  resonance  (EPR), 39-40 
or  MRI. 36,41-45  The  latter  two  methods  have  advantages  of  making  re¬ 
peated  measurements  of  p02  noninvasively,  and  MRI  has  the  further  ad¬ 
vantage  of  providing  dynamic  maps  of  p02,  which  can  reveal  tumor 
heterogeneity.46  Although  NIRS  does  not  quantify  p02,  it  can  indicate 
dynamic  changes  in  vascular  oxygenation  and  has  the  advantage  of  being 
entirely  noninvasive,  providing  real-time  measurements,  and  being  cost- 
effective  and  portable.  Our  recent  studies  have  revealed  the  need  for  NIR 
imaging  of  tumor  vasculatures  to  study  tumor  heterogeneous  response  to 
therapeutic  interventions  and  therapy.  In  the  following  sections,  we  briefly 
introduce  basic  algorithms  used  to  quantify  tumor  hemoglobin  oxygenation, 
followed  by  the  NIR  instrument  description,  system  calibration,  and  experi¬ 
mental  methods.  Then,  we  provide  several  representative  results  taken 


33  J.  Overgaard,  in  “Progress  in  Radio-Oncology  V”  (H.  D.  Kogelnik,  ed.),  p.  469.  Monduzzi 
Editore,  Bologna,  Italy,  1995. 

34  H.  B.  Stone,  J.  M.  Brown,  T.  Phillips,  and  R.  M.  Sutherland,  Radiat.  Res.  136,  422  (1993). 

35  R.  P.  Mason,  S.  Ran,  and  P.  E.  Thorpe,  J.  Cell.  Biochem.  87,  45  (2002). 

36  D.  Zhao,  L.  Jiang,  and  R.  P.  Mason,  Methods  Enzymol.  386,  378  (2004). 

37  O.  Thews,  D.  K.  Kelleher,  B.  Lecher,  and  P.  Vaupel,  Adv.  Exp.  Med.  Biol.  428,  123  (1996). 

38  F.  Steinberg,  H.  J.  Rohrborn,  K.  M.  Scheufler,  S.  Asgari,  H.  A.  Trost,  V.  Seifert,  D.  Stolke, 
and  C.  Streffer,  Adv.  Exp.  Med.  Biol.  428,  553  (1996). 

39  P.  E.  James,  J.  A.  O’Hara,  O.  Y.  Grinberg,  T.  Panz,  and  H.  M.  Swartz,  Adv.  Exp.  Med.  Biol. 
428,  97  (1996). 

40  J.  A.  O’Hara,  F.  Goda,  J.  F.  Dunn,  and  H.  M.  Swartz,  Adv.  Exp.  Med.  Biol.  411, 233  (1996). 

41  D.  Le,  R.  P.  Mason,  S.  Hunjan,  A.  Constantinescu,  B.  R.  Barker,  and  P.  P.  Antich,  Magn. 
Reson.  Imaging  15,  971  (1997). 

42  S.  Hunjan,  R.  P.  Mason,  A.  Constantinescu,  P.  Peschke,  E.  W.  Hahn,  and  P.  P.  Antich,  Int. 
J.  Radiat.  Oncol.  Biol.  Phys.  41,  161  (1998). 

43  D.  Zhao,  A.  Constantinescu,  L.  Jiang,  E.  W.  Hahn,  and  R.  P.  Mason,  Am.  J.  Clin.  Oncol.  24, 
462  (2001). 

44  D.  Zhao,  A.  Constantinescu,  C.-H.  Chang,  E.  W.  Hahn,  and  R.  P.  Mason,  Radiat.  Res.  159, 
621  (2003). 

45  D.  Zhao,  S.  Ran,  A.  Constantinescu,  E.  W.  Hahn,  and  R.  P.  Mason,  Neoplasia  5, 308  (2003). 

46  R.  P.  Mason,  A.  Constantinescu,  S.  Hunjan,  D.  Le,  E.  W.  Hahn,  P.  P.  Antich,  C.  Blum,  and 
P.  Peschke,  Radiat.  Res.  152,  239  (1999). 
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from  both  prostate  and  breast  tumors  under  hyperoxic  respiratory  inter¬ 
ventions,  using  both  one-channel  and  multichannel  NIR  systems,  as  well 
as  one-channel  and  three-channel  p02  fiber  optic  needle  probes.  At  the 
end,  we  wish  to  demonstrate  that  the  NIR  techniques  are  complementary 
with  tumor  p02  readings  and  can  be  used  as  a  new  prognostic  means  for 
cancer  therapy  prognosis  and  therapy  planning. 


Theory  and  Algorithms 

It  is  well  known  that  hemoglobin  concentrations  and  oxygen  saturation 
in  tissue  vasculature  can  be  determined  using  NIRS,  since  light  absorptions 
of  HbO  and  Hb  in  the  NIR  range  are  distinct.  As  with  our  previous 
work,19,20  we  assumed  that  HbO  and  Hb  are  the  major  significant 
absorbing  species  in  tissue  vasculature,  including  tumors,  within  the 
selected  NIR  range  of  700-900  nm.  Although  diffusion  theory  has  been  a 
well-accepted  theoretical  approach  to  mathematically  quantify  light-tissue 
interaction,6  7,10,47  no  analytical  solution  has  been  available  for  solid  tumors 
because  of  their  finite  size  and  high  heterogeneity.  While  our  future 
work  will  use  the  finite  element  method  to  obtain  numerical  solutions 
for  the  diffusion  equation,  our  current  approach  is  based  on  modified 
Beer-Lambert’s  law  to  account  for  light  scattering  in  tumor  tissue. 

For  a  nonscattering  medium,  Beer-Lambert’s  law  gives  rise  to  the 
following  expressions  for  the  relationships  between  the  absorption  coeffi¬ 
cient,  /xa,  and  the  extinction  coefficient  for  deoxyhemoglobin  (eHb)  and 
oxyhemoglobin  (eHbo): 

Ml58  =  2.3{£^[Hb]  +  4ibo[HbO]}  (1) 

M785  =  2.3  {e7H8b5  [Hb]  +  e7860  [HbO] }  (2) 

where  [HbO]  and  [Hb]  are  concentrations  of  HbO  and  Hb,  respectively, 
the  factor  of  2.3  results  from  the  different  definitions  of  /xa  and  t  in  relation 
to  the  incident  and  detected  optical  intensities.  The  conventional  defi¬ 
nitions  for  /ra  and  e  are  7  =  /0e~^aL  and  7  =  7olO_£CL,  respectively,  where 
70  and  7  are  the  incident  and  detected  optical  intensities  in  transmission 
measurement  from  a  nonscattering  medium,  C  is  the  concentration  of 
hemoglobin  measured  in  mM/L,  and  L  is  the  optical  pathlength  through 
the  medium  in  cm.  By  eliminating  7  and  70  in  these  two  expressions,  we 
arrived  at  a  relationship  of  =  2.3  eC. 


47  M.  S.  Patterson,  B.  Chance,  and  B.  C.  Wilson,  Appl.  Opt.  28,  2331  (1989). 
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Because  of  light  scattering  in  tissue,  Beer-Lambert’s  law  cannot  be 
applied  directly  to  biologic  tissue,  such  as  tumors.  By  taking  an  empiric  ap¬ 
proach  to  modify  Beer-Lambert’s  law,  we  arrived  at  //a  =  2.3  eC  = 
2.3/Lslog(/o//),  where  Ls  represents  the  optical  pathlength  in  a  scattering 
medium  and  is  no  longer  equal  to  the  physical  distance  between  the  source 
and  detector,  d.  In  a  highly  light-scattering  medium,  Ls  is  much  longer  than 
d.  Specifically,  in  the  tumor  study,  we  use  Io  and  I  to  represent  the  incident 
and  detected  light  intensities,  when  the  tissue  sample  is  without  and  with 
light  absorption  from  [HbO]  and/or  [Hb].  Then,  changes  in  absorption  co¬ 
efficient  of  the  tumor,  A/ra,  between  baseline  and  transient  conditions 
under  respiratory  intervention  can  be  expressed  as 

A^a  =  A^aT  -  FaB  =  2.31og(/Jg//;r)/L,s  (3) 

where  IB  and  IT  are  baseline  and  transient  light  intensities  of  the  measured 
optical  signals,  respectively. 

Using  the  transmitted  light  intensities  at  A)  =  758  nm  and  A2  =  785  nm 
and  manipulating  Eqs.  (1)— (3),  we  can  quantify  the  changes  of  tumor 
[HbO]  and  [Hb]  due  to  an  intervention  as  follows: 

A  [HbO]  =  -11.73  x  log(V/7~)  +  14.97  x  log(W7r)  (4) 

Ls  Ls 

A[Hb]  =  8.09  x  ■  -  6.73  x  (5) 

Ls  Ls 

where  L^1  and  Lf-  are  optical  path  lengths  between  the  source  and  de¬ 
tector  at  the  respective  wavelengths.  The  units  of  A[HbO]  and  A[Hb] 
in  Eqs.  (4)  and  (5)  are  m M.  The  constants  given  in  the  equations  were  com¬ 
puted  with  the  extinction  coefficients  for  oxygenated  and  deoxygenated 
hemoglobin  at  the  two  wavelengths  used.48 

In  principle,  and  L^2  in  Eqs.  (4)  and  (5)  are  variables,  depending  on 
the  actual  separation  of  source  and  detector,  as  well  as  the  optical  proper¬ 
ties  of  tumor.  Previously,  a  differential  pathlength  factor  (DPF),  (i.e.,  Ls  = 
d  x  DPF)  has  been  introduced  to  associate  Ls  with  d.49  The  DPF  values  of 
blood-perfused  tissues  have  been  studied  intensively  for  muscles50  and 


48  W.  G.  Zijlstra,  A.  Buursma,  and  W.  P.  Meeuwsen-van  der  Roest,  Clin.  Chem.  37,  1633 
(1991). 

49  D.  T.  Delpy,  M.  Cope,  P.  van  der  Zee,  S.  Arridge,  S.  Wray,  and  J.  Wyatt,  Phys.  Med.  Biol. 
33,  1433  (1988). 

50  M.  Ferrari,  Q.  Wei,  L.  Carraresi,  R.  A.  De  Blasi,  and  G.  Zaccanti,  J.  Photochem.  Photobiol. 
16,  141  (1992). 
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brains51  with  approximate  values  of  4-6  and  5-6,  respectively.  However, 
little  is  known  about  DPF  for  tumors,  although  a  DPF  value  of  2.5  has 
been  used  by  others.52  In  our  approach,  we  made  two  assumptions:  (1) 
DPF  is  the  same  for  both  A[HbO]  and  A[Hb]  at  785  nm  (i.e.,  DPF  = 
DPF™0  =  DPFji^).  This  assumption  was  based  on  the  fact  that  the  absorp¬ 
tion  difference  between  oxygenated  and  deoxygenated  blood  at  785  nm 
is  much  smaller  than  that  at  758  nm.  The  maximal  relative  error  caused  by 
this  assumption  in  tumor  oxygen  interventions  was  estimated  to  be  less 
than  12%,  and  detailed  justification  and  discussion  were  given  by  Liu 
et  al19  Furthermore,  we  associated  Ls  to  Ha  by  L  =  v/3/2dy//4//ia,  where 
n's  is  the  reduced  scattering  coefficient,  according  to  Sevick-Muraca  et  al.53 
and  Liu.54  The  second  assumption  was  /x'(758  nm)  =  /j' (785  nm)  because 
of  weak  wavelength  dependence  of  light  scattering  from  tissue.  With  all 
the  conditions  mentioned  and  further  derivation,  as  given  in  details  by 
Kim  et  al.20  we  arrived  at  the  final  expressions  for  A[HbO]  and  A[Hb]: 

-10.63  ■logf^a)  +  14.97 -log(  4a) 

A[HbO]  =  - - ^  (6) 


where  d  is  the  direct  physical  distance  between  source  and  detector. 
A[Hb]total  is  followed  by  adding  Eqs.  (6)  and  (7),  as 

A[Hb]total  =  A[HbO]  +  A[Hb]  (8) 

Equations  (6)  to  (8)  will  be  used  in  calculating  A[HbO],  A[Hb],  and 
A  [Hb] tolai  in  tissue  phantoms  and  tumors  during  gas  interventions  in  later 
sections. 

The  units  for  A[HbO],  A[Hb],  and  A[Hb]totai  in  Eqs.  (6)  to  (8)  are  mMI 
DPF,  which  still  depends  on  the  optical  properties  of  tumor  at  a  particular 
wavelength.  Because  our  focus  is  on  dynamic  changes  in  [HbO]  due  to  re¬ 
spiratory  challenges,  in  some  cases,  we  normalize  A[HbO]  at  its  maximal 


51  P.  van  der  Zee,  M.  Cope,  S.  R.  Arridge,  M.  Essenpreis,  L.  A.  Potter,  A.  D.  Edwards,  J.  S. 
Wyatt,  D.  C.  McCormick,  S.  C.  Roth,  E.  O.  R.  Reynolds,  and  D.  T.  Delpy,  Adv.  Exp.  Med. 
Biol.  316,  143  (1992). 

52  R.  G.  Steen,  K.  Kitagishi,  and  K.  Morgan,  J.  Neurooncol  22,  209  (1994). 

53  E.  M.  Sevick-Muraca,  B.  Chance,  J.  Leigh,  S.  Nokia,  and  M.  Maris,  Anal.  Biochem.  195, 330 
(1991). 

54  H.  Liu,  Appl.  Opt.  40,  1742  (2001). 
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value  (i.e.,  A[HbO]/A[HbO]max)  to  eliminate  the  effect  of  DPF  on  our 
results. 


Experimental  Methods 

One-Channel  and  Multichannel  NIR  System 

A  dual-wavelength  (at  758  and  785  nm),  one-channel  NIR  system 
(NIM,  Inc.,  Philadelphia,  PA)  uses  an  in-phase  and  quadrature-phase  chip 
(IQ  chip).  As  shown  in  Fig.  1,  the  NIR  system  starts  with  a  radiofrequency 
(RF)  source  to  modulate  the  light  intensities  of  two  laser  diodes  (LDX  and 
LD2)  at  140  MHz  through  a  time-sharing  system.  The  light  passes  through 
a  bifurcated  fiber  optic  probe,  is  transmitted  through  the  tumor  tissue,  and 
is  collected  by  a  second  fiber  bundle.  The  light  is  then  amplified  by  a  photo¬ 
multiplier  tube  (PMT),  demodulated  by  the  IQ  circuit,  and  filtered  by  a 
lowpass  filter  (LPF)  for  passing  only  the  direct  current  (DC)  components. 
The  signals  are  digitized  by  an  analog-to-digital  converter  (ADC)  and 
stored  in  a  laptop  computer.  The  measured  DC  signals  at  the  IQ  branches 


Reference  signal 


Fig.  1.  Experimental  setup  for  simultaneous  tumor  oximetry,  using  an  NIRS  system  with 
either  a  needle  electrode  or  a  three-channel  fiber-optic  FOXY  system.  The  NIR  system 
consists  of  two  laser  diodes  modulated  at  140  MHz,  two  fiber  bundles  for  light  delivery  and 
collection,  a  photomultiplier  tube  (PMT),  an  in-phase  and  quadrature-phase  demodulator  (IQ 
chip),  and  lowpass  filter  (LPF)  for  retrieving  amplitude  and  phase  information.  The  FOXY 
system  comprises  three  fiber-optic,  oxygen-sensing  probes  that  are  inserted  into  different 
regions  of  the  tumor.  Alternatively,  a  needle  electrode  probe  can  be  inserted  into  the  tumor 
for  single-channel  p02  reading. 
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Fig.  2.  Schematic  experimental  setup  for  three-channel  NIRS  experiments.  One  of  the 
three  detectors  was  placed  opposite  to  the  light  source  in  order  to  detect  light  in  transmission 
mode,  and  the  other  two  detectors  were  set  in  the  semireflection  mode  on  the  tumor  surface. 

at  the  measured  wavelengths,  /bc( A)  and  Qbc(A),  lead  to  the  quantification 
of  optical  amplitudes,  A  (A),  and  phase,  61(A),  that  have  passed  through  the 
tumor  tissues55: 


71(A)  = 

\A(^)dc  +  QWdc 

(9) 

0(A)  = 

arctanfe(A)M 

V  / 

(10) 

where  A  represents  the  respective  wavelengths  used  in  the  NIR  system. 
Then,  changes  in  light  intensity  through  the  tumor  caused  by  hyperoxic 
interventions  are  used  to  compute  changes  in  tumor  vascular  [HbO]  and 
[Hb],  Whether  the  readings  of  0(A)  can  be  linked  to  tumor  physiology 
remains  to  be  further  explored. 

Besides  a  one-channel  NIR  system,  we  have  also  used  a  multichannel 
NIR  system  (Fig.  2)  to  explore  and  reveal  intratumoral  vascular  heterogen¬ 
eity  by  having  simultaneous  readings  at  several  separations  of  the  source 
and  detectors.  Unlike  the  single-channel  system,  the  multichannel  system 
uses  two  laser  diodes  at  730  and  850  nm  with  constant  light  intensities 
(i.e.,  DC  light)  and  three  photo  detectors.  These  detectors  are  placed  along 
the  circumference  together  with  the  light  source,  with  three  different 
source-detector  separations  in  order  to  detect  different  depths  within  the 
tumor  volume,  as  demonstrated  in  Fig.  3. 

Following  the  same  procedures  as  those  described  in  the  previous 
section,  we  arrived  at  the  expression  of  A[HbO]  for  the  multichannel 
system  with  the  two  wavelengths  of  730  and  830  nm,  as  follows: 

55  Y.  Yang,  H.  Liu,  X.  Li,  and  B.  Chance,  Opt.  Eng.  36,  1562  (1997). 
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Light  source 


Fig.  3.  A  schematic  diagram  showing  the  locations  of  the  three  detectors  and  possible 
internal  tumor  volumes  interrogated  by  the  different  detectors. 


/  \  730  /  .  \  830 

-0.674  •  logfe)  +  1.117- log  (£) 

A[HbO]  = - K  ~d - Lll_  (11) 

where  d  is  the  direct  source-detector  separation  in  cm,  and  the  unit  of 
A[HbO]  is  mM/DPF. 

System  Testing  and  Calibration 

It  is  important  to  characterize  and  verify  linearity  of  the  NIRS  system 
before  any  animal  experiment  starts.  We  have  performed  two  critical 
system  tests:  (1)  to  quantify  the  output  voltage  range,  within  which  the 
NIR  system  has  a  linear  response,  and  (2)  to  perform  liquid  phantom  mea¬ 
surements  to  confirm  the  physiologic  readings  derived  from  the  optical 
NIRS  system. 

Linearity  Testing  Associated  with  Crosstalk  Between  the  Amplitude  and 
Phase.  In  principle,  changes  in  amplitude  of  an  FD  photon  migration 
system  should  be  independent  of  its  changes  in  phase  (i.e.,  the  crosstalk  be¬ 
tween  amplitude  and  phase  should  be  minimal  for  an  ideal  NIR  system): 
Flowever,  in  reality,  it  is  difficult  to  have  such  a  perfect  condition  for  an 
actual  instrument,  so  it  is  necessary  to  conduct  the  crosstalk  test.  We 
altered  the  optical  densities  (O.D.)  of  optical  filters  in  front  of  the  two  light 
sources  (758  and  785  nm)  to  attenuate  the  detected  amplitude,  so  we  could 
see  whether  changes  in  amplitude  and  phase  were  correlated.  The  test 
showed  that  phase-amplitude  crosstalk  exists,  if  the  measured  output  in¬ 
tensity  is  either  too  small  (<50  mV)  or  too  large  (>380  mV),  as  shown  in 
Fig.  4A  and  4B  for  758  nm  and  785  nm,  respectively.  We  also  plotted  O.D. 
values  of  the  filters  versus  log  (amplitude  of  the  output  signal)  to  see  the 
electrical  attenuation  of  the  system  output  in  response  to  the  optical 
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Fig.  4.  Relationships  between  the  output  amplitude  (mV)  and  phase  (in  degree)  of  the 
NIR  single-channel  system  at  (A)  758  nm  and  (B)  785  nm. 


Fig.  5.  Linearity  tests  for  the  NIR  system  between  the  optical  input  attenuation  and 
electrical  output  signals  at  758  nm  (open  circles)  and  785  nm  (solid  circles).  The  X  variable  is 
the  optical  density  (O.D.  =  log[/o//])  of  the  filters,  which  were  inserted  before  the  detector  to 
attenuate  the  light  intensity.  Cases  (A)  and  (B)  correspond  to  a  higher  (350  mV  for  758  nm) 
and  lower  (193  mV  for  758  nm)  output  voltage,  respectively.  Both  350  and  193  mV  were 
obtained  when  the  measurement  was  performed  without  any  optical  attenuation. 


attenuation.  It  was  found  that  as  long  as  the  output  is  in  the  range  100- 
350  mV,  good  linearity  is  achieved  for  the  system  response,  as  shown  in 
Fig.  5A  and  5B. 

Blood  Phantom  Study.  To  calibrate  the  NIRS  system  for  quantifying 
hemoglobin  concentrations,  we  used  human  blood  for  tissue  phantom  ex¬ 
periments.  Two  packets  of  Sigma  (St.  Louis,  MO)  P-3813  phosphate  buf¬ 
fered  saline  (PBS)  (pH  7.4)  powder  was  used  to  make  a  2-liter  buffer, 
and  100  ml  of  20%  intralipid  was  added  into  the  buffer  for  a  1%  intralipid 
solution,  as  shown  in  Fig.  6A.  Then  14  g  of  baking  yeast  was  mixed  with  tap 
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Gas  tube 


Fig.  6.  (A)  Experimental  setup  for  system  calibration  using  a  liquid  phantom  with  1% 
Intralipid  and  baker’s  yeast  in  saline  buffer.  The  NIRS  probes  were  placed  in  reflectance 
geometry,  while  the  gas  bubbler  was  placed  opposite  to  minimize  liquid  movement  effects. 
(B)  Simultaneous  dynamic  changes  of  A[HbO],  A[Hb],  and  A[Hb]toIai  in  the  phantom 
solution  measured  using  the  NIR  single-channel  system.  The  dark  solid  curve  is  for  A[HbO], 
the  lighter  solid  line  is  for  A[Hb]tolai,  and  the  gray  solid  line  shows  A[Hb]  in  the  phantom 
solution.  Blood  was  deoxygenated  by  the  yeast’s  oxygen  consumption  and  oxygenated  by 
bubbling  02.  During  the  oxygenation  and  deoxygenation  process,  A[Hb]t0ia]  was  maintained 
as  a  constant.20 
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water  before  being  added  into  the  2-liter  buffer  solution  to  deoxygenate 
the  tissue  phantom  solution.  Amount  of  blood  for  each  addition  to  the 
solution  was  3  ml  for  two  additions  after  the  yeast  was  well  mixed  in  the 
solution. 

When  the  blood  was  fully  deoxygenated  by  yeast,  pure  oxygen  was 
introduced  in  the  solution  to  oxygenate  the  blood.  After  the  blood  was  fully 
oxygenated,  oxygen  blowing  was  stopped  in  order  to  re-deoxygenate  the 
solution  by  yeast  again.  The  gas  tube  for  oxygen  delivery  was  placed  far 
from  the  NIRS  probes  to  minimize  any  liquid  movement  effects,  as  shown 
in  Fig.  6A.  Source  and  detector  probes  for  the  NIRS  system  were  placed  in 
reflection  geometry  with  a  direct  separation  of  ~3  cm.  The  solution  was 
stirred  constantly  with  a  magnetic  stirrer  to  maintain  homogeneity. 

After  getting  the  raw  amplitude  data  from  the  experiment,  we  used  Eqs. 
(6)  to  (8)  to  calculate  A[HbO],  A[Hb],  and  A[Hb]totai.  As  shown  in  Fig.  6B, 
A[Hb]totai  remained  constant,  as  expected,  whereas  A[HbO]  and  A[Hb] 
altered  in  opposite  directions  during  the  oxygenation  and  deoxygenation 
cycle.  This  demonstrates  that  the  two  assumptions  used  to  derive  Eqs.  (6) 
and  (8)  are  correct  and  necessary  to  compensate  the  differences  in  DPF 
caused  by  two  different  wavelengths.20  Notice  that  there  appears  to  be  a 
delay  between  the  gas  switch  and  responses  of  A[HbO]  and  A[Hb]  as  shown 
in  Fig.  6B.  This  delay  is  possibly  caused  by  the  entire  process  of  transporting 
oxygen  gas  from  the  tube  to  the  solution,  oxygenating/deoxygenating  the 
solution,  and  reaching  a  homogeneous  state  near  the  detectors. 

Animal  Models 

In  the  prostate  tumor  study,  Dunning  prostate  R3327-HI  rat  tumors 
(originally  obtained  from  Dr.  Peter  Peschke,  DKFZ  Heidelberg)  were  im¬ 
planted  in  pedicles  on  the  foreback  of  adult  male  Copenhagen  rats.56  In  the 
breast  tumor  study,  rat  mammary  adenocarcinomas  13762NF  (originally 
obtained  from  the  Division  of  Cancer  Therapeutics,  National  Cancer  Insti¬ 
tute  [DCT  NCI])  grown  in  the  hindlimb  of  adult  female  Fisher  344  rats 
(~200  g)  were  used.  Once  the  tumors  reached  approximately  1-2  cm  in 
diameter,  the  rats  were  anesthetized  with  0.2  ml  ketamine  hydrochloride 
(100  mg/ml;  Aveco,  Fort  Dodge,  IA)  and  maintained  under  general  gas¬ 
eous  anesthesia  with  isoflurane  in  air  (1.3%  isoflurane  at  1  dm3/min  air) 
through  a  mask  placed  over  the  mouth  and  nose.  Tumors  were  shaved 
to  improve  optical  contact  for  transmitting  light.  Body  temperature  was 
maintained  by  a  warm  water  blanket  and  was  monitored  by  a  rectally 

56  E.  W.  Hahn,  P.  Peschke,  R.  P.  Mason,  E.  E.  Babcock,  and  P.  P.  Antich,  Magn.  Reson. 

Imaging  11,  1007  (1993). 
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inserted  thermal  probe  (see  Fig.  1)  connected  to  a  digital  thermometer 
(Digi-Sense,  Model  91100-50;  Cole-Parmer  Instrument  Company,  Vernon 
Hills,  IL).  A  pulse  oximeter  (Model  8600;  Nonin,  Inc..  Plymouth,  MN)  was 
placed  on  the  hindfoot  to  monitor  arterial  oxygenation  (Sa02).  Tumor 
volume  V  (in  cm3)  was  estimated  as  V  =  47t/3 (L  +  W  +  H/6)J ,  where  L, 
W,  and  H  are  the  three  respective  orthogonal  dimensions. 

In  general,  the  source-detector  fiber  separation  was  about  1-1.5  cm  in 
transmittance  geometry,  and  thus  the  tumor  volume  interrogated  by  NIR 
light  can  be  estimated  as  follows.  The  radius  of  probe  is  0.4  cm,  which 
makes  the  volume  of  the  cylinder  0.5-0.75  cm3.  By  diffusion  approxima¬ 
tion,  the  optical  penetration  depth  is  about  one  third  of  the  direct  distance 
between  the  source  and  detector,  which  makes  an  additional  spherical 
volume  with  a  radius  from  0.5-0.75  cm.  Therefore  the  total  tumor  volume 
interrogated  by  NIR  light  can  be  estimated  in  the  range  of  1.0-2.5  cm3, 
depending  on  the  actual  source-detector  separation.20 


Measurements  of  p02  Using  Oxygen  Needle  Electrode  and  Fiber-Optic 
Fluorescence-Quenched  Probe 

Needle  Electrode.  We  used  a  needle-type  oxygen  electrode  for  local 
p02  readings  in  several  tumors.  Linear  two-point  calibrations  were  per¬ 
formed  with  air  (21%  02)  and  pure  nitrogen  (0%  02)  saturated  saline 
buffer  solutions  before  the  animal  experiments,  and  we  estimated  an  in¬ 
strumental  precision  of  2-3  mmHg.  After  calibration,  the  needle  electrode 
was  inserted  carefully  into  the  tumor  with  the  reference  electrode  being 
placed  rectally,  as  described  in  detail  previously.20  Both  the  sample  and  ref¬ 
erence  electrodes  were  connected  to  a  picoammeter  (Chemical  Microsen¬ 
sor;  Diamond  Electro-Tech,  Inc.,  Ann  Arbor,  MI)  and  polarized  at 
-0.75  V.  Measurement  points  of  p02  were  manually  recorded,  whereas 
the  NIRS  data  were  acquired  automatically.  Measurements  of  p02  and 
NIRS  were  initiated,  while  rats  breathed  air  for  ~10  min  to  demonstrate 
a  stable  baseline.  The  inhaled  gas  was  then  switched  to  carbogen  (5% 
C02  and  95%  02)  or  oxygen  for  15  min  and  switched  back  to  air. 

Fiber-Optic  Fluorescence-Quenched  Probe.  For  multiple  p02  readings, 
we  used  a  multichannel,  fiber-optic,  oxygen-sensing  system  (FOXY;  Ocean 
Optics,  Inc.).  We  calibrated  the  three  fluorescence-quenched,  optical  fiber 
probes  initially  and  then  inserted  them  into  different  regions  of  the  tumors 
to  monitor  the  changes  of  oxygen  tension  (Ap02)  in  response  to  the  re¬ 
spiratory  challenges  as  described  previously.24  The  probes  were  placed  in 
such  a  way  that  at  least  one  was  in  a  relatively  poorly  oxygenated  region 
(low  baseline  p02)  and  at  least  one  was  in  a  relatively  well-oxygenated 
region  (high  baseline  p02).  A  mean  ApQ2  of  the  tumor  was  obtained  by 


364 


GENERAL  METHODS 


[17] 


averaging  the  three  readings.  A  typical  FOXY  channel  used  a  pulsed  blue 
light-emitting  diode  (LED)  at  ~475  nm  and  was  coupled  into  one  branch 
of  a  bifurcated  optical  fiber  bundle.  Then,  the  fluorescent  light  was  propa¬ 
gated  to  the  FOXY  probe  tip.  Each  probe  had  two  300-/im-diameter 
optical  fibers  with  an  aluminum  jacket.  The  distal  end  of  the  probes  was 
coated  with  a  thin  layer  of  a  hydrophobic  sol-gel  material,  so  that  an 
oxygen-sensing  ruthenium  complex  was  effectively  trapped  and  protected 
from  water.  The  ruthenium  complex  at  the  probe  tips  was  excited  by  the 
blue  LED  and  emitted  fluorescence  at  ~600  nm.  If  the  excited  ruthenium 
complex  encountered  an  oxygen  molecule,  the  excess  energy  would  be 
transferred  to  the  oxygen  molecule  in  a  nonradiative  transition,  decreasing 
or  quenching  the  fluorescence  signal. 

The  fluorescence  response  of  the  ruthenium  complex  was  highly  tem¬ 
perature  dependent,  so  the  probe  calibration  was  accomplished  by 
streaming  gases  of  known  oxygen  concentrations  (100,  20.9,  10,  2,  and 
0%)  through  a  cylindrical  water  jacket  heated  to  37°.  Calibration  curves 
were  calculated  by  the  vendor-supplied  software,  using  the  second-order, 
polynomial  calibration: 

J=1  +K,[0]  +  K2[0]2  (12) 

where  I0  is  the  fluorescence  intensity  at  zero  concentration  (nitrogen),  I 
is  the  measured  intensity  of  fluorescence  at  a  pressure  p  of  oxygen,  [O] 
represents  the  oxygen  concentration  (related  to  p02),  and  Kr  and  K2  are 
the  first-  and  second-order  coefficients  and  are  automatically  supplied  by 
the  curve-fitting  routine  from  the  calibration  measurements.  After  the 
system  calibration,  the  oxygen  concentration  in  tumor-tissue  sample 
measurements  was  deduced  using  Eq.  (12). 

NIR  Measurements  of  Breast-Prostate  Tumors  Under  Interventions 

We  conducted  all  animal  experiments  in  a  darkened  room,  and  the 
measurements  were  initiated  while  the  rats  breathed  air  for  10  min  to  get 
a  stable  baseline.  The  inhaled  gas  was  switched  to  carbogen  (5%  C02  and 
95%  02)  or  pure  oxygen  for  15-20  min,  and  then  back  to  air  for  15  min. 
Sometimes,  repeated  cycles  were  taken  either  to  check  the  reproducibility 
or  to  minimize  preconditioning  effects.  Raw  amplitude  data  from  either 
the  one-channel  or  three-channel  detectors  were  recorded  simultaneously 
during  the  experiments,  processed,  and  displayed  through  the  computer  to 
obtain  tumor  A[ HbO]  and  A[Hb]totai.  The  dynamic  data  were  fitted  with  a 
single-exponential  and/or  double-exponential  expression,  using  Kaleida- 
graph  software  (Synergy  Software,  Reading,  PA),  to  obtain  time  constants 
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and  amplitudes.  Those  fitted  parameters  allowed  us  to  further  develop  a 
two-region  tumor  model,  which  will  be  discussed  in  a  later  section. 

In  many  cases,  we  performed  simultaneous  measurements  of  changes  in 
tumor  HbO  and  in  tumor  p02,  using  either  the  needle  electrode  or  the 
fiber-optic  FOXY  system  along  with  the  NIRS  readings,  as  shown  in 
Fig.  1,  for  both  breast  and  prostate  tumors. 


Results  and  Model  Development 

HbO  and  p02  Changes  in  Prostate  Tumors  Under  Carbogen  Intervention 

We  have  measured  relative  changes  of  [HbO],  [Hb]total,  and  tumor 
tissue  p02  from  several  Dunning  prostate  R3327-HI  tumors,  and  Fig.  7 
shows  two  representative  data  sets.  Figure  7A  shows  the  temporal  profiles 
of  A[HbO]  and  p02  in  a  Dunning  prostate  R3327-HI  tumor  (3.6  cm3) 
measured  simultaneously  with  the  NIRS  and  p02  needle  electrode  during 
carbogen  respiratory  challenge.  After  the  breathing  gas  was  switched  from 
air  to  carbogen,  A[HbO]  increased  rapidly,  whereas  A[Hb]totai  seemed  to 
have  much  smaller  responses.  In  this  case,  tumor  tissue  p02  increased  at 
a  much  slower  rate,  from  a  baseline  of  60  mmHg  to  ~80  mmHg  during 
the  entire  carbogen  intervention.  Figure  7B  is  obtained  from  another  pro¬ 
state  tumor  (3.1  cm3);  the  electrode  readings  showed  a  slower  p02  re¬ 
sponse  from  ~15  to  40  mmHg,  whereas  the  NIRS  response  was  biphasic, 
a  sharp  rise  in  A  [HbO]  followed  by  a  further  slow,  gradual  significant  in¬ 
crease  over  the  next  ~15  min.  In  this  case,  A[Hb]totai  showed  little  change, 
~2%  of  the  maximum  in  A[HbO].  The  biphasic  feature  of  A[HbO]  has 
been  a  commonly  observed  dynamic  characteristic  under  tumor  vascular 
oxygenation,  as  reported  previously.19,20,24 

HbO  and  p02  Changes  in  Breast  Tumors  Under  Carbogen  and 
Oxygen  Intervention 

For  a  representative  13762NF  breast  tumor  (3.2  cm3),  typical  time  pro¬ 
files  of  the  normalized  A  [HbO]  and  mean  Ap02  in  response  to  carbogen 
and  then  oxygen  intervention  are  shown  in  Fig.  8.  When  the  inspired  gas 
was  switched  from  air  to  carbogen,  the  normalized  A[HbO]  showed  a  sharp 
initial  rise  in  the  first  minute  ( p  <  0.0001),  followed  by  a  slower,  gradual, 
but  further  significant  increase  over  the  next  19  min  (p  <  0.001).  The  mean 
Ap02  profile  was  averaged  over  the  three  individual  Ap02  readings 
(FOXY),  and  it  increased  rapidly  by  about  50  torr  (=  mmHg)  within  8  min 
( P  <  0.0005)  and  also  continued  a  slower  and  gradual  increase  over  the 
next  12  min  (p  <  0.005).  Return  to  breathing  air  produced  a  significant 
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Fig.  7.  Simultaneous  dynamic  changes  of  A[HbO]  and  p02  in  R3327-HI  rat  prostate 
tumors  using  NIRS  and  p02  needle  electrode.  (A)  A  representative  tumor  (3.6  cm3)  showed  a 
rapid  NIR  response,  whereas  (B)  another  tumor  (3.1  cm3)  showed  a  clear  biphasic  feature  in 
A[HbO]  with  a  slow  p02  response.  The  unit  of  A[HbO]  is  mM/DPF,  where  DPF  is  equal  to 
the  optical  pathlength  divided  by  the  source-detector  separation.  Dotted  vertical  line  marks 
the  time  when  the  gas  was  changed. 
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Fig.  8.  Dynamic  responses  of  A[HbO]  and  mean  Ap02  to  hyperoxic  gas  intervention  in  a 
rat  breast  tumor  (3.2  cm3).  In  response  to  carbogen  breathing,  single-exponential  curve  fitting 
yielded  Ap02  =  42.68  jl  -  exp[-(t- 21.01  )/4.56])  +  16.66  (r  =  0.98),  and  biexponential 
fitting  resulted  in  AfHbO]  =0.373  (1  —  exp[— (t-20.36)/0.61]j  +  0.648  (1— exp[-(t— 20.36)/21][ 
( r  =  0.97).  The  smooth  black  curves  plotted  along  with  the  data  are  obtained  according  to  the 
previous  exponential  expressions,  respectively. 


decline  for  both  signals  (p  <  0.0001).  When  the  hyperoxic  intervention 
was  repeated  with  pure  oxygen,  the  dynamic  features  of  A[HbO]  and 
Ap02  remained  consistent,  except  that  the  biphasic  behavior  of  AfHbO] 
was  less  apparent.  A  detailed  investigation  on  AfHbO]  differences  in 
response  to  carbogen  and  oxygen  interventions  has  been  given  in  our 
earlier  study.24 

In  response  to  carbogen  intervention,  the  p02  profile  displayed  a 
single-phase  dynamic  behavior,  whereas  AfHbO]  showed  an  apparent 
biphasic  response.  These  dynamics  may  be  characterized  by  time  con¬ 
stants  of  single-  and  double-exponential  responses,  respectively.  A  single 
exponential  fitting  gives  rise  to  a  slow  Ap02  response  of  r(Ap02)  = 
4.56  ±  0.06  min  ( r  —  0.98).  A  double-exponential  expression  with  two  time 
constants,  r(  and  r2,  was  used  to  fit  the  normalized  AfHbO],  yielding  fast 
and  slow  time  constants  of  0.61  min  and  21  min,  respectively.  Figure  8  plots 
the  fitted  curves  along  with  the  experimental  data.  During  the  course  of  our 
study,  we  have  often  observed  biphasic  characteristic  of  AfHbO],  which 
motivated  us  to  develop  a  model  to  interpret  the  experimental  finding,  to 
be  described  in  a  later  section. 

The  advantage  of  using  FOXY  p02  probes  is  to  detect  distinct  hetero¬ 
geneity  in  tumor  p02.  The  individual  responses  at  different  locations  within 
the  tumor  to  the  hyperoxic  gas  were  diverse:  those  probes  that  indicated 
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Fig.  9.  Time  profiles  of  tumor  oxygen  tension  changes,  Ap02,  measured  with  three 
channels  of  a  FOXY  fiber-optic,  oxygen-sensing  system  with  respect  to  different  gas 
inhalations  for  a  breast  tumor  (4.6  cm3).  The  mean  signal  for  the  three  channels  was 
calculated  and  plotted  by  the  thicker  trace.  Modified  from  Gu  et  at ,24 


apparently  well-oxygenated  regions  usually  showed  a  large  and  rapid  re¬ 
sponse,  whereas  those  with  lower  baseline  p02  often  showed  little  change 
(Fig.  9).  By  studying  a  group  of  five  rats  with  both  carbogen  and  oxygen 
interventions,  we  obtained  a  distinct  correlation  between  the  maximal 
values  of  global  A[HbO]  and  the  mean  Ap02,  as  shown  in  Fig.  10.  Because 
of  heterogeneity  in  regional  p02,  the  standard  deviations  of  the  mean  p02 
values  were  quite  large. 

Tumor  Heterogeneity  Observed  by  a  Multichannel  NIRS  System 

To  study  tumor  heterogeneity,  the  multichannel  NIRS  system  also  has 
been  used  to  detect  changes  in  [HbO]  with  different  source-detector  sepa¬ 
rations.  Fig.  11A  shows  three  temporal  profiles  of  A[HbO]  obtained  from 
three  detectors  in  a  breast  tumor  (16.6  cm3),  with  a  source-detector  sepa¬ 
ration  of  1.5  cm  for  detector  1,  2.5  cm  for  detector  2,  and  2.8  cm  for  de¬ 
tector  3,  respectively  (Fig.  3).  The  measurement  uncertainties  are  plotted 
at  discrete  times  in  Fig.  11  A.  After  10  min  of  air  breathing  measurement 
as  the  baseline,  the  inhaled  gas  was  switched  from  air  to  carbogen, 
causing  a  sharp  increase  in  A[HbO]  (p  <  0.0001  after  1  min  from  gas 
switch),  followed  by  a  further  gradual,  but  significant,  increase  over  the 
next  15  min  (p  <  0.0001). 
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Fig.  10.  Correlation  between  mean  Ap02  and  AHbO  for  five  breast  tumors  (r  >  0.86)  with 
a  linear  relationship  of  ApOj  =  110.65  AHbO  +  4.1285.  Solid  circles  are  for  transition  air  to 
carbogen,  whereas  open  circles  are  for  transition  air  to  oxygen.24 


As  with  the  single-channel  NIRS  measurement,  the  biphasic  features 
appear  clearly  in  the  multichannel  results.  The  rising  parts  of  A[HbO]  from 
detectors  1, 2,  and  3  after  gas  switch  to  carbogen  are  shown  in  Fig.  11B-11D, 
respectively,  along  with  the  fitted  curves.  The  two  equations  used  for  the 
curve  fitting  are 

A[HbO]single  =  A{1  -  exp[-(t  -  t0)/r]}  (13) 

A[HbO]doublc  =  Aj{1  -  exp[-(t  -  t0)/ri]}  +  A2{1  -  exp[-(t  -  t0)/r2]} 

(14) 

respectively,  for  the  single-exponential  and  double-exponential  expres¬ 
sions.  It  is  clear  that  the  double-exponential  expression  gives  a  much  better 
fit,  as  confirmed  by  the  respective  chi  squared  (%2)  and  R 2  values.  Table  I 
summarizes  all  the  fitted  parameters  for  the  respective  detectors. 

Mathematical  Model  Development  to  Interpret  Data 

As  shown  in  the  last  section,  the  temporal  profiles  of  tumor  A[HbO] 
caused  by  respiratory  challenge  can  be  well  fitted  with  a  double-exponential 
equation,  represented  by  two  time  constants  (fast  and  slow).  To  understand 
these  time  constants  and  to  interpret  the  experimental  findings,  we  further 
developed  a  hemoperfusion  model,19  briefly  described  as  follows. 
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TABLE  I 

Summary  of  Vascular  Oxygen  Dynamics  from  Three  Detectors  from  Fig.  11 

Double-exponential  fitting  AHb02  =  Aj  (1  - 

-  exp(— t/n)]  +  A2[l 

-  exp(— t/r2)] 

Parameters 

Detector  1 

Detector  2 

Detector  3 

Separation:  d  (cm) 

1.5 

2.5 

2.8 

Ai  (mAf/DPF) 

0.037  ±  0.001 

0.0125  ±  0.0002 

0.0134  ±  0.0002 

tj  (min) 

0.24  ±  0.01 

0.30  ±  0.02 

0.27  ±  0.01 

A2  (mM/DPF) 

0.020  ±  0.001 

0.0130  ±  0.0004 

0.0060  ±  0.0002 

t2  (min) 

8.27  ±  0.72 

9.87  ±  0.84 

7.00  ±  0.58 

x2 

0.005 

0.001 

0.0005 

R2 

0.95 

0.96 

0.95 

71/72  =  Ai/A2 

1.85  ±  0.07 

0.96  ±  0.05 

2.23  ±  0.08 

n/r2 

0.029  ±  0.004 

0.030  ±  0.004 

0.039  ±  0.005 

/1//2  =  (Ai/A2)/(ri/T2) 

64  ±  9 

32  ±  5 

58  ±  9 

Following  an  approach  used  to  measure  regional  cerebral  blood  flow 
(rCBF)  with  diffusible  radiotracers,57-59  we  made  an  analogy  to  evaluate 
tumor  hemodynamics,  using  the  respiratory  intervention  gas  as  a  tracer. 

In  response  to  respiratory  intervention,  a  sudden  small  change  intro¬ 
duced  in  arterial  02  saturation  (Sa02)  results  in  an  increase  in  arterial 
HbO  concentration  (AHbOartery).  This  increase  in  HbOartery  can  be  con¬ 
sidered  as  an  intravascular  tracer.60  Following  Kety’s  method  and  assuming 
that  changes  in  dissolved  02  are  negligible,60  we  have 

,  /  a  TTLpjvasculatureX 

(AHbOvaseulature)  =  / 1  AHbOartery  -  - J  (15) 


57  S.  S.  Kety,  Pharmacol.  Rev.  3,  1  (1951). 

58  S.  S.  Kety,  Israel  J.  Med.  Sci.  23,  3  (1987). 

59  H.  Watabe,  M.  Itoh,  V.  Cunningham,  A.  A.  Lammertsma,  P.  Bloomfield,  M.  Mejia, 
T.  Fujiwara,  A.  K.  P.  Johes,  T.  Johes,  and  T.  Nakamura,  J.  Cerebr.  Blood  Flow  Metab.  16, 
311  (1996). 

60  A.  D.  Edwards,  C.  Richardson,  P.  van  der  Zee,  C.  Elwell,  J.  S.  Wyatt,  M.  Cope,  D.  T.  Delpy, 
and  E.  O.  R.  Reynolds,  J.  Appl.  Physiol.  75,  1884  (1993). 


Fig.  11.  (A)  Dynamic  changes  of  [HbO]  measured  at  three  detectors  from  a  rat  breast 
tumor  (volume:  16.6  cm3).  Dotted  vertical  lines  mark  the  points  when  the  gas  was  changed. 
The  rising  parts  of  A[HbO]  obtained  from  three  detectors  were  fitted  using  both  single¬ 
exponential  and  double-exponential  expressions,  and  (B)  to  (D)  are  from  detectors  1  to  3, 
respectively. 
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where  /  represents  blood  perfusion  rate,  and  7  is  defined  as  a  vasculature 
coefficient  of  the  tumor.  The  coefficient,  7,  is  the  ratio  of  HbO  con¬ 
centration  change  in  the  vascular  bed  to  that  in  veins  and  equals 
(AHbOvasculature)/(AHbOvein). 

In  Eq.  (15),  whereas  /  and  7  are  constants,  AHbOvasculature  is  a  time- 
dependent  variable.  In  principle,  AHbOvascula,ure  can  be  solved  rigorously 
given  a  constant  input,  H0,  for  AHbOartery  after  time  0.  Our  previous  study 
demonstrates  that  changes  in  arterial  HbO  (Sa02)  are  much  faster  than 
that  in  the  vascular  bed.19  Then  we  arrive  at  Eq.  (16): 

AHbOvasculature(r)  =  7  x  H0  x  (1  -  e-to)  (16) 

Equation  (16)  indicates  that  the  change  in  oxygenated  hemoglobin  concen¬ 
tration  in  tumor  vasculature,  AHbOvasculature(r),  depends  on  the  blood 
perfusion  rate,  /,  the  arterial  oxygenation  input,  H0,  and  the  vasculature 
coefficient  of  the  tumor,  7. 

Our  NIR  instrument  indeed  measures  changes  in  vascular  HbO  concen¬ 
tration  (i.e.,  A[HbO]'asLUlalure).  Then,  Eq.  (16)  can  be  used  to  interpret  the 
NIR  measurements,  since  it  gives  the  same  exponential  form  as  Eqs.  (13) 
and  (14),  used  to  fit  the  experimental  data.  The  comparison  between  Eqs. 
(13)  and  (16)  reveals  that  the  measured  time  constant  is  closely  associated 
with  the  blood  perfusion  rate,  /,  and  the  vasculature  coefficient  of  the 
tumor,  7  in  the  measured  volume.  Furthermore,  if  the  measured  volume  in¬ 
volves  two  distinct  regions,  then  the  detected  signal  will  involve  two  differ¬ 
ent  blood  perfusion  rates,  fi  and  /2,  and/or  two  different  vasculature 
coefficients,  7j  and  72.  In  this  case,  it  is  reasonable  to  assume  that  the  mea¬ 
sured  signal  results  from  both  of  the  regions,  as  illustrated  in  Fig.  12A. 
Consequently,  Eq.  (16)  can  be  modified  with  a  two-exponential  expression 
and  two  time  constants: 

AHbOvaseulalure(f)  =  7l  x  H0  x  (1  -  e~^')  +  l2  x  H0  x  (l  -  «-*/“») 

=  Aj  x  (1  -  +  A2  x  (1  - 

(17) 

where  fx  and  7,  are  the  blood  perfusion  rate  and  vasculature  coefficient 
from  region  1,  respectively,  f2  and  72  are  from  region  2,  and 
Ai  =  71  x  Ho  and  A?  =  72  x  Ho.  The  two  time  constants  are  equal  to 
t\  —  71//1  and  t\  =  72/ fii  respectively.  When  Ai,  A2,  and  the  two  time  con¬ 
stants  are  determined  from  our  NIR  measurements,  we  arrive  at  the  ratios 
for  the  two  vasculature  coefficients  and  the  two  blood  perfusion  rates: 

7i=Ai.  /1  =  A]  /A2 

72  A2  ’  fi  ti/t2 


(18) 
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B  Light  delivery 

Well-perfused 
region 

Poorly  perfused 
region 


Light  detectors 

Fig.  12.  (A)  Schematic  diagram  to  show  a  tumor  model  with  two  vascular  perfusion 
regions,  along  with  the  light  patterns  propagating  within  the  tumor  tissue.  “A”  represents  a 
portion  of  the  detected  signal  interrogating  the  well-perfused  region,  and  “B”  represents 
another  portion  of  the  detected  signal,  mainly  passing  through  the  poorly  perfused  region. 
Our  assumption  is  that  the  total  detected  signal  is  a  sum  of  A  and  B.  (B)  A  more  realistic 
tumor  tissue  model  with  the  center  poorly  perfused  (leading  to  signal  A)  and  the  periphery 
well  perfused  (leading  to  signal  B).  It  demonstrates  possible  tumor  volumes  interrogated  by 
the  multichannel  detectors.  The  total  detected  signal  is  still  a  sum  of  A  plus  B. 


In  this  way,  by  quantifying  these  two  ratios,  we  are  able  to  obtain  in¬ 
sight  into  the  tumor  vasculature  and  blood  perfusion.  For  example,  a  rather 
uniform  vascular  structure  will  result  in  a  ratio  of  71/72  near  1,  whereas  the 
coexistence  of  two  time  constants  possibly  reveals  two  mechanisms  of  re¬ 
gional  blood  perfusion  in  the  tumor.  A  large  time  constant  implies  a  slow 
perfusion  through  a  poorly  perfused  area,  whereas  a  small  time  constant  in¬ 
dicates  a  fast  perfusion  through  a  well-perfused  area.  In  the  meantime,  we 
can  determine  the  ratio  of  the  perfusion  rates  in  these  two  areas  after 
obtaining  the  ratios  of  amplitudes  and  time  constants  from  the  two- 
exponential  curve  fitting.  In  this  way,  we  have  developed  a  hemodynamic 
model  for  the  NIR  data  interpretation,  which  allows  us  to  associate  tumor 
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blood  oxygenation  dynamics  with  regional  blood  perfusion  and  vascular 
structures  of  the  tumor  within  the  measured  volume. 


Discussions 

NIRS  is  noninvasive  and  provides  a  real-time  assessment  of  changes  in 
tumor  vascular  hemoglobin  oxygenation.  In  this  chapter,  we  basically  pro¬ 
vide  demonstration  of  the  ease  and  utility  of  NIRS  studies  of  tumors.  By 
switching  the  inhaled  gas  from  air  to  carbogen,  the  NIRS  measurement 
produces  a  rapid  biphasic  elevation  in  A[HbO],  The  rapid  time  constant 
is  in  the  range  of  seconds  to  a  minute,  whereas  the  slow  component  (10 
to  50  times  slower)  continues  for  many  minutes  (see  Figs.  7B  and  8).  The 
elevated  oxygenation  process  is  completely  reversible  upon  returning  air 
breathing,  but  still  present  10  to  20  min  after  the  baseline  inhalation  in 
many  cases  (see  Figs.  7B  and  8).  The  high  reproducibility  of  results  suggests 
that  one  can  apply  repeated  interventions  to  explore  the  efficacy  of  inter¬ 
ventions  designed  to  alter  tumor  vascular  oxygenation  (e.g.,  vascular 
targeting  agents).  Our  recent  study  has  also  shown  that  tumor  response 
to  oxygen  was  much  more  rapid  and  fit  well  to  a  monoexponential  curve. 
For  almost  every  tumor  in  a  group  of  seven  breast  tumors,  the  time  to 
reach  80%  of  maximum  elevation  in  A[HbO]  was  longer  with  carbogen 
intervention  than  with  oxygen  breathing.24 

With  single-channel  and  multichannel  NIRS,  we  have  measured  rela¬ 
tive  changes  in  [Hb]totai  and  [HbO]  in  breast  and  prostate  rat  tumors  in  re¬ 
sponse  to  hyperoxic  respiratory  intervention.  We  observed  that  respiratory 
challenge  caused  A[HbO]  to  rise  promptly  and  significantly  in  both  breast 
and  prostate  tumors.  However,  the  total  concentration  of  hemoglobin  did 
not  always  behave  consistently.  The  dynamic  changes  of  tumor  oxygen¬ 
ation  can  be  modeled  by  a  two-exponential  expression  with  a  fast  and  slow 
time  constant.  Based  on  the  model,  we  suggest  that  the  NIRS  measurement 
can  “see”  two  vascular  mechanisms  during  tumor  oxygenation  under  hy¬ 
peroxic  gas  intervention.  By  fitting  the  two-exponential  model  with  the 
NIR  experimental  data,  we  can  determine  the  two  time  constants  and 
their  corresponding  amplitudes,  leading  to  the  relationships  between  the 
two  perfusion  rates  and  between  the  vascular  structures,  as  expressed  in 
Eq.  (18). 

With  this  tumor  model,  we  are  able  to  obtain  information  on  the  blood 
perfusion  of  tumor:  a  large  time  constant  usually  represents  a  slow 
blood  perfusion,  whereas  a  small  time  constant  indicates  a  fast  blood  perfu¬ 
sion.  A  combination  of  well-perfused  and  poorly  perfused  mechanisms 
in  tumor  vasculature  will  result  in  coexistence  of  two  time  constants. 
Indeed,  some  tumor  lines  have  been  reported  with  only  20-85%  of  vessels 
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perfused,61  and  it  is  known  that  tumor  structures  and  oxygen  distribution 
are  highly  heterogeneous.46,62  Therefore  it  is  highly  possible  that  our  NIRS 
readings  often  detect  a  mixture  of  both  well-perfused  and  poorly  perfused 
regions  in  the  tumor,  depending  on  the  locations  of  sources  and  detectors 
of  the  NIR  system. 

In  the  process  of  developing  the  model,  we  introduced  a  vasculature  co¬ 
efficient,  7.  We  expect  that  7  depends  on  (1)  oxygen  consumption  and  (2) 
capillary  density  of  the  tumor.  Further  studies  are  underway  to  quantify 
more  about  this  coefficient  and  to  confirm  our  speculations. 

As  an  example  to  demonstrate  the  direct  use  of  the  newly  developed 
model,  the  multichannel  NIR  results  are  fitted  with  the  model,  and  the  cor¬ 
responding  fitted  parameters  are  listed  in  Table  I.  It  is  reasonable  to  expect 
that  the  central  region  of  a  solid  tumor  is  close  to  necrosis  (i.e.,  a  poorly 
perfused  region),  whereas  the  tumor  peripheral  region  is  highly  vascular¬ 
ized.  Then,  Fig.  12B  can  be  used  to  schematically  illustrate  possible  tumor 
volumes  interrogated  by  the  multichannel  detectors.  Close  inspection  of 
Table  I  reveals  that  the  tumor  structure  ratios  of  71/72  from  the  three  re¬ 
spective  detectors  are  significantly  different,  whereas  the  perfusion  ratios 
of  fi/f2  from  detectors  1  and  3  are  almost  doubled  with  respect  to  that  from 
detector  2.  Besides  the  results  presented  here,  we  have  observed  similar 
AHbO  profiles  and  fitted  parameters  from  other  breast  tumors.  In  such  a 
way,  we  can  obtain  more  details  of  tumor  heterogeneity  by  having  more 
source-detector  pairs  in  the  measurements.  Indeed,  we  are  now  in  the  pro¬ 
cess  of  developing  a  multichannel  NIR  imaging  system,  along  with  the  im¬ 
aging  reconstruction  algorithm,  so  that  tumor  heterogeneity  can  be  studied 
with  better  spatial  resolution. 

The  time  constants  are  not  source-detector  separation-sensitive.  Equa¬ 
tions  (6),  (7),  and  (11)  demonstrate  that  AHbO  and  AHb  are  proportional 
to  lid,  where  d  is  the  source-detector  separation.  This  indicates  that  a  dif¬ 
ferent  d  value  will  only  stretch  or  compress  an  entire  temporal  profile  of 
AHbO,  but  it  does  not  change  the  transient  behavior  of  time  response. 
The  same  argument  can  apply  to  DPF.  Currently  we  have  grouped  DPF 
into  the  unit  of  [HbO]  for  simplicity.  If  the  DPF  value  is  larger  than  1, 
the  values  of  A[HbO]  and  A[Hb]lotal  will  decrease  by  a  factor  of  DPF. 
However,  this  decrease  does  not  affect  the  time  constants,  ti  and  T2  the 
dynamic  responses  of  tumor  A  [HbO]  to  respiratory  intervention. 


61  H.  J.  J.  A.  Bernsen,  P.  F.  J.  W.  Rijken,  T.  Oostendorp,  and  A.  J.  van  der  Kogel,  Br.  J. 
Cancer  71,  721  (1995). 

62  R.  P.  Mason,  P.  P.  Antich,  E.  E.  Babcock,  A.  Constantinescu,  P.  Peschke,  and  E.  W.  Hahn, 
Ini.  J.  Radiat.  Oncol.  Biol.  Phys.  29,  95  (1994). 
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The  simultaneous  measurements  of  A[HbO]  and  Ap02,  shown  in  this 
chapter,  demonstrate  the  compatibility  of  the  NIRS  system  with  the  needle 
electrode  and  the  FOXY  fiber-optic  oxygen-sensing  system,  without  inter¬ 
ference.  All  three  systems  are  relatively  inexpensive  and  provide  real-time 
measurements,  but  the  needle  electrode  and  the  three-channel  FOXY 
system  monitor  p02  (viz.,  Ap02)  in  specific  locations,  whereas  the 
NIRS  system  provides  global  measurements.  There  are  a  few  advantages 
of  the  FOXY  system  over  the  needle  electrode:  (1)  it  allows  multiple  loca¬ 
tions  to  be  interrogated  simultaneously  and  (2)  it  is  much  easier  to  use  than 
electrodes,  particularly  in  calibration  and  stability.  However,  the  FOXY 
probes  may  not  provide  accurate  absolute  p02  readings. 

It  is  well  known  that  measurements  of  tumor  p02  have  prognostic  value 
in  the  clinic.28,63  In  this  chapter,  we  have  shown  a  linear  correlation  between 
tumor  AHbO  and  mean  Ap02  (Fig.  10)  when  a  tumor  undergoes  a  hyper- 
oxic  respiratory  intervention,  and  thus  the  noninvasive  NIRS  monitoring 
could  have  a  potential  value  for  clinical  practice.  The  major  deficiency  in 
our  current  approach  is  lack  of  spatial  resolution,  so  implementation  of  an 
NIR  imaging  system  for  tumor  monitoring  is  our  current  effort. 

NIRS  provides  a  global  assessment,  in  contrast  to  blood  oxygenation 
level  dependent  (BOLD)-MRI,  which  can  provide  high-resolution 
images.64,63  Although  the  latter  approach  can  show  heterogeneity  in  both 
temporal  and  spatial  response,  the  results  are  often  summarized  to  show 
mean  values  only.  As  such,  an  a  priori  global  measurement  can  provide 
similar  insight  into  dynamic  tumor  physiology  or  drug  pharmacodynamics, 
while  being  more  cost-effective,  portable,  and  easier  to  implement  and 
operate.  For  example,  several  groups  of  researchers  have  previously  used 
NIR  studies  of  tumors  to  observe  the  influence  of  chemotherapy,52 
pentobarbital  overdose,52  ischemic  clamping,66  and  infusion  of  perfluoro- 
carbon  blood  substitute.67  All  of  these  investigations  demonstrate  potential 
versatility  of  the  NIRS  application  for  diverse  future  studies. 

Furthermore,  the  setup  in  our  NIRS  uses  transmission  mode  between 
the  source  and  detector,  which  we  believe  probes  large  and  deep  portions 
of  the  tumor,  including  the  periphery  on  each  side,  as  well  as  the  center  (see 
Fig.  3).  This  is  in  contrast  to  the  setups  used  by  most  other  investigators, 


63  M.  Hockel  and  P.  Vaupel,  J.  Natl.  Cancer  Inst.  93,  266  (2001). 

64  S.  P.  Robinson,  D.  R.  Collingridge,  F.  A.  Howe  et  al.,  NMR  Biomed.  A 2,  98  (1999). 

65  X.  Fan,  J.  N.  River,  M.  Zamora  et  al.,  Int.  J.  Radiat.  Oncol.  Biol.  Phys.  54,  1202  (2002). 

66  F.  Steinberg,  H.  J.  Rohrbom,  K.  M.  Scheufler,  T.  Otto,  and  C.  Streffer,  Adv.  Exp.  Med. 
Biol.  428,  69  (1996). 

67  H.  D.  Sostman,  S.  Rockwell,  A.  L.  Sylva,  D.  Madwed,  G.  Cofer,  H.  C.  Charles, 
R.  Negro-Villar,  and  D.  Moore,  Magn.  Reson.  Med.  20,  253  (1991). 
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who  apply  reflectance  mode,21,66,68  which  predominately  detects  the  behav¬ 
ior  of  the  peripheral  vasculature.  Indeed,  the  dynamic  results  measured 
from  breast  tumors  reported  by  Hull  el  al .21  are  consistent  with  the  fast 
component  that  we  have  observed  and  attributed  to  the  well-perfused 
regions  of  our  tumors. 

Given  the  evidence  for  intratumoral  heterogeneity  from  MRI46,69  and 
histology,70  it  will  be  important  to  advance  the  NIRS  system  to  an  NIR  im¬ 
aging  system  so  as  to  study  not  only  dynamic,  but  also  spatial  aspects  of 
blood  oxygenation  in  tumor  vasculature.  In  the  meantime,  we  believe  the 
preliminary  results  reported  here  are  proof  of  principle  for  NIR  imaging 
of  tumor  vascular  oxygenation,  laying  a  foundation  for  more  extensive 
studies  to  correlate  NIR  imaging  measurements  with  tumor  heterogeneity 
and  heterogeneous  responses  to  various  tumor  therapeutic  interventions 
and  treatments. 

Conclusions 

In  conclusion,  we  believe  that  NIRS  presents  a  new  potential  imaging 
modality  to  examine  tumor  vasculature  rapidly,  noninvasively,  and  cost- 
effectively.  Ease  of  implementation  and  operation  permit  rapid  application 
to  accessible  tumors  in  cancer  patients.  The  inherent  compatibility  of  fiber- 
optics  technology  and  light  with  other  modalities,  such  as  electrodes20,24 
and  MRI,71  will  facilitate  multiparametric  multimodality  investigations  of 
tumor  heterogeneity  and  vasculature  in  the  near  future. 

In  summary,  we  have  demonstrated  in  this  chapter  that  the  NIR  tech¬ 
nology  can  provide  an  efficient,  real-time,  noninvasive  means  for  monitor¬ 
ing  vascular  oxygenation  dynamics  in  tumors  during  hyperoxic  respiratory 
challenge.  Concentration  changes  in  HbO  measured  from  both  breast  and 
prostate  tumors  often  display  a  very  prompt  rise,  followed  by  a  gradual  per¬ 
sistence  throughout  the  intervention.  By  developing  a  hemoperfusion 
model  with  two  exponential  terms  and  fitting  the  model  to  the  AHbO  data, 
we  have  recognized  two  perfusion  mechanisms  in  tumor  vasculature  and 
quantified  the  ratios  of  the  two  perfusion  rates. 

Furthermore,  we  have  also  obtained  tumor  p02  measurements  using  a 
needle  electrode  or  multichannel,  fiber-optic,  FOXY  probe  in  simultaneous 
conjunction  with  the  noninvasive  NIRS  measurement.  The  comparative 

68  M.  Kragh,  B.  Quistorff,  and  P.  E.  Kristjansen,  Eur.  J.  Cancer  37,  924  (2001). 

69  S.  P.  Robinson,  F.  A.  Howe,  L.  M.  Rodrigues,  M.  Stubbs,  and  J.  R.  Griffiths,  Semin.  Radial. 

Oncol  8,  198  (1998). 

70  B.  M.  Fenton,  Radiother.  Oncol.  44,  183  (1997). 

71  G.  Gulsen,  H.  Yu,  J.  Wang,  O.  Nalcioglu,  S.  Merritt,  F.  Bevilacqua,  A.  J.  Durkin,  D.  J. 

Cuccia,  R.  Lanning,  and  B.  J.  Tromberg,  Technol.  Cancer  Res.  Treatment  1,  497  (2002). 
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results  exhibit  a  linear  correlation  between  AHbO  and  Ap02  of  the  tumors 
under  hyperoxic  gas  intervention,  suggesting  that  the  NIRS  approach  could 
have  a  good  potential  value  in  the  clinic.  Finally,  the  newly  developed 
tumor  hemodynamic  model  allows  us  to  reveal  tumor  heterogeneities  at 
different  tumor  locations  based  on  the  multichannel  NIRS  results.  Through 
this  chapter,  we  lay  a  foundation  for  an  NIR  imaging  technique  to  be  fur¬ 
ther  developed  to  facilitate  investigations  of  tumor  heterogeneity  and  vas¬ 
cular  perfusion.  Such  a  noninvasive  imaging  approach  can  enhance  our 
understanding  of  the  dynamics  of  tumor  oxygenation  and  the  mechanism 
of  tumor  physiology  under  baseline  and  perturbed  conditions. 


Acknowledgments 

This  work  was  supported  in  part  by  the  Department  of  Defense  Breast  Cancer  Research 
grants  BC990287  (HL)  and  BC000833  (YG),  and  NIH  R01  CA79515  (NCI)/EB002762 
(NIBIB)  (RPM).  We  are  grateful  to  Vincent  Bourke  for  his  collaborative  work  on 
multichannel  p02  measurements  and  Dr.  Anca  Constantinescu  for  her  assistance  with  all 
the  tumor  investigations.  We  also  gratefully  acknowledge  Dr.  Britton  Chance  for  his  technical 
support  on  the  multichannel  NIR  system. 


[  1 8]  Measuring  Changes  in  Tumor  Oxygenation 

By  Dawen  Zhao,  Lan  Jiang,  and  Ralph  P.  Mason 


Introduction 

Significance  of  pC>2  in  Oncology 

It  has  long  been  appreciated  that  hypoxic  tumor  cells  are  more  resistant 
to  radiotherapy.1  Indeed,  a  3-fold  increase  in  radio  resistance  may  occur 
when  cells  are  irradiated  under  hypoxic  conditions  compared  with  oxygen 
pressure  p02  >  15  torr  for  a  given  single  radiation  dose.  However,  recent 
modeling  has  indicated  that  the  proportion  of  cells  in  the  range  0-20  torr 
may  be  most  significant  in  terms  of  surviving  a  course  of  fractionated  radio¬ 
therapy.2  Certain  chemotherapeutic  drugs  also  present  differential  efficacy, 
depending  on  hypoxia.3,4  Increasingly,  there  is  evidence  that  hypoxia  also 
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3  B.  Teicher,  J.  Lazo  et  al..  Cancer  Res.  41,  73  (1981). 

4  A.  C.  Sartorelli,  Cancer  Res.  48,  775  (1988). 
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Abstract:  This  study  introduces  an  experimental  approach  to  estimate 
percentage  of  hemoglobin  content  and  volume  sampled  by  near  infrared 
spectroscopy  (NIRS).  Carbogen  (5%  0O2,  95%  02)  respiratory  intervention 
was  used  to  induce  physiological  changes  in  a  group  of  six  Fisher  rat  breast 
tumors.  Changes  in  total  hemoglobin  concentration,  A[Hb]totai,  and  in  total 
blood  volume,  AVT-biood,  of  the  tumors  were  measured  by  NIRS  and  19F 
magnetic  resonance  spectroscopy  of  perflubron,  respectively.  The  ratio  of 
A[FIb]tota/ AV T-biood  was  used  to  calculate  the  fraction  of  hemoglobin 
contents  sampled  by  NIRS.  The  results  showed  that  the  mean  value  of 
estimated  fractions  is  within  a  range  of  15-30%  of  total  hemoglobin 
content  in  the  tumor  tissues.  Based  on  the  results,  we  suggest  that  NIRS 
does  not  sample  the  entire  hemoglobin  volume  of  the  tissue  vasculature,  but 
is  more  sensitive  to  microvasculature.  This  study  helps  to  understand  the 
blood  vascular  volume  sampled  by  NIRS,  and  demonstrates  that  the  low 
cost,  portable  NIRS  system  may  be  a  reliable,  non-invasive,  real-time, 
monitoring  tool  for  changes  in  tumor  blood  contents. 
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1.  Introduction 

In  biomedical  research,  optical  spectroscopy  and  imaging  can  potentially  provide  rapid, 
economical,  and  non-invasive  diagnostic  links  between  crucial  tissue  characteristics  and 
cancer  detection  or  diagnosis  [1-7].  Near  infrared  (NIR)  light  in  the  range  of  700  nm  to  900 
nm  has  an  optimal  penetrating  depth  in  tissue  (up  to  15  centimeters)  since  photon  transport  in 
tissue  is  dominated  by  light  scattering  rather  than  absorption  in  the  NIR  region.  Therefore,  the 
use  of  NIR  light  enables  us  to  sample  large  tissue  volumes,  such  as  the  breast,  the  brain,  and 
skeletal  muscles  as  well  as  tumors  deep  in  the  tissue.  In  the  NIR  spectral  range,  two 
important  endogenous  tissue  chromospheres,  i.e.,  oxygenated  and  deoxygenated  hemoglobin 
(Hb02  and  Hb,  respectively),  display  oxygen-dependent  absorption  features  [8,9]. 
Experimental  determination  of  absorption  quantities  of  tissue  can  provide  quantification  of 
several  physiological  parameters,  such  as  hemoglobin/myoglobin  concentrations,  vascular 
oxygenation,  and  blood  volume  changes.  Such  quantification  could  be  important  for  cancer 
prognosis  and  treatment  monitoring,  as  we  have  demonstrated  in  our  previous  reports  [10-12]. 

Because  of  large  scattering  of  NIR  light  in  tissue,  however,  it  is  not  clear  how  much  of  the 
total  vascular  volume  within  a  tissue  is  sampled  by  the  NIR  signals.  It  has  been  reported  that 
the  NIR  detection  is  most  sensitive  to  microvessel  density,  i.e.,  arterioles,  capillaries,  and 
venules  [13-14].  The  opaque  microstructures  in  tissue  also  have  contributions  to  the 
absorption  [15].  The  measured  absorption  is  actually  a  volume -weighted  average  of 
absorption  of  the  sampled  vessels  and  the  surrounding  tissue.  Thus,  the  measured  (apparent) 
chromophore  concentration  may  be  less  than  the  actual  concentration.  Moreover,  the  blood 
flow  measured  by  NIR  spectroscopy  (NIRS)  was  reported  2~3  times  lower  than  the 
plethysmographic  flow  [16-17].  These  reports  motivated  us  to  investigate  what  fraction  of 
blood  volume  in  tissue  is  sampled  by  NIR  .  In  particular,  tumor  tissue  is  recognized  as  highly 
heterogeneous.  Better  understanding  and  accurate  determination  of  tumor  vascular  blood 
volume  sensed  by  the  NIRS  could  provide  crucial  information  for  assessing  tumor  progress 
and  guiding  tumor  therapy  and  treatment. 

Several  investigators  have  shown  that  vascular  volume  can  be  estimated  using  19F 
Magnetic  Resonance  spectroscopy  (MRS)  following  infusion  of  a  perfluorocarbon  blood 
substitute  emulsion  [18-19].  Our  previous  studies  have  demonstrated  that  NIRS  can  be  used 
to  monitor  changes  in  total  (i.e.,  oxygenated  and  deoxygenated)  hemoglobin  concentrations 
within  tumor  vasculature.  Thus,  correlations  between  MRS  and  NIRS  measurements  may 
allow  us  to  investigate  the  relationship  between  the  blood  volumes  sampled  by  these  two 
respective  methods.  Specifically,  in  this  paper,  we  will  report  our  measurements  of  changes  in 
total  hemoglobin  concentration,  A[Hb]totai,  and  changes  in  tumor  blood  volume,  AVT-biood, 
induced  by  hyperoxic  gas  intervention  in  rat  breast  tumors  using  NIRS  and  19F  MRS  of 
perflubron  (formerly  called  perfluorooctyl  bromide  or  PFOB),  respectively.  The  ratio  of  rat 
tumor  A[Hb]totai/AVT-biood  will  be  estimated  to  quantify  the  fractions  of  tumor  blood  volume 
sampled  by  NIRS  and  MRS  in  tumors. 


#6289 -$15.00  US 
(C)  2005  OSA 


Received  1 1  January  2005;  revised  25  February  2005;  accepted  28  February  2005 
7  March  2005  /  Vol.  13,  No.  5  /  OPTICS  EXPRESS  1726 


2.  Experimental  materials  and  methods 

2.1  Near  Infrared  spectroscopy  for  the  measurement  of  A[Hb]totai 

A  homodyne,  frequency-domain  NIRS  system  (NIM,  Philadelphia,  PA)  was  used  in  our 
experiments,  as  described  in  detail  previously  [10-12].  Briefly,  the  amplitude-modulated  light 
at  140  MHz  from  two  NIR  laser  diodes  at  758  nm  and  785  nm  was  projected  on  one  side  of 
the  tumor  through  a  delivery  fiber  bundle.  The  diffused  light  through  the  tumor  was  collected 
and  propagated  to  a  photomultiplier  tube  (PMT)  by  a  second  fiber  bundle.  The  signal  from  the 
PMT  was  demodulated  through  an  In-phase  and  Quadrature -phase  circuit,  and  the  amplitude 
and  phase  were  recorded,  as  shown  in  Fig.  1.  Based  on  modified  Beer-Lambert’s  law  [10,  20- 
21]  and  the  recorded  amplitude,  changes  in  oxygenated,  deoxygenated,  and  total  hemoglobin 
concentrations,  A[Hb02],  A[Hb]  and  A[Hb]totai,  due  to  respiratory  intervention  were  estimated, 
as  described  below.  In  this  NIRS  system,  only  amplitude  signal  was  used  for  the  calculation 
since  phase  signal  was  very  noisy. 


Fig.  1.  Schematic  diagram  of  experimental  setup  for  the  NIRS  system.  A  3  mm-diameter  fiber 
bundles  delivered  and  detected  the  laser  light  through  the  rat  tumor  in  transmittance  geometry.  PMT 
represents  a  photomultiplier  tube.  I/Q  is  an  in-phase  and  quadrature  phase  demodulator  for  retrieving 
amplitude  information. 


2.2  ,9F  MRS  of perflubron  for  the  measurement  of  tumor  blood  volume  Vj-biood 


An  Omega  CSI  4.7  T  superconducting  magnet  system  (Acustar™,  Braker  Instrument,  Inc., 
Fremont,  CA)  was  used  for  the  measurement  of  tumor  blood  volume.  An  emulsion  of 
perflubron  (Oxygent™,  Alliance  Pharmaceutical  Corp.,  San  Diego,  CA)  (2  ml)  was  infused 
into  tumor-bearing  rats  i.v.  as  a  blood  volume  indicator.  The  tumors  were  placed  within  a 
frequency- tunable  ( 1 H,  l9P),  single -turn,  solenoid  coil,  together  with  a  sealed  capillary 
containing  sodium  trifluoroacetate  (TFA),  which  was  used  as  an  external  standard  reference 
for  calibrating  tumor  blood  volume. 

Under  fully  relaxed  conditions,  integration  of  19F  signal  from  the  tumor  was  linearly 
proportional  to  the  total  number  of  19F  nuclear  spins  of  perflubron  in  the  tumor,  which,  in 
turn,  was  linearly  proportional  to  the  total  blood  volume  in  the  tumor,  assuming  that  the 
perflubron  emulsion  had  reached  an  equilibrium  state  with  the  blood  throughout  the  tumor.. 
To  ensure  full  relaxation  a  repetition  time  (TR)  of  30  seconds  was  used  [22].  The  spectral 
peaks  of  perflubron  were  integrated  in  the  data  post-processing.  Following  the  tumor 
measurement,  the  rat  was  removed,  leaving  the  reference  TFA  capillary  in  the  original 
position  within  the  RF  coil.  A  small  amount  of  blood  (0.5  ml)  was  drawn  from  the  rat  tail  vein 
and  placed  in  the  RF  coil  without  disturbing  the  reference  TFA  capillary.  Another  quantitative 
19F  spectrum  was  then  acquired.  Thus,  the  tumor  blood  volume  can  be  calculated  based  on  the 
following  Eq. : 


Vt  blood  —  Vs  blood  •  ( 


It  blood  x  ,1s  TFA 


-)(- 


)• 


XS  blood 


(1) 


where  Vt  blood  and  Vs  blood  were  the  tumor  blood  volume  and  blood  sample  volume  (ml), 
respectively,  It  blood  and  Is  blood  were  the  integrated  MRS  signal  of  19F  taken  from  the  rat 
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tumor  and  blood  sample,  respectively,  and  Ibiood  tfa  and  ITumor _tfa  were  the  integrated  19F 
MRS  signals  from  the  TFA  capillary  for  the  two  respective  measurements. 

2.3  Animal  model  and  protocols 

Mammary  adenocarcinomas  13762NF  were  implanted  in  skin  pedicles  [23]  on  the  forebacks 
of  female  Fisher  344  rats  (—150  g).  Once  the  tumors  reached  1~2  cm  diameter,  rats  were 
anesthetized  with  ketamine  hydrochloride  (100  mg/ml,  i.p.)  and  maintained  under  general 
gaseous  anesthesia  with  1.3%  isoflurane  in  air  (1  dmVmin).  Tumors  were  shaved  to  improve 
optical  contact  for  NIR  light  transmission.  Tumor  volume,  V,  was  estimated  using  an 
ellipsoid  approximation  as  V=(7t/6).a.b.c  from  the  three  orthogonal  diameters,  i.e.,  a,  b,  c, 
measured  with  a  caliper. 

In  this  study,  breast  tumor  bearing  Fisher  rats  were  challenged  with  carbogen  (95%  O2  + 
5%  C02)  using  a  sequence  of  air-carbogen-air.  The  A[Hb]totai  values  in  response  to  the 
carbogen  inhalation  were  monitored  by  NIRS.  For  the  rats  that  exhibited  a  significant  change 
in  [Hb] total  induced  by  carbogen  (n=6),  they  were  re -anesthetized  the  next  day  after  the  NIRS 
measurement  and  infused  with  2  ml  of  perflubron  emulsion  for  the  tumor  blood  volume 
determination.  The  19F  MRS.  experiment  was  performed  on  the  rat  30  min  after  the  emulsion 
injection  allowing  sufficient  time  for  the  perflubron  emulsion  to  reach  an  equilibrium  state 
within  the  tumor  blood  vasculature.  The  same  inhaled  gas  sequence  was  repeated  during  the 
MRS  measurements. 

3.  Algorithms  for  calculations  of  A[Hb]totai  and  fraction  of  sampled  blood  volume 

3.1  Calculations  for  tumor  A[Hb]tomi 

Based  on  Beer-Lamberf  s  law,  optical  attenuation,  presented  by  Optical  Density  (OD),  can  be 
expressed  as  a  function  of  concentration  of  chromospheres  (24-28), 

Optical  Density  (OD)=  Log(A(/A)  =ecl  (2) 

where  A0  and  A  are  light  intensities  of  the  incident  and  transmitted  light,  respectively,  £  is  the 
extinction  coefficient  of  chromophore,  c  is  the  concentration  of  chromophore,  and  /  is  the 
optical  path-length  through  the  measured  sample.  When  the  measured  sample  has  a  mixture 
of  chromophores,  e.g.,  oxygenated  and  deoxygenated  hemoglobin  in  tumor  tissue,  the 
changes  in  oxygenated  and  deoxygenated  hemoglobin  concentrations  (i.e.,  A[Hb02]  and 
A[Hb])  can  cause  a  change  in  optical  density,  which  can  be  presented  as  [26-28]: 

AOD(A)  =  logfAfAd  =  {em(A)A[Hb]  +  £Hb02(A)A[Hb02 1}  l ,  (3) 

where  AOD(A)  represents  a  change  in  optical  density  at  the  specific  wavelength,  A.  AB  and  Ar 
correspond  to  light  intensities  measured  under  the  baseline  and  transient  conditions.  £Hb(A) 
and  £/iho2( A)  are  extinction  coefficients  at  wavelength  A  for  molar  concentrations  of 
deoxygenated  hemoglobin  and  oxygenated  hemoglobin,  respectively.  By  employing  two 
wavelengths  in  Eq.  (3),  both  A[Hb02]  and  A[Hb]  can  be  determined  by  measuring  the  AOD 
values  at  the  two  specific  wavelengths,  provided  that  the  values  for  £nb(A)  and  £Hb02(A)  are 
known: 


A  [HbQ2] 


£Hb  (h  )AOD(A\ )  -  £p[b{A\  )AOD(A2 ) 

1  i£Hb(h  )£Hb02  Ul )  “  £Hb ( h  )£Hb02  (^2 )] 


(4) 


Aj- Hb j  _  £HbQ2  (h  )A OD(A\ )  -  £HbQ2(^\  )A OD(A2 ) 
1  i£Hb  Ul  )£Hb02  (h  )  -  £Hb(A 2  )£Hb02^\ )] 
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Note  that  in  principle,  /  represents  the  optical  path  length  between  the  source  and  detector. 
While  1  is  simply  the  physical  separation,  d,  between  the  source  and  detector  through  a  non¬ 
scattering  medium,  exact  quantification  of  /  for  an  intact  tissue  or  organ  is  complex  because 
of  light  scattering  in  tissue.  Since  1  is  in  proportion  to  the  separation,  d,  we  can  associate  /  to  d 
as  l=DPF*d ,  where  DPF  is  a  differential  path  length  factor  to  account  for  light  scattering 
[20].  It  has  been  well  accepted  [9,  10,  20]  that  together  with  DPF,  Eq.  (3)  can  be  treated  as 
modified  Beer-Lamberf  s  law;  and  consequently,  Eqs.  (4)  and  (5)  can  be  used  to  quantify 
changes  in  [Hb]  and  [FLbOJ  in  highly  scattering  media,  such  as  in  intact  tissue  or  organs. 

To  be  consistent  with  our  previous  work,  we  adopt  in  this  paper  the  8  values  published  by 
Zijlstra  et  al.  [29].  We  had  to  interpolate  the  8  values  at  the  two  wavelengths  employed  in  our 
study,  i.e.,  eHb(758  nm)  =  1.418,  eHb02(75  8  nm)  =  0.6372,  eHb(785  nm)  =  1.111,  and  eHb02(785 
nm)  =  0.766,  all  in  mVT'cnf Note  that  a  factor  of  4  has  been  multiplied  for  each  of  the  e’s  at 
the  respective  wavelengths  to  account  for  light  absorption  from  4  hemes  per  hemoglobin 
molecule  [30]  since  the  extinction  coefficients  published  in  the  field  of  biochemistry  were 
expressed  on  a  per  heme  basis  [26-30]. 

Because  of  the  interpolation  of  e  values,  necessary  calibration  is  needed  to  assure  the 
accuracy  of  A[Hb02 ],  A[F[b],  and  A[Flb]totai  calculation  as  given  in  Eqs.  (4)-(5).  We 
conducted  a  set  of  phantom  experiments,  as  previously  described  in  ref.  1 1 .  The  empirical 
calibration  has  resulted  in  two  correction  factors  (3)  and  ff.  where  |3i=1.103  and  (32=0.90  3  5,  as 
expressed  below: 


gffi>(785wm)  xA6>£)(758„m) -£Hb(75$nm)x  AOD(7S5nm) 

A[Hb(h  ]  =  -p - & - T , 

l[£Hb(7%5nm)£Hb02(75%nm)  ~  £Hb^^nm)£Hb02^J^nm)\ 

£Hb02(™5"m)  x  AO£)(758 _ eHhol (758 nm) x AOD(785/?w) 

_ _ P2 _ 

1  i£Hb  (75%nm)£Hb02  (785  nm)  -  £Rb (785 nm)£Hb02  (758 nm)] 


(6) 


(7) 


After  substituting  all  of  the  parameters  into  Eqs.  (6)  and  (7),  we  arrived  at  Eqs.  (8)  and  (9) 
to  quantify  changes  in  hemoglobin  concentration: 


A[HbO:] 


-2.658xAOZ)(758wm)  +  3.743xAC>Z)(785wm) 
dxDPF 


(8) 


A[Hb] 


2.238  x  AOD(758  nm )-1.683x  AOD(785  nm ) 
dxDPF 


(9) 


where  1  has  been  replaced  by  dxDPF.  Notice  that  the  difference  in  the  coefficients  between 
these  equations  and  those  published  earlier  [10-12]  results  from  the  factor  of  4  in  lvalues. 
A[Flb]totai  can  also  be  obtained  by  adding  Eqs.  (8)  and  (9): 

Aruui  Arui  n  7  4-  \ruK\  -0.42 x OD(758 nm)  +  2.06 x OD(785nm) 

A[Hb]tota,  =  A[HbO:]  +  A[Hb\  = - - - .  (10) 

d  xDPF 

In  principle,  units  of  A[Hb02 ■],  A[Hb],  and  A[Hb]totai  are  in  mM  with  £  in  cnf'.mM'1.  Since 
DPF  is  so  far  an  unknown  parameter  for  tumors,  we  include  it  within  the  unit  as  mM/DPF  for 
the  relative  measurements.  Thus,  Eq.  (10)  becomes  Eq.  (11),  and  it  permits  us  to  obtain 
A[F[b],otai  with  a  unit  of  mM/DPF  for  breast  tumors  under  carbogen  intervention. 
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(11) 


A[Hb]total  =  A[HbO:J  +  A[Hb ] 


— 0.42  xOD(758nm)  +  2.06xO£>(785/?w) 
d 


3.2  Calculations  for  the  fraction  of  vascular  blood  volume  sampled  by  NIRS 

The  determination  of  A[Hb]totai  can  further  lead  to  estimation  of  changes  in  total  hemoglobin 
content  within  a  rat  tumor,  A[CHb] tumor,  since  total  hemoglobin  content  of  a  tumor,  [CHt] tumor, 
is  directly  associated  with  hemoglobin  concentration  and  tumor  physical  volume.  Namely, 
[Cm] mmor  determined  by  NIRS  should  be  a  product  of  tumor  vascular  hemoglobin 
concentration,  [Hb]totai  in  mole/liter=T03  mM,  hemoglobin  molecular  weight,  M  in 
gram/mole,  and  the  tumor  physical  volume,  VT -physical  in  cm,  and  can  be  expressed  as 

/ C ffb] ’ tttmor_NIRS  [Hb] total  XK4 XC T-phvsical,  ( 1  2u) 

Or  Aj  C l ji] tun'.nr  MRS  A[Hb] total  X^XVi-physica!,  (12b) 

where  [Cm] tumor_NiRS  and  A[ CHb]tumor_MRs  are  in  gram,  and  A[CHb] tumor _nirs  is  induced  by 
carbogen  intervention. 

On  the  other  hand,  the  total  tumor  blood  volume,  VT-biood  in  liters,  can  be  quantified  from 
the  19F  MRS  measurement  with  and  without  carbogen  inhalation.  The  corresponding 
hemoglobin  content  determined  by  the  MRS  can  be  written  as: 

[Cffb] tumor _MRS  Vr-blood^R  ,  (13a) 

Cl  A  [ Cm]  tumor_MRS  AVT.bloodxK,  (13b) 

where  K  is  hemoglobin  concentration  in  gram/liter  for  liquid  blood,  VT-biood  can  be  determined 
using  Eq.  (1)  based  on  the  MRS  measurements,  and  AVT.biood  results  from  respiratory 
intervention. 

To  quantify  the  fraction  of  tumor  blood  volume  sampled  by  NIRS  and  MRS  in  tumors,  we 
introduce  a  variable,  y  to  associate  [CHb] ,Umor_NiRs  with  [CHb]  tumor _mrs  as: 

[C ffb] tumor_NIRS  Y [Cffb] tumor _MRS  •  (14a) 

Since  yis  expected  to  be  a  constant  for  a  given  tumor,  it  follows  that 

A  [ C ffb] tumor_NIRS  Y A[ C ffb] tumor_MRS  •  (14b) 

Then,  we  can  obtain  Eq.  (15)  by  taking  the  ratio  of  Eqs.  (13b)  and  (12b)  to  determine  the 
percentage  of  tumor  hemoglobin  content,  sampled  by  NIRS  and  by  MRS: 


y_  A[  CHb]  tumor  NIRS  *  jqqqx  _  M  xVT- physical  ^  A[Hb]total  ^  ]QQQ/0 

A[ C ] tumor  jrffis  K  AVT_biood 


(15) 


where  K  is  often  given  as  150  g/1  [30],  M  is  taken  as  68000  g/mol  [31],  A[Hb]totai  and  AVT-biood 
are  obtained  from  our  NIRS  and  MRS  experimental  measurements,  respectively,  and  VT-Phvsicai 
is  measured  from  the  tumor  physical  volume.  With  Eq.  (15),  we  are  able  to  compare  and 
estimate  the  fraction  of  hemoglobin  content  or  blood  volume  of  the  tumor,  y,  sampled  by  near 
infrared  spectroscopy  and  19F  magnetic  resonance  spectroscopy. 

4.  Experimental  results 

Changes  of  total  hemoglobin  concentration,  A[Hb],otai:  and  blood  volume,  VT_bi00d,  of  six 
selected  rat  tumors  were  monitored  by  NIRS  and  19F  MRS  of  perflubron  on  consecutive  days, 
respectively.  Figure  2  shows  the  time  course  profiles  of  A[Hb]totai  and  VT-biood  from  a 
representative  breast  tumor  (V=2.6  cm3)  with  respect  to  carbogen  intervention.  When  the 
inhaled  gas  was  switched  from  air  to  carbogen,  A[Hb]totai  increased  significantly  (yKO.0001) 
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from  a  baseline  value  of  0.0005  ±  0.0015  mM/DPF  to  a  maximal  value  of  0.020  ±  0.001 
mM/DPF  over  the  period  of  carbogen  intervention,  with  a  maximal  A[Hb],otai  of  0.0185 
mM/DPF.  After  the  gas  was  switched  back  to  air,  a  significant  drop  (p<  0.0001)  of  A[Hb]totai 
occurred,  followed  by  a  plateau  at  0.009  ±  0.001  mM/DPF.  A  similar  temporal  pattern  in 
response  to  carbogen  intervention  was  found  for  VT-bi0od  derived  from  the  MRS  measurement. 
Specifically,  Vr-biood  increased  significantly  (p<0.0001)  from  a  baseline  of  0.809  ±  0.004  cm3 
to  a  maximum  of  0.837  ±  0.004  cm3  during  carbogen  intervention,  having  a  maximal  AVr-biood 
of  0.028  cm3. 


Fig.  2.  Temporal  profile  of  A[Hb]Totai  and  Vr-biood  for  a  representative  rat  breast  tumor  (2.6  cm3) 
with  respect  to  carbogen  (CB)  intervention,  monitored  by  NIRS  and  19F  MRS  of  perflubron, 
respectively.  Error  bars  indicate  the  standard  deviation. 


By  substituting  A[Hb]totai,  AV T.blood,  K,  M,  and  the  corresponding  VT_Physica,=2.6  cm3  into 
Eq.  (15),  we  obtain 


7  = 


Hb\tumor  NIRS 
- - - X 1 00% 

[0 Hb^ tumor  MRS 


///>] tumor  NIRS 
- - - X 100% 

Hb] tumor  MRS 


68000  [g/mo/e]x2.6[cm3]  0.0185  [mM /DPF\  77.9% 

- — - - - - - -x - - - - - -100%  = - . 

150  [g/1]  0.028  [cm3]  DPF 


(16) 


Typical  values  of  DPF  for  the  brain  and  muscle  have  been  reported  in  the  range  of  4-6  and 
3-4,  respectively  [32-35],  independent  of  the  source-detector  separation  when  d  is  larger  than 
2.5cm  [36-37].  However,  the  DPF  values  for  tumors,  DPFlumor,  are  not  well  studied.  It  is  still 
reasonable  to  assume  DPFtumor  to  be  in  the  range  of  2-4  because  of  1)  higher  blood  content,  2) 
denser  blood  vasculature,  and  3)  the  finite  size  of  solid  tumors  in  this  study.  Thus,  the 
percentage  of  tumor  vascular  content  or  volume  sampled  by  NIRS  over  that  by  19F  MRS  is 
within  the  range  of  1 9%~39%,  based  on  the  result  of  Eq.  (16),  for  a  given  DPFtumor  of  4  and  2, 
respectively. 

Data  for  all  six  rat  breast  tumors  are  shown  in  Table  1.  The  percentage  of  tumor  vascular 
hemoglobin  content  sampled  by  NIRS  and  by  MRS  is  within  the  range  of  17-48%  if 
considering  a  DPF,umor  of  2,  and  8-24%  if  considering  a  DPFtumor  of  4.  The  mean  value  of  the 
percentage,  y,  is  about  15±6%  and  30±12%,  depending  on  the  DPF  values  of  actual  tumors. 
No  obvious  correlation  was  found  between  percentage  values  and  tumor  volume. 
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Table  1.  Percentage  of  tumor  vascular  blood  volume  sampled  by  NIRS  over  that  by  19F  MRS 


Tumor  physical 
volume  (cm3) 

A[Hb]totai  from 
NIRS  (mM/DPF) 

AVT_blood  from 

MRS  (cm3) 

Y  (%) 

with  DPF=2 

Y  (%) 

with  DPF=A 

1.6 

0.026 

0.056 

16.8 

8.4 

1.2 

0.062 

0.067 

25.2 

12.6 

1.9 

0.053 

0.081 

28.2 

14.1 

2.5 

0.038 

0.045 

47.9 

23.9 

2.6 

0.019 

0.028 

40.0 

20.0 

2.2 

0.022 

0.056 

19.6 

9.8 

Mean  ±  S.D. 

0.037  ±0.018 

0.056  ±0.018 

30  ±  12 

15  ±  6 

*S.D.  indicates  standard  deviation. 


5.  Discussion 

In  this  study,  changes  in  total  hemoglobin  concentration,  A[Hb]totai,  and  in  total  blood  volume, 
AV T-bioodj  of  rat  breast  tumors  in  response  to  carbogen  intervention  were  determined  by  NIRS 
and  by  19F  MRS  of  perflubron,  respectively.  The  ratio  of  A[  HbJtotai/AFr-iw  was  used  to 
calculate  the  fraction  of  vascular  hemoglobin  contents  sampled  by  NIRS  and  MRS.  Since  19F 
MRS  of  perflubron  samples  all  blood  volume  in  the  tumor  [18-19],  the  MRS  measurement 
may  be  considered  as  a  “gold  standard”  for  the  complete  total  blood  volume  within  the  tumor, 
provided  that  the  entire  rat  tumor  was  placed  within  the  RF  coil.  Then,  this  experimental 
approach  provides  evidence  that  the  vascular  blood  content  or  volume  sampled  by  the  NIR 
system  is  about  15-30%  of  the  total  blood  volume  (Table  1).  This  estimation  should  be 
applicable  not  only  to  tumor  tissues  but  also  to  other  kinds  of  tissues,  providing  us  with 
important  information  for  better  understanding  of  blood  vascular  volume  sampled  by  NIRS. 
Moreover,  the  consistent  trends  between  A[Hb]totai  and  AVT-biood  during  hyperoxic  gas 
intervention  (Fig.  2)  demonstrate  that  the  low  cost,  portable  NIRS  system  provides  a  non- 
invasive,  real-time,  monitoring  tool  to  detect  changes  in  tumor  blood  content. 

Is  it  reasonable  to  believe  that  NIRS  samples  only  1 5-30%  of  the  total  blood  content  or 
volume  in  tissue  vasculature?  It  has  been  reported  that  NIR  light  is  most  sensitive  to 
microvessel  densities,  such  as  arterioles,  capillaries,  and  venules  [13-14]  and  surface  areas  of 
microstructures  in  tissue  [15].  The  size  of  micro  blood  vessels  is  often  in  the  range  of  -100 
pm,  which  is  in  the  same  order  of  magnitude  as  the  average  scattering  length  of  NIR  light  in 
tissue  [15,  31].  So  NIR  light  is  scattered  more  frequently  by  the  small  blood  vessels  than  large 
vessels,  whereas  large  vessels  have  larger  chances  to  terminate  the  NIR  light  due  to  high 
absorption.  Thus,  multiple  scattering  from  micro  vasculature  permits  better  sensitivity  to  be 
sampled  by  NIR  light.  Moreover,  small  vessels  have  lower  hematocrit  compared  to  big 
vessels  because  of  reduction  of  hematocrit  with  the  decrease  of  vessel  diameters  [38].  Thus, 
based  on  our  experimental  measurements,  it  is  reasonable  to  conclude  that  the  NIR  system 
samples  15-30%  portion  of  total  hemoglobin  contents  in  tissue  vasculature,  mainly  sensing 
the  signals  in  microvessels. 

DPF  is  a  very  important  parameter  for  quantifying  NIR  measurements.  As  mentioned 
earlier,  many  investigators  have  conducted  research  on  the  estimation  of  optical  path  length  of 
the  brain  and  muscle  of  both  human  subjects  and  animal  models  [20,  32-37].  However,  the 
NIR  light  propagated  in  tumor  tissue  is  different  from  that  in  brain  and  muscle.  Due  to  tumor 
angiogenesis,  tumor  tissues  possess  more  blood  vessels  and  have  higher  light  absorption  than 
the  normal  brain  and  muscle.  Furthermore,  the  solid  tumors  in  the  study  have  a  finite  size, 
where  “photon  escape”  from  the  measured  volume  can  take  place  [39].  The  term  of  “photon 
escape”  has  been  used  to  refer  to  the  situation  where  the  sample  size  under  study  is  finite,  and 
a  large  amount  of  light  may  escape  from  the  sample  and  never  be  detected,  leading  to  a 
shorter  path  length  in  comparison  with  the  path  length  existing  in  a  large  volume  of  tissues. 
Steen  et  al  took  DPF  =  2.5  for  their  tumor  study  [40].  Therefore,  we  believe  that 
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DPFtumor=2~4  is  a  reasonable  range  used  in  our  estimation  of  hemoglobin  contents  for  breast 
tumors. 

The  errors  in  the  percentage  measurement  may  be  potentially  reduced  by  conducting 
measurements  of  NIRS  and  19F  MRS  simultaneously,  eliminating  experimental  variations 
between  the  two  kinds  of  measurements.  One  might  also  apply  'H  MRI  of  vascular  volume 
markers  such  as  super  paramagnetic  iron  oxide  particles  (SPIOs)  [41-42]  in  place  of  the  19F 
NMR  approach  used  here.  Since  agents,  such  as  Combidex,  are  in  clinical  use,  correlation 
studies  with  NIRS  could  allow  examination  of  percentage  of  blood  content/volume  sampled 
by  the  NIRS  system  in  human  tumors.  Also,  the  uncertainty  in  both  the  size  measurement  and 
volume  calculation  of  the  tumors  should  be  reduced  since  they  all  contribute  critically  to  the 
error  in  the  estimated  fraction  of  blood  vascular  contents. 

In  our  study,  not  all  tumors  showed  significant  changes  in  [Hb]totai  accompanying 
challenge  with  hyperoxic  gas.  This  is  mainly  due  to  both  heterogeneity  and  irregularity  of 
breast  tumors,  which  have  leaky,  broken,  and  dis-functioning  blood  vasculatures.  The 
inconsistency  of  changes  in  tumor  [Hb]total  induced  by  hyperoxic  gas  interventions  has  been 
observed  throughout  our  research  studies  over  several  years.  Thus,  we  believe  it  is  reasonable 
to  select  rat  breast  tumors,  which  show  significant  changes  in  [Hb]totai,  for  this  study  in  order 
to  obtain  estimated  fraction  of  vascular  blood  contents  sampled  by  NIRS. 

In  conclusion,  we  have  introduced  a  new  experimental  approach  to  estimate  the 
percentage  of  tumor  blood  content  or  volume  sampled  by  the  NIRS  system  by  combining  the 
measurements  of  NIRS  with  19F  MRS.  The  results  gave  rise  to  an  estimated  range  of  1 5 — 30%, 
suggesting  that  NIRS  does  not  sample  the  entire  blood  hemoglobin  content/volume  of  the 
tissue,  but  it  senses  more  toward  microvasculature.  The  study  helps  better  understand  blood 
vascular  volume  sampled  by  NIRS,  and  demonstrates  that  the  low  cost,  portable  NIRS  system 
can  be  a  reliable,  non-invasive,  real-time,  monitoring  tool  on  changes  in  tumor  blood 
contents. 
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Abstract:  We  demonstrate  the  ability  to  simultaneously  investigate  tumor  oxygen  dynamics  by 
Near-infrared  spectroscopy  and  Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for 
Dynamic  Oxygen  Mapping  in  breast  tumors,  with  close  correlations  between  the  two  techniques. 
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1.  Introduction 

Tumor  oxygenation  has  been  well  recognized  as  a  crucial  factor  governing  the  efficacy  of  radiotherapy  [1,2]. 
Accordingly,  accurate  evaluation  of  tumor  oxygenation  could  provide  a  better  understanding  of  tumor  response  to 
therapy.  Because  of  its  importance,  many  techniques  have  been  developed  for  monitoring  tumor  oxygenation  [3].  In 
recent  years,  near-infrared  spectroscopy  (N1RS)  has  been  developed  as  a  promising  non-invasive  and  real-time 
technique  to  quantify  tumor  vascular  oxygenation  [4,5].  Comparatively,  magnetic  resonance  imaging  (MRI)  offers 
alternative  noninvasive  approaches  to  monitor  tumor  oxygenation,  with  its  strength  of  high  spatial  resolution  but 
weakness  in  low  temporal  resolution.  By  using  Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for  Dynamic 
Oxygen  Mapping  (FREDOM),  MRI  can  provide  quantitative  pCL  readings  at  multiple  specific  locations 
simultaneously  with  millimeter  spatial  resolution.  It  is  of  interest  to  correlate  NIRS  and  FREDOM  in  order  to 
improve  and  complement  their  individual  strength  for  tumor  physiology  and  pathology  under  therapeutic 
interventions. 

2.  Materials  and  Methods 

2.1  Animal  Preparation  and  experimental  setup 

Ten  adult  female  Fisher  344  rats  were  planted  with  mammary  adenocarcinomas  13762NF  in  surgically  prepared 
foreback  pedicles.  When  the  tumors  reached  ~1  cm  in  diameter,  the  rats  were  anesthetized  with  ketamine 
hydrochloride  i.p.  (1.5ml;  lOOmg/ml;  Aveco,  Fort  Dodge,  1A)  and  were  maintained  under  general  gaseous 
anesthesia  (air  and  1%  isoflurane;  Baxter  International  Inc.,  Deerfield,  1L).  Tumor  hair  was  trimmed  for  the  ease  of 
optical  contact  for  transmitting  NIR  light  and  FIFB  injection.  The  rats  were  placed  in  the  magnet  on  their  side,  and 
body  temperature  was  maintained  using  a  warm  water  blanket.  A  Flelmholtz  coil  with  a  diameter  of  3  cm  was 
designed  to  facilitate  simultaneous  MRI  and  NIRS  measurements.  We  measured  tumor  oxygenation  parameters 
simultaneously  by  FREDOM  and  NIRS  during  respiratory  challenge  with  oxygen  or  carbogen.  The  response  was 
characterized  with  both  the  signal  magnitude  and  rate  of  oxygenation. 

2.2  NIRS  for  measuring  changes  in  tumor  vascular  oxygenation  (AHb02) 

We  used  a  homodyne,  frequency-domain  system  (NIM,  Philadelphia,  PA)  to  monitor  global  changes  of  oxy-  and 
deoxy-hemoglobin  concentration  of  tumors  in  response  to  different  gas  interventions.  Light  from  785  nm  and  758 
nm  laser  diodes  was  modulated  at  140  MHz  and  was  coupled  into  a  22 -foot,  bifurcated  fiber  bundle.  The  probe  tips 
were  made  with  a  non-metallic  material  and  in  good  contact  on  the  surface  of  tumor.  After  being  absorbed  and 
scattered  in  the  tumor  tissue,  the  transmitted  light  was  collected  by  another  fiber  bundle  and  amplified  by  a 
photomultiplier  tube  (PMT).  An  In-phase  and  Quadrature  (IQ)  demodulator  chip  was  used  to  demodulate  the 
measured  signal  from  PMT.  More  details  about  this  system  have  been  described  previously  [4,6]. 

2.3  Mathematical  model  for  the  blood  oxygenation  dynamics  of  tumors 

In  our  previous  study,  we  derived  a  simplified  model  to  relate  the  time  constants  to  tumor  blood  perfusion  rates  [4]. 
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—ft  —fit  —ft  — fit 

A[Hb02(t)rcula,ure  =  nH.il- e^'  ]  +  r2H0[ l-e»  ]  =  A,[  1-e*  ]  + A2[ l-e*  ].  (1) 

where  H0  is  the  arterial  oxygenation  input,//  and  yi  are  the  blood  perfusion  rate  and  the  vasculature  coefficient  in  a 
well-perfused  region,  and  f2  and  yi  are  the  blood  perfusion  rate  and  the  vasculature  coefficient  in  a  poorly  perfused 
region.  Also,  Ai=yiH0,  A2=  y2H0,  and  the  two  time  constants  in  eq.  (1)  are  equal  to  Ti=yi  If)  and  x2=y2/ /).  Since  A], 
A2  and  the  two  time  constants  can  be  determined  from  our  NFR  dynamic  measurements,  we  arrive  at  the  ratios  for 
the  two  vasculature  coefficients  and  the  two  blood  perfusion  rates,: 

n_  =  A_  f_  =  a2  (2) 

Yl  A2  ’  fl  Tl /r2  ’ 

Therefore,  by  studying  tumor  oxygenation  dynamics  with  NFRS  and  obtaining  the  time  constants  together  with  their 
amplitudes,  we  are  able  to  gain  important  information  on  regional  blood  perfusion  and  vascular  structures  of  the 
tumor  within  the  measured  volume. 

2.4  FREDOM  for  measuring  oxygen  tension  of  tumors  (p02) 

Simultaneous  MR1  experiments  were  performed  on  an  Omega  CSI  4.7T  horizontal  bore  system  equipped  with 
actively  shielded  gradients  (Acustar™,  Bruker  Instrument,  Inc.  Fremont,  CA).  Before  the  experiment,  HFB  (50  pi, 
99.9%,  Aldrich  Chemical  Co.,  St.  Louis,  MO)  was  administered  along  two  or  three  tracks  in  central  and  peripheral 
regions  of  the  tumors  in  a  single  plane  (transverse  to  the  rat’s  tumor,  and  in  the  region  of  N1R  photon  pathway) 
using  a  Hamilton  syringe  with  a  fine  sharp  32  G  needle.  The  tumor  was  placed  inside  a  Helmholtz  coil  designed  to 
facilitate  simultaneous  MR1  and  N1RS  measurements.  Tl-weighted  images  were  acquired  with  TR=200  ms,  TE=10 
ms,  and  128x128  points  over  a  40x40  mm  area,  giving  312x312  pm  resolution.  Following  'H  MR  imaging,  the 
corresponding  19F  MR  imaging  (188.252  MHz)  was  performed  to  show  the  distribution  of  HFB  in  the  tumor. 
FREDOM  used  an  MBEST  spin  echo  planar  imaging  sequence  with  pulse  burst  saturation  recovery  (PBSR)  signal 
preparation,  a  variable  recovery  time  and  spin  echo  acquisition.  The  Alternated  R1  Delays  with  Variable 
Acquisitions  to  Reduce  Clearance  effects  (ARDVARC)  protocol  was  applied  to  improve  the  precision  of  R1 
measurements  [7].  R1  was  calculated  on  a  pixel-by-pixel  basis  using  nonlinear  least-squares  data  fitting  by  the 
Gauss-Newton  method.  Readings  for  p02  were  estimated  using  the  equation  p02  =  (Rl-0.0835)/0. 001876,  as 
described  in  detail  previously  [8].  Three  baseline  p02  maps  were  generated  over  24  mins  for  all  tumors,  and  then  the 
rats  were  exposed  to  inhaled  gas  interventions.  Five  FREDOM  readings  were  taken  during  each  gas  breathing 
period.  A  total  of  23  maps  were  obtained  during  a  complete  five-phase  imaging  session  over  about  2.5  hours. 

3.  Results 


air  oxygen  air  carbogen  air 


Fig.  1  Time  course  of  AHbC>2  and  pC>2  in  response  to  respiratory  challenge  for  tumor  No.  1,  simultaneous  measured  by 
NIRS  (curve)  and  FREDOM  (symbols). 

Figure  1  shows  simultaneous  readings  of  mean  p02  and  relative  changes  of  A[Hb02]  of  a  breast  tumor  in  response  to 
respiratory  challenge  measured  by  FREDOM  and  NIRS.  This  tumor  showed  a  rapid  response  to  the  hyperoxic  gas 
interventions  with  a  significant  increase  in  A[Hb02].  The  p02  readings  increased  significantly,  but  more  slowly. 
These  changes  can  be  characterized  by  their  magnitude  of  response  and  kinetics.  Namely,  eq.(l)  can  be  used  to  fit 
A[Hb02],  whereas  the  p02  data  can  be  quantified  by  a  rising  rate,  which  is  the  slope  of  best-linear-fit  in  the  p02  time 
course  during  oxygen  or  carbogen  intervention.  Compared  with  Aj/ti  and  A2/t2  by  fitting  the  AHb02  readings,  the 
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PO2  rising  rate  has  a  significant  positive  correlation  with  Ai/ij  in  the  well-perfused  region,  but  no  correlation  with 
that  in  the  poorly-perfused  region  within  the  tumor,  as  shown  in  Figure  2. 
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Fig.  2  A)  pOj  rising  rate  vs  Ai/ii  with  a  positive  correlation  between  the  two  variables  (R>0. 7, p<0. 02).  B)  pOi  rising  rate  vs 
A2/X2  with  no  correlation  between  them. 


Furthermore,  Figure  3  shows  that  a  good  correlation  existed  between  the  tumor  perfusion  rate  ratio  (fi/f2)  and  tumor 
p02,  where  the  latter  values  are  the  average  oxygen  tensions  determined  in  the  non-hypoxic  regions  of  tumor  during 
the  hyperoxic  interventions. 


Fig.3  Tumor  pO,  vs  perfusion  rate  ratio  (///jG)  (R>0.9,  p<0.001).  The  pCF  readings  are  the  average  of  maximal  p02  values  in 
selected  voxels,  which  have  pCb  increased  to  more  than  10  torr  during  oxygen  (solid  points)  or  carbogen  intervention  (open 
points). 

4.  Discussion 

This  study  shows  that  the  two  techniques,  i.e.,  NIRS  and  MRI,  are  consistent  and  complementary  with  one  another  for  tumor 
oximetry.  It  also  demonstrates  the  feasibility  to  correlate  tumor  vascular  oxygenation  with  tumor  tissue  pCte,  indicating  the 
potential  applications  of  the  NIRS  technique  for  clinical  cancer  practice  and  providing  valuable  insights  into  tumor 
patho-physiology  and  response  to  therapy. 
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Abstract:  Dynamic  changes  of  oxygenated  and  deoxygenated  hemoglobin  concentrations  in  response  to 
hyperoxic  gas  interventions  on  rat  breast  tumors  were  simultaneously  investigated  by  near  infrared 
spectroscopy  and  BOLD  (blood  oxygenation  level  dependent)  contrast  MR  imaging. 
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1.  Introduction 

Hyperoxic  gases  have  been  investigated  as  potential  contrast  agents  in  magnetic  resonance  (MR)  imaging. 
Manipulation  of  oxygen  (CL)  content  in  the  inhaled  gas  can  provide  qualitative  and  quantitative  information  on 
regional  oxygenation.  During  hypoxia,  the  concentration  of  deoxyhemoglobin  ([Hb] )  is  increased,  causing  decreased 
signal  intensity  (SI)  on  T2*-weighted  MR  images1'2.  Oxygenation,  on  the  other  hand,  increases  local  T2*,  resulting 
in  increased  SI3.  Inhalation  of  02  or  carbogen  (95%02,  5%C02)  in  tumor-bearing  rats  alters  blood  deoxy¬ 
hemoglobin  concentration,  resulting  in  SI  changes  on  T2*  weighted  images  of  the  measured  tumors.  Therefore, 
blood  oxygenation  level  dependent  (BOLD)  contrast  MR  imaging  has  been  utilized  to  detect  hemodynamic  changes 
in  the  brain  and  tumor  tissues4'5.  It  has  also  used  to  evaluate  the  physiologic  parameters  of  tissue  metabolic 
activities,  such  as  cerebral  blood  flow  (CBF)  and  02  consumption.  Results  of  BOLD  MR  imaging  for  CBF  and  02 
consumption  are  compatible  to  those  obtained  with  animal  models  and  PET6  7 

While  MR  imaging  provides  a  noninvasive  in  vivo  imaging  modality  for  mapping  tumor  hemodynamic 
changes,  the  signal  from  BOLD  contrast  MRI  is  also  influenced  by  variations  in  blood  volume  (,  hemoglobin 
concentration)  and  blood  flow.  This  confusion  may  limit  its  wide  clinic  application.  Near  Infrared  spectroscopy 
(NIRS)  provides  an  alternative  noninvasive,  real  time,  monitoring  means  for  tumor  vascular  hemoglobin 
concentration.  Moreover,  the  signals  from  NIRS  are  independent  of  blood  flow.  Therefore,  simultaneous 
measurements  by  NIRS  and  BOLD  MRI  will  allow  us  to  correlate  the  signals  from  both  non-invasive  modalities 
and  to  better  interpret  the  signals  of  BOLD  MRI.  In  this  way,  the  ability  of  understanding  BOLD  signal  better  will 
increase  the  potential  of  BOLD  MRI  to  be  used  for  clinical  applications,  particularly  for  cancer  research. 

In  this  study,  we  performed  simultaneous  measurements  of  dynamic  oxygenated  and  deoxygenated 
hemoglobin  concentrations,  A[Hb02]  and  A[Hb],  using  NIRS  and  BOLD  contrast  MRI,  in  response  to  carbogen 
inhalation  on  rat  breast  tumor  vasculature.  The  average  signals  from  the  BOLD  MR  imaging  were  correlated  well 
with  the  tumor  global  A[Hb02]  signals  obtained  from  the  NIRS  and  BOLD  MR  imaging  revealed  the  heterogeneity 
of  tumor  hypoxia. 

2.  Material  and  Methods 

2.1  Near  Infrared  Spectroscopy  for  global  measurements  of  tumor  A[Hb02]  and  A[Hb]totai 

The  changes  of  oxygenated,  deoxygenated  hemoglobin  concentrations  and  total  hemoglobin  in  tumor  tissue, 
A[HbO?],  A[Hb]  and  A[Hb]totai,  caused  by  external  interventions  were  determined  by  a  homodyne  frequency-domain 
(NIRS  system  (NIM,  Philadelphia,  PA)8'9'10.  Briefly,  the  RF-modulated  light  at  the  wavelengths  of  758  nm  and  785 
nm  was  projected  on  the  tumor  surface  through  a  delivery  fiber  bundle.  The  light  diffusing  through  the  rat  tumor 
tissue  was  collected  and  propagated  to  a  photomultiplier  tube  (PMT)  through  another  fiber  bundle.  The  signal  from 
the  PMT  was  demodulated  through  an  In-phase  and  Quadrature -phase  circuit,  and  the  amplitude  and  phase  were 
recorded  in  real  time.  Based  on  modified  Beer-Lambert’s  law  and  recorded  amplitudes,  the  values  of  A[Hb02], 
A[Hb]  and  A[Hb]tota|  due  to  respiratory  intervention  were  derived  previously9. 


2.2  Magnetic  resonance  imaging  system 

BOLD  contrast  was  measured  by  using  a  Varian  IVOVA  4.7  T  super-conducting  magnet  system  single  turn  solenoid 
RF  coil.  The  coil  was  tuned  at  proton  resonance  frequency  (200.13  MHz  at  4.7  T).  The  paramagnetic 
deoxyhemoglobin  creates  local  magnetic  field  variations  in  the  proximity  of  blood  vessels,  producing  additional 
phase  dispersion  of  water  proton  signals  in  the  surrounding  tissue.  Gradient  recalled  echo  (GRE)  'H  MR  imaging 
sequences  are  sensitive  to  these  local  susceptibility  variations  near  the  vessels;  hence,  an  increase  in  concentration  of 
deoxyhemoglobin,  due  to  oxygen  desaturation  or  blood  flow  changes,  leads  to  a  signal  attenuation  in  gradient  echo 
or  T2  (apparent  spin-spin  relaxation  time)-weighted  image.  The  parameters  for  data  acquisition  of  T2  -  weighted 
image  were:  TR  (Recycle  time)=300ms,  TE  (Echo  time)=30ms,  Flip  angle=20  °. 

2.3  Animal  model  and  protocol 

Mammary  adenocarcinomas  13762NF  were  implanted  in  skin  pedicles  on  the  forebacks  of  female  Fisher  344  rats 
(—150  g).  Once  the  tumors  reached  1~2  cm  diameter,  rats  were  anesthetized  with  ketamine  hydrochloride  (100 
mg/ml,  i.p.)  and  maintained  under  general  gaseous  anesthesia  with  1.3%  isoflurane  in  air  (1  dm  Vmin).  Tumors  were 
shaved  for  better  optical  contact  and  placed  in  the  frequency-tunable  (*H/19F),  single-turn  solenoid  RF  coil.  The 
delivering  and  detecting  fiber  bundles  of  the  N1RS  were  touched  firmly  to  the  surface  of  the  shaved  tumors  (without 
compression  of  tumor)  for  the  measurement  of  A[Hb],otai  and  A[Hb02].  BOLD  MR1  was  performed  simultaneously 
with  N1RS  during  hyperoxic  gas  interventions.  Specifically,  two  experiments  were  setup  in  each  of  the  tumor 
measurements:  the  BOLD  MRI,  was  first  performed  on  a  single  thick  slice  equivalent  to  the  volume  interrogated  by 
N1R  and  then  from  three  thinner  slices  (2  mm)  were  selected  to  study  the  heterogeneity  of  tumor  oxygenation. 

3.  Results 

3.1  Consistency  of  tumor  vascular  oxygenations  obtained  from  N1RS  and  BOLD  MRI 


Figure  1  Simultaneous  measurements  with  BOLD  MRI  and  NIRS  of  a  rat  breast  tumor  (No.  2,  4.8  cm  3)  for  its 
vascular  oxygenation  status  with  respect  to  carbogen  intervention,  (a)  Average  intensity  of  BOLD  MR  images 
within  the  tumor  region,  (b)  The  changes  in  oxygenated  hemoglobin  concentration,  A[Hb02],  measured  by  NIRS. 


The  whole  tumor  was  first  observed  as  one  thick  slice  for  the  simultaneous  measurements.  Figure  1  shows  the  time 
course  profiles  of  average  SI  from  BOLD  MR  images  (a)  and  A[Hb02]  obtained  from  NIRS  (b).  Analysis  of  ROl 
(region  of  interest)  on  whole  tumors  gave  rise  to  the  mean  Sis  of  each  MR  image  (Fig.  1(a)).  During  the  first  20 
minutes,  when  the  rat  breathed  air  as  baseline,  the  mean  intensity  of  BOLD  signals  is  about  0.0162±0.0005. 
Carbogen  intervention  after  the  baseline  breathing  produced  a  significant  increase  (p<0.001)  in  BOLD  MRI 
intensity,  to  about  0.0189±0.0008.  The  signal  decreased  significantly  to  0.011  ±0.00 1  when  the  inhaled  gas  was 
switched  back  to  air  (p<0.001).  A  similar  tendency  in  A[Hb02]  was  observed  by  NIRS  (Fig.  1(b)),  indicating  that 
both  BOLD  MRI  and  NIRS  provide  consistent  results  for  oxygenation  dynamics  in  tumor  vasculature. 


3.2  Heterogeneity  of  tumor  oxygenation  dynamics 


Figure  2  (a)  Relative  signal  intensity  of  BOLD  MR  images  at  different  tumor  regions  in  response  to  carbogen 

intervention  for  a  representative  breast  tumor  (3.2  em  3).  (b)  Typical  cross-sectional  intensity  images  of  BOLD 
MRI  under  three  different  gas  inhalations  for  the  same  breast  tumor. 

The  heterogeneity  of  tumor  oxygenation  with  respect  to  carbogen  intervention  for  a  representative  breast  tumor  (3.2 
cm* 2 3 4 5 6 7 8 9 10)  was  demonstrated  in  Fig.2.  Several  locations  exhibit  increased  signals,  whereas  others  remain  unchanged, 
when  the  breathing  gas  was  switched  to  carbogen  and  then  back  to  air. 

4.  Summary 

We  have  demonstrated  the  ability  to  correlate  NIR  and  1H  MRI  investigations  simultaneously.  This 
preliminary  study  demonstrates  that  the  global  signals  of  breast  tumors  obtained  by  N1RS  and  BOLD  MRI  in 
response  to  respiratory  intervention  are  in  good  agreement.  These  two  techniques  may  be  complementary  to  one 
another  for  in  vivo  measurements  of  tumor  oxygenation  dynamics.  MRI  can  be  applied  to  deep  tumors  with  special 
pulse  sequences,  whereas  NIRS  can  be  used  to  monitor  acute  responses  of  tumors  to  different  therapeutic 
interventions.  The  correlation  of  NIRS  and  BOLD  MRI  techniques  could  provide  a  new  prognostic  tool  for 
assessing  efficacy  of  tumor  therapies  and  could  help  us  better  interpret  BOLD  MRI  signals.  Furthermore,  the 
heterogeneity  of  tumor  vascular  oxygenation  dynamics  to  external  interventions  can  be  investigated  directly  by 
BOLD  MRI,  while  the  development  of  NIR  imaging  for  tumor  oximetry  has  also  great  potential  to  reveal  tumor 
vascular  heterogeneities. 
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Abstract:  Correlation  between  total  hemoglobin  concentration  and  blood  volume  in  solid  mammary 
tumors  was  explored  using  both  N1RS  and  19F  MRS  of  PFOB,  respectively,  while  studying  dynamic 
changes  in  tumor  vascular  volume  induced  by  a  vasoactive  agent,  carbogen.  A  close  correlation 
between  the  results  obtained  from  the  two  methods  indicates  the  potential  utility  of  N1RS  for  non- 
invasive  in  vivo  quantification  of  tumor  blood  volume. 
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1.  Introduction 

Tumor  blood  volume  is  an  important  physiological  parameter  that  affects  not  only  tumor  growth  and 
survival,  but  also  efficacy  of  cancer  treatments  that  alter  oxygen  delivery  and  chemotherapeutic  agents  to  tumor 
cells.  The  ability  to  measure  tumor  blood  volume  may  provide  a  valuable  diagnostic  tool  for  cancer  treatment 
planning.  Longitudinal  measurements  could  reveal  insight  into  tumor  development  and  perhaps  provide  a  prognostic 
indicator  for  tumor  therapy.  For  many  years,  the  most  commonly  used  methods  of  measuring  tumor  blood  volume 
involved  the  use  of  radioactive  labeled  tracers  and  also  required  the  excision  of  the  tumor1,2. 

Recently,  several  new  techniques  have  been  developed  for  in  vivo  measurements  associated  with  tumor 
blood  volume,  including  near  infrared  spectroscopy  (N1RS)3,  position  emission  tomography  (PET),  and  19F  NMR  of 
perfluorocarbons  (PFCs)4.  Since  N1RS  is  entirely  noninvasive,  inexpensive,  portable,  and  real-time,  it  would  be  a 
useful  and  desirable  technique  in  cancer  research  and  practice  if  direct  readouts  of  tumor  blood  volume  in  vivo  could 
be  obtained.  Flowever,  the  current  N1RS  can  measure  quantitative  values  of  only  total  hemoglobin  concentration, 
[Hbtotai],  and  oxygenated  hemoglobin  concentration,  [  HbCF].  Few  studies  have  reported  on  how  to  quantitatively 
convert  or  correlate  tumor  Flbtotai  to  tumor  blood  volumes.  In  this  paper,  we  report  our  preliminary  investigation  on 
the  correlation  between  the  Hbtotai  and  blood  volume  in  solid  mammary  tumors  using  both  N1RS  and  19F  MRS  of 
PFOB,  respectively,  while  a  vasoactive  agent,  carbogen,  was  employed  to  induce  dynamic  changes  in  tumor 
vascular  volume.  A  close  correlation  between  the  results  obtained  from  the  two  methods  supports  the  ability  of 
N1RS  for  in  vivo  quantification  of  tumor  blood  volumes.  Also,  the  vascular  density  of  breast  tumors  was  obtained  as 
a  function  of  the  tumor  volume. 

N1RS  is  based  on  the  endogenous  chromophores,  i.e.,  oxygenated  and  deoxygenated  hemoglobin.  The 
absorption  of  light  by  these  intrinsic  chromophores  can  be  altered  to  determine  important  physiological  indicators, 
such  as  Hbtotai  changes.  Others  have  shown  that  perfluorocarbon  emulsions  provide  a  valid  indicator  of  tumor 
vascular  volume  for  some  time  after  infusion.  19F  NMR  can  be  applied  to  measure  signal  intensity,  which  is  related 
to  blood  concentration.  In  this  study,  a  perfluorooctylbromide  (PFOB)  (CgF17Br)  emulsion  was  used  as  a  tumor 
blood  volume  indicator  based  on  the  earlier  work  of  Baldwin  et  at.  PFOB  is  extremely  inert  chemically,  and 
biochemically,  and  ideally  suited  for  in  vivo  biological  studies. 

2.  Materials  and  Methods 
2.1  Animal  Preparation 

Mammary  adenocarcinomas  13762NF  were  implanted  in  skin  pedicles  on  the  forebacks  of  female  Fisher  344  rats 
(~250  g).  Once  the  tumors  reached  ~  2  cm  diameter,  tumor  [Hb]Totai  and  [Hb02]  were  measured  by  dual-wavelengths 
N1RS  in  anesthetized  rats  (isoflurane/air).  On  the  next  day,  the  rats  were  re-anesthetized  and  intravenously  infused 
with  2  ml  of  a  PFOB  emulsion  by  tail  vein  injection.  The  19F  MRS  measurements  were  performed  30  min  later  to 


give  sufficient  time  for  the  PFOB  emulsion  to  reach  an  equilibrium  state  within  the  blood  stream.  To  study  the 
dynamic  change  of  blood  volume,  the  inhaled  gas  was  alternated  in  the  order  of  air  ->  carbogen  ->  air  ->  carbogen  - 
>  air. 

2.2  Near  Infrared  Spectroscopy 

[Hb02]  and  [Hbtota|]  can  be  determined  with  multi-wavelength  NIR  spectroscopy  based  on  their  absorption 
properties.  The  instrument  setup  was  described  elsewhere3,  but  the  algorithms  used  in  this  paper  to  obtain  A[  Hbtotai] 
were  modified  in  the  light  of  in  vitro  laboratory  calibration.  The  equation  for  A[Hbtotai]  is  given  below: 

A[Hb  totaJ  =[-2.12*Ln(AB/AT)758  +  8.25*Ln(AB/AT)785]  /  L,  (1) 

where  AB  =  baseline  amplitude;  AT  =  transient  amplitude;  L  =  optical  pathlength  between  source/detector;  the 
constants  given  in  the  equations  were  computed  with  the  extinction  coefficients  for  oxy-  and  deoxy-hemoglobin  at 
the  two  wavelengths  used.  In  principle,  the  unit  of  A[Hb]totai  is  mM,  and  L  is  equal  to  the  source-detector  separation 
multiplied  by  a  Differential  Pathlength  Factor  (DPF)5.  Little  is  known  about  the  DPF  for  tumors,  though  a  DPF 
value  of  2.5  has  been  used  by  others6.  Since  our  focus  is  on  dynamic  changes  of  tumor  [Hbtotai]  with  respect  to 
carbogen  intervention,  we  include  the  DPF  in  the  unit,  i.e.,  modifying  eq.  (1)  as  follows: 

A[HbtotaJ  =  [-2. 12*Ln(AB/AT)758  +  8.25*Ln(AB/AT)785]  /  d,  (2) 

where  d  is  the  direct  source-detector  separation  in  cm,  and  the  unit  of  A[Hbtotai]  in  eq.  (2)  is  mM/DPF. 

2.3  19 F  MRS  of  PFOB 

The  rats  were  placed  on  their  sides  in  an  animal  bed  and  the  body  temperature  was  maintained  at  37°  C.  A  frequency 
tunable  ( 1 H  l9F )  single  turn  solenoid  coil  was  placed  around  the  tumors  and  a  sealed  capillary  containing  sodium 
trifluoroacetate  (TFA)  used  as  an  external  standard  for  the  purpose  of  quantifying  tumor  blood  volume.  The  rats 
were  positioned  in  the  magnet  (Omega  CSI  4.7  T,  Bruker  Instrument,  Inc.,  Fremont,  CA)  with  the  tumors  at  the 
isocenter  of  the  magnet.  Spectra  were  acquired  with  a  long  repetition  time  (TR  =  30  s)  to  ensure  that  changes  in  T1 
due  to  variable  oxygenation  would  not  interfere  with  volume  measurements  (ref  Shukla)7.  The  rats  were  removed 
from  the  magnet  and  the  RF  coil,  with  the  reference  TFA  capillary  left  in  the  original  position.  0.5  ml  blood  sample 
drawn  by  tail  vein  was  put  into  the  RF  coil  without  disturbing  the  reference  TFA  capillary.  The  RF  coil  was  then 
placed  back  into  the  magnet  in  the  same  position  as  with  the  rats.  A  quantitative  19F  spectrum  was  acquired.  Given  a 
set  of  data  acquisition  parameters,  the  integration  of  19F  signal  from  the  tumor  was  linearly  proportional  to  the  total 
number  of  19F  nuclear  spins  of  PFOB  in  the  tumor,  which,  in  turn,  was  linearly  proportional  to  the  total  blood 
volume  in  the  tumor,  assuming  that  the  PFOB  emulsion  had  reached  an  equilibrium  state  with  the  blood  throughout 
the  tumor. 

The  tumor  blood  volume  is  given  as  vx  blood  =  Vs  blood  •  (lT-blood) .  (Is-TFA  ) ,  (2) 

Is_blood  It_TFA 

where  Is  blood  and  It  blood  were  the  integration  of  NMR  signal  from  19F  in  the  blood  sample  and  in  the  tumor,  Vs  blood 
and  VT  biood  were  the  blood  sample  volume  and  total  tumor  blood  volume  (ml),  respectively,  flood  tfa  and  humor  tfa 
were  the  integration  of  NMR  signal  from  19F  in  the  reference  TFA  capillary  acquired  with  the  blood  sample  and  the 
tumor. 


3.  Results  and  Discussion 
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Fig  1.  (A)  Time  course  profiles  of  tumor  vascular  A[HbO:]  and  A[Hb]Totai  determined  by  NIRS; 

(B)  tumor  blood  volume  Vmiojd  measured  by  19F  MRS  of  PFOB  for  breast  tumor  No.  1  (4.2  cm3). 


The  rat  with  mammary  tumor  first  breathed  air,  and  the  values  of  [Hb,ota|]  and  [Hb02]  in  the  tumor  were 
measured  for  baseline.  The  inhaled  gas  was  alternated  from  air  to  carbogen,  and  the  dynamic  changes  in  [Hb,ota|]  and 
blood  volume  of  the  tumor  were  monitored  by  N1RS  and  19F  MRS  of  PFOB,  respectively.  Figure  1(A)  shows  the 
time  course  profiles  of  A[Hb02]  and  A[HbTotai]  in  response  to  carbogen  intervention  measured  by  N1RS  for  a 
representative  breast  tumor  No.  1  (4.2  cm1 2 3 4 5 6).  Figure  1(B)  shows  the  tumor  blood  volume,  VT  b|00d,  as  a  function  of 
time  as  measured  by  19F  MRS  of  PFOB  for  the  same  breast  tumor.  In  spite  of  fluctuation,  the  overall  trend  of  VxBiood 
was  similar  to  that  of  A[FlbTota|].  The  correlation  between  the  A[HbTotai]  measured  by  N1RS  and  VT  blood  measured  by 
19F  MRS  of  PFOB  is  shown  in  Figure  2,  which  shows  reasonable  consistency  between  the  results  obtained  from  the 
two  methods. 


Tumor  Blood  Volume  (cm3) 


Fig.  2  Correlation  between  mean  A[Hb]x0tai  measured  by  NIRS 
and  mean  Vx  biood  measured  by  19F  MRS  of  PFOB  for  each 
gas  switch  for  breast  tumor  No.  1  (4.2  cm3). 


Fig.  3.  Relationship  between  tumor  vascular  density 
versus  tumor  volume  for  a  group  of  four  breast  tumors. 


Furthermore,  a  number  of  breast  tumors  were  investigated  with  19F  MRS  of  PFOB,  and  the  results  are 
compared  and  shown  in  Figure  3.  The  tumor  physical  volumes  were  calculated  from  the  three  dimensions  of  the 
tumor.  We  can  calculate  the  vascular  density  of  the  tumors  as  the  ratios  of  the  tumor  blood  volume  (measured  with 
air  inhalation)  to  the  tumor  physical  volume.  The  relationship  between  the  tumor  vascular  density  and  the  tumor 
physical  volume  is  plotted  in  Figure  3.  The  strong  negative  linear  relationship  suggested  that  tumors  became 
increasingly  poorly  perfused,  and  thus,  poorly  oxygenated  as  they  grew  in  size. 

4.  Summary 

In  summary  with  this  preliminary  study,  the  agreement  between  the  results  obtained  by  19F  MRS  of  PFOB 
and  NIRS  suggests  that  NIRS  is  a  valid  technique  for  in  vivo  measurement  of  changes  in  [Hbiotai]  and  tumor 
vascular  volume.  I9F  MRS  of  PFOB  can  be  applied  to  deep  tumors  with  special  pulse  sequences  and  NIRS  can  be 
used  to  monitor  the  acute  response  of  tumors  to  drug  intervention.  The  correlation  of  NIRS  and  19F  MRS  of  PFOB 
techniques  could  provide  a  new  prognostic  tool  for  assessing  efficacy  of  tumor  therapy  and  could  reveal  the 
relationship  between  percent  tumor  blood  volume  and  other  tumor  kinetic  and  physiological  parameters. 
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ABSTRACT 

1  The  goal  of  this  study  was  to  investigate  the  effect  of 
photothermal  laser  irradiation  on  rat  breast  tumor  (DMBA- 
4)  vascular  contents.  An  805-nm  diode  laser  was  used  in  our 
experiment  with  a  power  density  ranging  from  0.32  to  1.27 
W/cm2.  The  dynamic  changes  of  oxygenated  hemoglobin  and 
total  hemoglobin  concentrations,  A[Hb02]  and  A[Hb]total,  in 
rat  tumors  during  photothermal  irradiation  were  noninva- 
sively  monitored  by  a  near-infrared  spectroscopy  system.  A 
multichannel  thermal  detection  system  was  also  used  simul¬ 
taneously  to  record  temperatures  at  different  locations  within 
the  tumors.  Our  experimental  results  showed  that:  (1)  pho¬ 
toirradiation  did  have  the  ability  to  induce  hyperthermic  ef¬ 
fects  inside  the  rat  breast  tumors  in  a  single  exponential 
trend;  (2)  the  significant  changes  (P  <  0.005)  of  A[Hb02]  and 
A[Hb]total  in  response  to  a  low  dosage  of  laser  irradiation  (0.32 
W/cm2)  have  a  single  exponential  increasing  trend,  similar  to 
that  seen  in  the  tumor  interior  temperature;  and  (3)  the 
increase  in  magnitude  of  A[Hb02]  is  nearly  two  times  greater 
than  that  of  A[Hb]total,  suggesting  that  photoirradiation  may 
enhance  tumor  vascular  oxygenation.  The  last  observation 
may  be  important  to  reveal  the  hidden  mechanism  of 
photoirradiation  on  tumors,  leading  to  improvement  of  tumor 
treatment  efficiency. 


INTRODUCTION 

Laser  phototherapy,  a  nonsurgical  modality  for  cancer  treatment,  is 
gaining  widespread  acceptance  because  of  precise  energy  delivery 
into  the  tumor  tissue  (1-3).  The  laser  energy  delivered  to  the 
targeted  tumors  can  induce  localized  photomechanical,  photo¬ 
chemical  and  photothermal  reactions,  thus  killing  tumor  cells.  In 
general,  photochemical  reactions  may  cause  the  change  of 
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chemical  bonds  and  form  toxic  radicals,  such  as  release  of  singlet 
oxygen,  leading  to  the  death  of  organized  tissues.  Photomechanical 
reactions  may  induce  tissue  stress,  resulting  in  the  disruption  of 
tissue  cells  and  ejection  of  material.  Photothermal  reactions  may 
induce  hyperthermia  and  coagulation,  causing  cell  destruction  (4). 
In  particular,  thermal  cytotoxicity  effects  are  found  to  be  more 
profound  under  acidotic  (low  pH)  conditions,  which  are  often 
present  in  poorly  oxygenated  tissues  (5-7),  such  as  in  tumors. 
Many  research  studies  have  reported  the  existence  of  acute/chronic 
hypoxic  regions  within  the  majority  of  solid  tumors  (8-10).  As 
a  result,  cytotoxicity  induced  by  selective  photothermal  irradiation 
can  be  more  effective  in  tumor  tissue  than  in  normal  tissue.  To 
further  enhance  the  desirable  photothermal  effects,  selective 
absorbers  are  often  used  in  the  targeted  tumor  tissue  along  with 
laser  irradiation  to  cause  selective  and  localized  photon-tissue 
interactions  (11,12).  In  addition,  administration  of  an  immunoad- 
juvant  has  been  proved  to  increase  systemic  cancer  cure  and  long¬ 
term  resistance  to  cancer  of  the  same  origin,  when  combined  with 
the  selective  photochemical  treatment  (13,14).  Chen  et  al.  have 
developed  a  laser  immunotherapy  by  combining  laser  irradiation, 
a  laser-absorbing  dye  and  an  immunoadjuvant  to  improve  cancer 
therapy  efficiency  (1 1,12,15,16).  The  previous  experimental  results 
have  demonstrated  that  this  unique  approach  has  positive  effects  on 
both  primary  and  metastatic  tumors. 

Although  the  underlying  mechanism  of  the  laser  irradiation  has 
been  intensively  studied  (17-20),  the  hemodynamic  and  physio¬ 
logical  properties  inside  the  tumors  during  phototherapy  are  poorly 
understood.  Measuring  the  reduction  of  tumor  size  and  survival 
rates  is  the  current  practice  to  evaluate  the  therapeutic  outcome  of 
therapy.  Obviously,  noninvasive,  real-time  monitoring  of  tumor 
hemodynamic  characteristics  before,  during  and  after  therapy  is 
highly  desirable.  Such  a  technique  may  help  reveal  the  tumor 
physiology  and  therapeutic  process  caused  by  the  therapy, 
providing  treatment  prognoses  and  guidance  to  optimize  the  light 
dosimetry  so  as  to  improve  the  therapeutic  outcome.  Near-infrared 
spectroscopy  (NIRS)  has  been  demonstrated  in  our  recent  studies 
to  be  such  a  noninvasive  means  in  monitoring  tumor  vascular 
oxygenation  during  therapeutic  interventions  (21-23). 

In  this  study,  we  investigated  the  effects  of  photothermal  therapy 
on  tumor  vascular  contents  and  tumor  temperature.  An  805-nm 
diode  laser  was  used  for  photothermal  irradiation  of  mammary 
tumors  in  rats.  The  dynamic  changes  caused  by  the  laser  irradiation 
in  oxygenated,  deoxygenated  and  total  hemoglobin  concentrations, 
A[Hb02],  A[Hb],  and  A[Hb]total,  in  the  tumor  vasculature  were 
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Figure  1.  Experimental  setup  for  the  simultaneous  NIRS  and  thermal 
measurements  of  breast  tumors  under  photothermal  therapy.  PMT 
represents  a  photomultiplier  tube;  I&Q  is  an  in-phase  and  quadrature- 
phase  demodulator  for  retrieving  amplitude  and  phase  information.  The 
thermal  detection  system  comprises  multiple  thermal  sensing  probes,  three 
of  which  are  inserted  into  different  regions  of  the  tumor,  and  one  probe  was 
used  to  measure  room  temperature.  The  power-adjustable  laser  light  (805 
nm)  is  vertically  delivered  to  tumor  surface  with  a  2-cm  beam  diameter  for 
photoirradiation. 

monitored  by  an  NIRS  system.  A  multichannel,  thermal  monitor 
was  also  used  simultaneously  to  record  the  temperatures  at  three  dif¬ 
ferent  locations  inside  the  breast  tumors.  The  correlations  between 
the  optical  irradiation  dosages,  tumor  temperatures,  and  tumor  vas¬ 
cular  contents  were  studied.  The  dynamic  features  of  tumor  tem¬ 
perature  and  vascular  oxygenation  were  quantified  to  reveal  the 
dynamic  effects  within  the  tumors  under  the  laser  irradiation. 

MATERIALS  AND  METHODS 

NIRS.  The  changes  of  oxygenated,  deoxygenated  and  total  hemoglobin 
concentrations,  A[Hb02],  A[Hb]  and  A[Hb]totai,  in  tumor  tissue  caused  by 
laser  irradiation  were  determined  by  a  near-infrared,  homodyne,  frequency- 
domain  spectroscopy  system  (NIM,  Philadelphia,  PA),  which  has  been 
described  in  detail  elsewhere  (21-23).  Briefly,  the  amplitude-modulated 
light  at  140  MHz  at  wavelengths  of  758  nm  and  785  nm  illuminated  the 
tumor  surface  through  a  light-delivery  fiber  bundle  with  output  powers  of 
9  mW  and  11  mW  at  785  nm  and  758  nm  wavelengths,  respectively, 
measured  at  the  tip  of  the  bundle  by  an  optical  power  meter.  The 
experimental  setup  is  shown  in  Fig.  1.  The  805 -nm  laser  light  was  vertically 
applied  on  the  tumor  surface.  The  delivering  and  detecting  fiber  bundles  of 
the  NIRS  system  were  placed  horizontally  and  in  good  contact  with  the 
surface  of  the  tumors  in  a  transmittance  mode,  without  compressing  the 
tumors.  The  diffused  light  through  the  tumor  tissue  was  collected  and 
filtered  by  a  sharp  low-pass  optical  filter,  which  cut  off  the  optical  signal 
above  805  nm  from  the  laser  irradiation,  and  then  propagated  into 
a  photomultiplier  tube  (PMT)  through  the  detecting  fiber  bundle.  The  signal 
detected  by  the  PMT  was  demodulated  through  an  in-phase  and  quadrature- 
phase  (I&Q)  circuit,  and  the  amplitude  and  phase  of  the  signal  were 
recorded. 

Algorithms  for  calculations  of  A[HhO 2/,  A [Hh],  and  A[Hh]totai.  On  the 
basis  of  the  Beer-Lambert  law,  light  attenuation,  represented  by  optical 
density  (OD),  through  a  sample  is  proportional  to  the  concentration  of 
hemoglobin  inside  the  sample  (24-28), 

OD  =  Log(A0/A)  =  8  cl  (1) 

where  Aq  and  A  are  light  intensities  of  the  incident  and  transmitted  light, 
respectively,  £  is  the  extinction  coefficient  of  hemoglobin,  c  is  the 
concentration  of  hemoglobin,  and  /  is  the  length  of  light  path  through  the 
measured  sample.  When  the  measured  sample  has  a  mixture  of  oxygenated 
and  deoxygenated  hemoglobin,  Eq.  (1)  can  be  further  evolved  to  (28-30), 

OD(A.)  =  {£HbM[Hb]  +  8Hb02(^)[Hb02]}l  (2) 

where  OD(X)  is  the  optical  density  at  wavelength  X,  £Hb(X)  and  8Hb02(^)  are 
the  extinction  coefficients  at  wavelength  X  for  molar  concentrations  of 
deoxygenated  hemoglobin,  [Hb],  and  oxygenated  hemoglobin,  [Hb02], 
respectively,  assuming  ferrihemoglobin  is  minimal.  By  using  two  wave¬ 
lengths,  both  [Hb02]  and  [Hb]  can  be  determined  by  measuring  the  OD 


values  at  the  two  specific  wavelengths,  provided  that  the  values  for  £Hb(^) 
and  £Hb02(^)  are  known: 


^Hb (A,2)OD(A,i)  —  £Hb(Al  )OD(A,2) 

^[£Hb(^2)£Hb02(^l)  —  £Hb(^l) ^02(^2)]  ’ 

^Hb02  ( 'Xl )  OD  (  A,  j )  —  £Hb02(^“l)OD(X,2) 

^[£Hb(^l)£Hb02(^2)  —  £Hb(7<2)£Hb02(^l )] 


It  follows  that  changes  in  [Hb]  and  [Hb02]  can  be  consequently  given  as: 


pjbO  ]  —  £Hb(^2)AOD(^i)  —  £Hb(^l)AOD(?l2) 

^[£Hb(^2)£Hb02(^l)  —  ^Hb (^1  )£Hb02(^2)]  ’ 

_  £m>o2(A2)AOD(Xi)  —  sHb02(Ai)AOD(A,2) 

^[£Hb(^l)£Hb02(^2)  — 


(5) 

(6) 


where  AOD(X)  represents  a  change  in  optical  density  at  the  specific 
wavelength,  X,  and  equals  log^e/Ax).  AB  and  At  correspond  to  light 
intensities  measured  under  the  baseline  and  transient  conditions. 

Note  that  in  principle,  /  represents  the  optical  path  length  between  the 
source  and  detector.  Whereas  /  is  simply  the  physical  separation,  d,  between 
the  source  and  detector  through  a  nonscattering  medium,  exact  quantifi¬ 
cation  of  /  for  an  intact  tissue  or  organ  is  complex  because  of  light  scattering 
in  tissue.  Since  /  is  in  proportion  to  the  separation,  d,  we  can  associate  /  to 
d  as  /  =  DPF*d,  where  DPF  is  a  differential  path  length  factor  to  account 
for  light  scattering.  It  has  been  well  accepted  that  together  with  DPF, 
Eq.  (2)  can  be  treated  as  modified  Beer-Lambert  law;  and  consequently, 
Eqs.  (5)  and  (6)  can  be  correctly  used  to  quantify  changes  in  [Hb]  and 
[Hb02]  in  highly  scattering  media  (29,30),  such  as  in  intact  tissue  or  organs. 

To  be  consistent  with  our  previous  work,  we  adopted  in  this  paper  the  s 
values  published  by  Zijlstra  et  al.  (31).  We  had  to  interpolate  the  £  values  at 
the  two  wavelengths  used  in  our  study,  followed  by  certain  corrections  due 
to  the  interpolation  errors  by  phantom  calibration  measurements  (22). 
Specifically,  we  included  two  factors,  pi  and  p2,  for  the  calibrated 
algorithm,  as  given  below: 


A[Hb02] 


A[Hb] 


£Ht(7pi5mn)  X  OD(758  nm)  -  £Hb(758  nm)  X  OD(785  nm) 
/[Snb(785  nm)£Hbo2(758  nm)  -  £Hb(758  nm)£Hbo2(785  nm)]  ’ 

(7) 

sag nm>  X  OD(758  nm)  -  £Hbo2(758  nm)  X  OD(785  nm) 
/[£Hb(785  nm)£Hbo2(758  nm)  -  £Hb(758  nm)£Hbo2(785  nm)]  ’ 

(8) 


where  £Hb(758  nm)  =  1.418,  £Hb02(758  nm)  =  0.6372,  £Hb(785  nm)  = 
1.111,  and  £Hb02(785  nm)  =  0.766,  all  in  mM-1cm-1.  Note  that  a  factor  of 
four  has  been  multiplied  for  each  of  the  £s  at  the  respective  wavelengths  to 
account  for  light  absorption  from  four  hemes  per  hemoglobin  molecule  (34)  ?2 
since  the  extinction  coefficients  published  in  the  field  of  biochemistry  were 
expressed  on  a  heme  basis  (24-28).  Furthermore,  we  have  used  pt  =  1.103 
and  p2  =  0.9035  according  to  our  phantom  study  (22).  After  substituting  all 
of  the  parameters  into  Eqs.  (7)  and  (8),  we  have  arrived  at  the  final 
equations  to  quantify  changes  in  hemoglobin  concentration: 

a  riiur>  i  —2.658  •  OD(758  nm)  +  3.743  •  OD(785  nm) 

A[Hb02J  — - — - ,  (9) 

r  n  2.238  •  OD(758  nm)  -  1 .683  •  OD(785  nm)  ,  N 
A[Hb]  = - - - — - - - -.  (10) 


A[Hb]total  can  also  be  obtained  by  adding  Eqs.  (9)  and  (10), 

A[Hb]totll  =  A[Hb02]  +  A[Hb] 

_  — 0.42-OD(758  nm)  +  2.06-OD(785  nm) 


Notice  that  the  units  for  A[Hb02],  A[Hb],  and  A[Hb]total  are  mM; 
however,  the  values  obtained  from  Eqs.  (9)  to  (1 1)  are  scaled  by  a  factor  of 
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DPF.  Since  DPF  is  so  far  an  unknown  parameter  for  tumors,  we  include  it 
within  the  unit  as  mM/DPF. 

Laser  photothermal  irradiation  system.  The  system  used  for  the 
photothermal  therapy  consisted  of  an  NIR  diode  laser,  DIOMED  25 
(DIOMEDICS,  The  Woodlands,  TX),  with  an  emitting  peak  wavelength  at 
805  nm  and  a  maximum  power  output  of  25  W.  The  laser  light  was  coupled 
into  an  optical  fiber,  fitted  with  a  microlens  at  the  tip  (Pioneer  Optics, 
Windsor  Locks,  CT)  to  ensure  a  uniform  beam  density,  and  delivered  onto 
the  tumor  surface.  The  laser  beam  diameter  projected  on  the  tumor  surface 
was  set  as  2  cm  in  all  the  measurements  by  using  the  visible  aiming  light  in 
the  laser  system.  After  the  beam  diameter  was  adjusted,  the  aiming  light 
was  turned  off  during  NIR  measurements. 

Multichannel  thermal  monitor.  A  multichannel  thermal  monitor  (Omega, 
Stamford,  CT),  incorporating  an  OM-700  Omega  engineering  data 
acquisition  and  control  unit,  was  used  simultaneously  with  the  NIRS 
observation  to  monitor  the  internal  temperatures  at  several  locations  within 
the  tumors.  Specifically,  three  thermal  probes  were  inserted  at  the  top 
(immediately  below  the  overlaying  skin),  middle  and  bottom  positions  of 
the  tumors,  and  one  probe  was  located  outside  the  tumor  for  room 
temperature  recording.  The  data  of  each  probe  were  recorded  continuously, 
and  the  middle  thermal  probe  recorded  the  internal  tumor  temperature  near 
the  center  of  the  optical  field  of  NIRS. 

Animal  model.  The  transplantable,  metastatic  mammary  tumor  cells 
(DMBA-4)  (33-35)  were  implanted  in  one  of  the  inguinal  fat  pads  of 
female  Wistar  Furth  rats  (Harlan  Sprague  Dawley,  Indianapolis,  IN).  Once 
the  tumors  reached  0.5  to  0.8  cm  in  diameter,  the  rats  were  anesthetized  for 
the  experiments.  Hairs  around  tumors  were  removed  for  better  probe 
contact  in  the  NIRS  measurement  and  for  more  efficient  tumor  exposure  to 
the  laser  irradiation.  The  tumor  diameters  along  the  three  major  orthogonal 
axes  (a,  h,  c )  were  measured  to  determine  the  volume  of  the  tumors  by 
using  an  ellipsoid  approximation  with  the  formula  of  V  =  (n/6)abc. 

In  this  study,  16  rat  breast  tumors  were  used  in  two  different  treatment 
protocols.  Ten  rat  tumors  were  first  exposed  to  a  single  laser  irradiation  with 
a  laser  power  of  0.32  W/cm2  (1  W  laser  power  in  a  2-cm-diameter  beam) 
for  10  min.  Three  or  four  days  later,  four  of  the  treated  tumors,  along 
with  other  six  untreated  tumors,  were  exposed  to  repeated  laser  irradia¬ 
tions  with  a  power  density  of  0.32,  0.64,  0.96  and  1.27  W/cm2  (1  W,  2  W, 
3  W,  and  4  W  in  2-cm-diameter  beam)  for  10  min  with  a  recovery  period 
of  10  min  between  two  adjacent  irradiations.  Tumor  temperatures,  tumor 
vascular  A[HbC>2],  and  A[Hb]totai  were  measured  simultaneously  in  all 
the  experiments. 

Statistic  analysis.  A  Student’s  t-test  was  applied  in  data  analysis  to 
determine  whether  the  paired  samples  are  significantly  different.  We 
considered  P  <  0.05  as  the  criteria  for  significant  difference  between  the 
two  sample  groups.  In  addition,  a  mean  and  standard  deviation  (mean  ±  sd) 
were  usually  used  to  express  the  measured  data  and  their  dispersions. 

RESULTS 

Dynamic  responses  of  A[Hb02],  A[Hb]total 
and  temperature  to  photothermal  treatment 

Using  the  experimental  setup  shown  in  Fig.  1  and  Eqs.  (9)  and 
(1 1),  tumor  vascular  A[Hb02],  A[Hb]total  and  internal  temperatures 
were  quantified  simultaneously.  Figure  2  shows  the  temporal 
profiles  of  A[Hb02],  A [Hb]total  and  the  internal  temperatures  of 
a  rat  tumor  under  photothermal  irradiation  with  a  laser  power 
density  of  0.32  W/cm2.  The  experimental  data  consisted  of  four 
temporal  periods,  as  labeled  in  the  Fig.  In  period  1,  without 
irradiation,  the  thermal  probes  recorded  the  baseline  temperatures 
at  three  different  locations  inside  the  tumor  for  about  6  min.  In 
period  2,  the  NIRS  system  was  turned  on  while  the  thermal 
recording  continued,  and  the  recorded  thermal  readings  showed  no 
obvious  perturbation  from  the  NIRS.  After  the  NIRS  signals  were 
stabilized,  a  10-min  laser  irradiation  started.  During  the  laser 
treatment,  the  tumor  temperatures,  A[Hb02],  and  A[Hb]total 
showed  significant  elevation  (period  3  in  Fig.  2a).  Tumor 
temperatures  increased  significantly  ( P  <  0.0001)  with  the 
irradiation.  The  tumor  volumes  near  the  tumor  surface  had  more 
significant  thermal  effects:  the  temperature  at  the  top  (T-top)  varied 


from  25.3  ±  0.1°C  to  30.3  ±  0.1°C  (mean  ±  sd),  while  the 
temperatures  at  the  middle  and  bottom  (T-middle  and  T-bottom)  of 
the  tumor  were  elevated  from  26.1  ±  0.1°C  to  29.8  ±  0.1°C  and 
from  27.5  ±  0.1°C  to  28.9  ±  0.1°C,  respectively.  The  maximal 
changes  in  T-top,  T-middle  and  T-bottom  caused  by  the  laser 
irradiation  were  5.0  ±  0.1°C,  3.7  ±  0.1°C  and  1.4  ±  0.1°C, 
respectively.  Similarly,  both  A[Hb02]  and  A[Hb]total  displayed 
significant  increases  over  the  entire  irradiation  period  (P  <  0.005), 
with  the  mean  elevations  of  0.04  ±  0.01  mM/DPF  and  0.02  ± 
0.01  mM/DPF,  respectively.  Notice  that  the  maximal  values  of 
A[Hb02]  and  A[Hb]total  shown  in  Fig.  2a  are  ~0.1  and  0.05  mM/ 
DPF.  When  the  laser  beam  was  turned  off  in  period  4,  A[Hb02], 
A[Hb]total  and  temperatures  all  decreased,  with  the  temperatures 
reaching  the  baselines  quickly  and  A[Hb02]  and  A[Hb]total 
declining  more  gradually  and  not  quite  reaching  the  baselines 
within  the  10-min  recovery  period. 

To  quantify  the  dynamic  behaviors  observed  in  Fig.  2a,  single¬ 
exponential  functions  were  applied  to  T-top,  T-middle,  A[Hb02], 
and  A[Hb]totai  on  both  rising  and  falling  periods,  as  shown  in  Fig. 
2b, c.  The  time  constants  obtained  from  such  fittings  reveal  the 
dynamic  responses  of  signals  to  the  initiation  and  termination  of 
photothermal  irradiation.  As  shown  in  Fig.  2b,  T-top  (time  constant 
of  4.5  ±  0.5  min)  has  a  faster  response  to  laser  irradiation  than 
T-middle  (time  constant  as  14.4  ±  0.6  min).  The  time  constants  for 
A[Hb02]  and  A[Hb]total  are  3.85  ±  0.03  min  and  3.95  ±  0.06  min, 
respectively.  In  response  to  the  termination  of  laser  irradiation, 
T-top  and  T-middle  have  a  time  constant  of  3.9  ±  0.5  min  and 
12.3  ±  0.6  min,  respectively  (Fig.  2c).  It  is  clearly  seen  that  the 
dynamic  thermal  responses  of  the  tumor  to  laser  irradiation  and  to 
the  termination  of  laser  irradiation  are  well  matched  both  at  the 
tumor  surface  and  within  the  tumor.  However,  the  time  constants 
for  A[Hb02]  and  A[Hb]total  to  return  to  the  baseline  are  much 
longer  in  the  falling  period,  with  18.5  ±  0.5  min  and  24.3  ± 
0.8  min,  respectively,  also  shown  in  Fig.  2c. 

Similar  protocols  and  measurements  were  performed  on  the 
other  nine  tumors.  Figure  3  shows  the  average  increases  of  the 
temperature  at  different  locations.  Temperatures  near  the  tumor 
surface  have  maximal  increases,  whereas  local  temperatures  near 
the  bottom  of  the  tumors  have  minimal  increases  observed  from  all 
the  10  tumors.  The  average  elevations  for  T-top  and  T-bottom  are 
7  ±  3°C  and  3  ±  2°C,  respectively.  The  thermal  readings  from  the 
middle  position  (T-middle)  represent  the  average  temperatures  in 
the  tumor  volumes  probed  by  the  NIRS  beams,  and  the  average 
increase  in  T-middle  for  the  10  rats  is  5  ±  2°C.  The  statistical 
analysis  indicates  significant  differences  between  T-top,  T-middle, 
and  T-bottom  (P  <  0.05). 

To  investigate  dynamic  behaviors  of  the  tumors  in  response  to 
the  initiation  and  termination  of  laser  irradiation,  the  time  constants 
of  T-top,  T-middle,  A[Hb02],  and  A[Hb]total  obtained  from  the 
rising  and  falling  periods  are  summarized  in  Table  1.  The  results 
demonstrate  several  points:  (1)  an  average  dynamic  change  in 
T-top  is  significantly  (P  <  0.05)  faster  than  that  in  T-middle  in  the 
rising  period,  with  a  time  constant  of  3.7  ±  1.1  min  for  T-top  and 
6.4  ±3.2  min  for  T-middle.  A  similar  significant  difference  (P  < 
0.01)  also  holds  in  the  falling  period  with  a  time  constant  of  3.8  ± 
1.2  min  for  T-top  and  6.6  ±  2.7  min  for  T-middle,  respectively.  (2) 
The  dynamic  responses  of  A[Hb02],  A[Hb]total,  and  T-middle  have 
no  significant  differences  during  the  initiation  of  photoirradiation 
(P  >  0.4).  (3)  The  average  time  constant  of  T-top  in  the  rising 
period  is  similar  (P  >  0.9)  to  that  in  the  falling  period,  and  so  is 
that  of  T-middle  (P  >  0.7).  However,  (4)  the  dynamic  response  of 
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Figure  2.  (a)  Temporal  profiles  of  the  changes  in  tumor  oxygenation  and 
total  hemoglobin  concentrations,  A[Hb02]  and  A[Hb]totai,  and  three 
intratumor  temperatures  during  photothermal  intervention  in  a  representative 
DMBA-4  rat  breast  tumor  (2.4  cm3).  The  profile  is  divided  into  four 
periods,  and  the  laser  irradiation  was  applied  only  in  period  3.  (b) 
Respective  single-exponential  curve  fittings  for  the  rising  period  of  T-top, 
T-middle,  A[Hb02]  and  A[Hb]tota|  displayed  in  (a).  The  top  .v-axis  is 
used  for  the  T-top  and  T-middle  curves,  whereas  the  bottom  .v-axis  is 
for  A[Hb02]  and  A[Hb],ota|.  Single  exponential  fitting  yields:  T-top  = 
5.76*(1  -  exp[— {f  -  147}/4.49])  +  24.95  with  R  =  0.99;  T-middle  = 


T-top 


T-middle 


T-bottom 


Positions  in  tumors 

Figure  3.  Average  changes  in  temperature  from  the  10  rats  for  T-top, 
T-middle  and  T-bottom. 


A[Hb02]  in  the  falling  period  (a  time  constant  of  28  ±  22  min)  is 
significantly  slower  ( P  <  0.01)  than  that  in  the  rising  period  (a  time 
constant  of  5.3  ±  2.8  min),  and  so  is  A[Hb]total  response,  with 
a  time  constant  of  6.2  ±  3.7  min  and  41  ±  32  min  for  the  rising 
and  falling  periods,  respectively  (P  <  0.01). 

Furthermore,  although  A[Hb02]  and  A[Flb]total  have  similar 
dynamic  trends  in  response  to  photoirradiation,  a  difference  be¬ 
tween  their  magnitudes  exists.  A  strong  correlation  (R2  =  0.95)  be¬ 
tween  A[Flb02]  and  A[Hb]totai  has  been  observed,  as  shown  in  Fig. 
4,  where  the  amplitude  in  A[Hb02]  is  nearly  twice  that  of  A[Hb]total 
for  all  10  tumors  (with  a  slope  of  the  fitted  curve  being  0.54). 


Effect  of  the  optical  power  density  on  tumor 
A[Ffb02],  A[Hb]total  and  temperatures 

Repeated  laser  irradiation  was  applied  to  the  10  tumors,  4  of  which 
were  used  in  the  previous  single-irradiation  treatment  3  or  4  days 
earlier.  Figure  5a  shows  the  profiles  of  the  tumor  vascular 
A[Hb02],  A[Hb]totai  and  tumor  temperatures  in  response  to 
irradiations  under  power  densities  of  0.32,  0.64,  0.96  and  1.27 
W/cm2,  with  an  irradiation  period  of  10  min  and  a  recovery  period 
of  10  min  between  each  treatment.  Tumor  temperatures  (T-top, 
T-middle  and  T-bottom)  increased  significantly  (P  <  0.0001)  from 
the  baseline  to  different  maximum  values  under  different  laser 
powers.  The  maximal  changes  in  temperature  in  the  middle  of  the 
tumor  (T-middle)  due  to  consecutive  irradiation  displayed  a  strong 
linear  dependence  on  the  laser  power  density,  as  shown  in  Fig.  5b. 
However,  A[Hb02]  and  A[Hb]total  showed  interesting  behaviors  in 
response  to  the  thermal  irradiation  at  different  light  dosages.  Under 
relatively  low  laser  power  irradiations  (0.32,  0.64,  0.96  W/cm2), 
both  A[Hb02]  and  A[Hb]totai  increased  significantly  (P  <  0.0001) 
with  respect  to  the  baselines,  and  they  both  declined  significantly 
under  a  high  laser  power  irradiation  (1.27  W/cm2),  as  shown  in 


7.28*(1  -  exp[-{f  -  14}/14.3])  +  25.98  with  R  =  0.99:  A[HbOJ  = 
0.093*(1  -  exp[-{t  -  14}/3.85])  +  0.013  with  R  =  0.98;  A[Hb]total  = 
0.044*(1  —  exp[— {t  —  14}/3.95])  +  0.012  with  R  =  0.98;  (c)  Respective 
single-exponential  curve  fittings  for  the  falling  period  of  T-top,  T-middle, 
A[Hb02]  and  A[Hb]totai  displayed  in  (a).  Single  exponential  fitting 
yields:  T-top  =  -4.48*(1  -  exp[-{f  -  25}/3.86])  +  29.23  with  R  =  0.98"; 
T-middle  =  -3.83*(1  -  exp[-{f  -  25}/12.29])  +  30.42  with  R  =  0.99; 
A[Hb02]  =  -0.163*(1  -  exp[-{t  -  25J/18.5])  +  0.10  with  R  =  0.94; 
A[Hb]total]  =  -0.1 14*(1  -  exp[— {f  -  25}/24.3])  +  0.06  with  R  =  0.92. 
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Table  1.  Time  constants  (minutes)  of  T-top,  Tmiddle.  A[Hb02],  and  A[Hb]total  with  respect  to  photoirradiation  and  tennination 


Tumor 

size(cm3) 

Rising  part  (in  response  to  irradiation) 

Falling  part  (in  response  to  termination) 

T-top 

Tmiddle 

A[Hb02] 

A[Hb]total 

T-top 

Tmiddle 

A[HbO,] 

A[Hb]total 

2.4 

4.49 

14.3 

3.85 

3.95 

3.86 

12.29 

18.5 

24.3 

2.5 

2.47 

6.77 

1.99 

1.53 

3.17 

8.27 

50 

74 

0.68 

4.83 

4.83 

1.76 

1.66 

3.87 

4.45 

20 

25 

0.53 

4.61 

7.62 

1.84 

1.91 

4.56 

8.37 

3.63 

2.75 

1.2 

3.28 

5.17 

6.9 

9.8 

4.20 

6.30 

80 

120 

3.3 

4.70 

6.65 

7.34 

7.83 

3.58 

4.45 

11.3 

29.9 

1.9 

1.89 

3.85 

10.0 

12.0 

1.30 

3.43 

28 

37 

3.4 

2.69 

3.13 

6.76 

7.81 

2.96 

4.58 

35 

45 

2.1 

3.25 

4.31 

7.32 

7.58 

4.52 

5.65 

23.5 

30 

3.0 

4.50 

7.80 

5.60 

8.40 

5.50 

8.30 

15.5 

23.4 

Average 

3.7 

6.4 

5.3 

6.2 

3.8 

6.6 

28. 

41 

SD* 

1.1 

3.2 

2.8 

3.7 

1.2 

2.7 

22. 

33 

*SD  =  standard  deviation. 


Fig.  5b.  Similar  to  the  feature  observed  in  single  irradiation, 
A[Hb02]  values  had  greater  magnitude  than  A[Hb]total,  whereas 
their  dynamic  trends  were  consistent. 

To  investigate  the  accumulative  effects  on  A[Hb02]  and 
A[Hb]total  from  several  thermal  irradiations  given  10  min  apart 
(Fig.  5a),  we  had  also  taken  a  direct  reading  for  A[Hb02]  and 
A[Hb]total  with  a  high  laser  power  of  1.27  W/cm2  without  any 
preirradiation.  Figure  5c  shows  an  example  of  this  case,  demon¬ 
strating  that  A[Flb02]  indeed  decreases  after  an  initial  increase  for 
a  few  minutes,  whereas  A[Flb]total  has  less  trend  to  drop  within 
the  10-min  treatment.  The  similarity  for  the  ratio  of  A[Hb02]  to 
A[Hb]totaj  close  to  a  factor  of  2  still  exists. 

Temperature  data  for  all  10  tumors  under  different  power 
irradiations  are  plotted  in  Fig.  6a.  The  maximum  changes  in  T- 
rniddle  increase  linearly  with  the  optical  radiation  dosages,  with 
a  correlation  of  R2  >  0.91  for  all  the  cases  and  having  an  average 
slope  of  13  ±  3  (°C/[W/cm2]).  To  distinguish  the  tumors  pretreated 
and  nonpretreated  with  a  single  low-power  irradiation  (0.32 
W/cm2)  3  or  4  days  before,  the  slopes  of  the  linear  regression  for 
each  of  the  tumors  shown  in  Fig.  6a  were  classified  into  two 
groups,  i.e.  without  pretreatment  and  with  pretreatment  (Fig.  6b). 
A  Student’s  ?-test  showed  that  there  is  no  significant  difference 
( P  =  0.95)  between  the  two  groups  of  slopes. 

DISCUSSION 

In  this  study,  we  have  investigated  vascular  A[Hb02],  A[Hb]total 
and  internal  temperature  in  response  to  photothermal  therapy 
simultaneously  in  rat  breast  tumors  by  using  an  NIRS  system  and 
a  thermal  monitor,  respectively.  The  results  obtained  in  this  study 
clearly  demonstrate  the  compatibility  of  the  NIRS  system  with  the 
temperature  monitoring  system,  without  noticeable  interference 
between  the  two  systems  (Fig.  2).  The  NIRS  system  offers  a  real¬ 
time  measurement  of  A[Hb02]  and  A[Hb]totai  in  tumor  vasculature 
within  the  optical  field  of  the  NIR  probes,  while  the  multichannel 
thermal  monitor  records  the  temperatures  at  specific  locations 
within  the  tumors  in  real  time.  The  compatibility  of  the  two 
systems  permits  simultaneous  determinations  of  thermal  and 
vascular  characteristics  of  treated  tumors  under  photothermal 
therapy,  providing  valuable  insight  into  the  dynamic  relationship 
between  the  tumor  vascular  contents  and  the  thermal  effect  in 
treated  tumor  tissues. 


Thermal  treatment  to  tumors  has  been  demonstrated  to  cause 
cellular  cytotoxicity,  and  thus,  to  destroy  the  tumor  cells  (36-38). 
Laser  irradiation  applied  in  this  study  proves  its  ability  to  induce 
photothermal  effects  inside  the  rat  breast  tumors.  The  single 
exponential  increase  in  tumor  temperature  obtained  in  our  results 
indicates  that  an  extended  exposure  period  to  the  laser  beam  may 
gradually  stabilize  the  tumor  temperature  at  a  certain  level, 
depending  on  both  the  laser  power  density  and  the  individual 
tumors  (as  seen  in  Figs.  2a  and  5a).  Manipulation  of  the  laser 
power  may  help  adjust  the  internal  tumor  temperature,  resulting  in 
optimal  tissue  destruction.  Although  a  strong  linear  relation 
between  the  maximum  temperature  increase  within  the  tumors 
and  the  irradiation  power  density  has  been  observed  for  all  the 
tumors,  the  actual  dependence  of  tumor  temperature  on  the  laser 
power  density  does  vary  from  tumor  to  tumor  (Fig.  6).  Such 
intertumor  variability  demonstrates  the  necessity  of  real-time 
monitoring  for  tumor  temperature  under  photothermal  therapy. 
Moreover,  the  tumor  temperatures  near  the  surface  (T-top)  respond 
to  the  photoirradiation  more  rapidly  and  with  a  greater  increase 
than  those  within  deeper  tumor  tissue  (T-middle)  (as  seen  in  Figs. 
2,  3.  and  5a,  and  Table  1).  Such  large  thermal  effects  can  cause 
tissue  damage  at  the  tumor  surface.  To  avoid  such  thermal  damage. 


Figure  4.  Relationship  between  AfHbOJ  and  A[Hb]totai  for  the  10  rat 
tumors  under  single  laser  irradiation  treatment  (0.32  W/cm2  for  10  min) 
with  a  linear  fit  of  A[Hb]total  =  0.54*A[HbO2]  —  0.0039  (the  solid  curve) 
and  a  correlation  of  R2  =  0.96.  The  dashed  line  represents  the  situation  of 
A[Hb02]  =  A[Hb],otal. 
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without  pretreatment  with  pretreatment 

Figure  6.  (a)  Correlations  between  the  maximum  changes  in  T-middle 
reached  at  the  end  of  each  photothermal  treatment  cycle  and  the  optical 
irradiation  power  density  from  the  10  treated  rats.  The  numbers  after  the 
symbol  “#”  are  labels  for  different  tumors.  The  numbers  inside  the 
parentheses  indicate  the  tumor  volume,  with  unit  of  cm3.  Note  that  there  are 
four  tumors  that  had  preirradiations  3—4  days  earlier,  and  their  tumor 
volumes  have  been  possibly  changed  in  comparison  to  those  listed  in 
Table  1.  (b).  The  slopes  in  (a)  for  all  the  tumors  were  classified  into  two 
groups:  without  pretreatment  and  with  pretreatment.  A  f-test  between  these 
two  groups  of  slopes  gives  rise  to  P  =  0.95. 


Time(mln) 

Figure  5.  (a)  The  temporal  profiles  of  vascular  A[Hb02].  A[Hb]tota|  and 
temperatures  of  a  tumor  (#6,  3.0  cm3)  in  response  to  repeated  laser 
irradiations  with  different  laser  power  densities  (0.32,  0.64,  0.96  and  1.27 
W/cm2)  for  a  10-min  treatment  cycle  with  a  10-min  recovery  period 
between  each  treatment,  (b).  The  correlations  between  the  changes  in 
A[Hb02]  (solid  triangles).  A[Hb]totai  (squares),  T-middle  (open  triangles), 
and  laser  irradiation  power  densities  in  the  tumors.  The  thick  straight  line  is 
the  best-fitting  curve  for  the  tumor  temperature  increases  as  function  of 
laser  dosage,  (c)  The  temporal  profiles  of  vascular  A[Hb02]  and  A[Hb]totai 
of  a  tumor  irradiated  with  a  single  laser  power  density  of  1.27  W/cm2. 

intratumoral  injection  of  selective  absorbers,  such  as  indocyanine 
green,  has  been  used  by  Chen  et  al.  to  selectively  manipulate  the 
temperature  within  the  tumor  without  severe  damage  at  the  tumor 
surface  (1 1,12,39). 

It  is  seen  clearly  in  this  study  that  tumor  vascular  contents,  i.e. 
oxygenated  and  deoxygenated  hemoglobin,  are  greatly  affected  by 
photothermal  therapy.  When  the  photothermal  irradiation  with 
a  low  power  density  (0.32  W/cm2)  was  applied  to  the  breast 
tumors,  both  A[HbC>2]  and  A[Hb]totai  increased  significantly  (P  < 
0.005)  in  all  the  tumors  (Figs.  2a,  4  and  5a).  The  dynamic 
behaviors  of  A[Hb02]  and  A[Hb]totai  in  response  to  photo¬ 
irradiation,  quantified  by  the  time  constants,  are  in  the  same  order 
as  those  of  tumor  interior  temperature,  T-middle,  while  their 


behaviors  in  response  to  the  termination  of  photoirradiation  are 
much  slower  than  those  of  T-middle  (Fig.  2  and  Table  1).  In 
contrast  to  the  lower  laser  power  irradiation,  higher  laser  power 
densities  could  cause  interesting  responses  in  tumor  vascular 
contents.  As  seen  in  Fig.  5b, c,  A[Hb02]  started  to  decrease  when 
the  laser  power  approached  a  certain  level,  either  with  or  without 
preirradiation.  Intriguingly,  although  A[Hb02]  has  the  same 
dynamic  trend  as  A[Hb]total  in  response  to  photoirradiation,  the 
magnitude  of  A[Hb02]  is  more  profound  than  that  of  A[Hb]total 
(Figs.  4,  5). 

The  analysis  on  the  dynamic  behaviors  of  tumor  vascular  contents 
and  interior  tumor  temperature  in  response  to  photoirradiation 
allows  us  to  better  understand  the  mechanism  of  photothermal 
therapy.  The  time  constants  of  A[Hb02]  as  well  as  of  A[Hb]total 
have  the  same  order  as  those  of  interior  tumor  temperature.  We  may 
speculate  that  an  increase  in  tumor  temperature  may  lead  to  vessel 
dilation  within  tumor  vasculature,  and  thus,  resulting  in  an  increase 
in  blood  inflow.  Such  an  increase  in  blood  inflow  will  give  rise 
to  an  increase  in  A[Hb]total  as  well  as  A[Hb02]. 

Besides  thermal  effects,  it  is  reasonable  to  expect  other 
photoactivation  processes  involved  in  the  laser  irradiation.  It  is 
known  that  direct  exposure  of  vascular  tissue  to  visible  light 
causes  vasodilation  (40),  whereas  the  thermal  effect  from  visible 
light  is  minor.  Photoactivation  has  also  been  reported  to  cause 
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upregulation  of  nitric  oxide,  and  thus,  produce  vasodilation  (41 
43).  The  increase  of  tumor  blood  contents  due  to  laser  irradiation 
observed  in  our  study  may  be  attributed  to  the  vessel  dilation 
induced  by  the  direct  effect  of  heating  as  well  as  photoactivation. 
Consistently,  the  single-exponential  trend  of  both  A[Hb]total  and 
A[Hb02]  in  response  to  laser  irradiation  indicated  the  existence 
of  vessel  dilation.  The  fact  that  the  time  constants  of  A[Hb02]  and 
A[Hb]total  in  response  to  the  laser  termination  are  longer  than  those 
in  response  to  the  laser  initiation  (Table  1,  Figs.  2a,  5a)  implies 
a  longer  settling  time  needed  for  the  recovery  of  dilated  vessels. 

In  the  laser  thermal  treatment,  we  realized  that  A[Hb02]  as  well 
as  A[Hb]totai  did  not  completely  return  to  the  baseline  during  each 
of  the  recovery  periods  (Figs.  2a,  5a).  Figure  5c  was  taken  from 
a  tumor  that  was  exposed  to  a  single  high-power  laser  irradiation 
(1.27  W/cm2),  displaying  somewhat  different  behavior  in  compar¬ 
ison  to  that  under  consecutive  repeated  treatments:  the  A[Hb02] 
signals  increased,  followed  by  an  obvious  decrease.  These  phe¬ 
nomena  lead  us  to  a  speculation  that  the  previous  laser  exposure 
“preconditions”  the  treated  tumor,  with  respect  to  consecutive 
treatments,  and  affects  the  response  of  tumor  to  subsequent  high 
irradiation  exposure,  although  the  single  lower-power  irradiation 
pretreated  3  or  4  days  earlier  has  no  obvious  influence  on  the  later 
treatment  (Fig.  6b).  It  is  reasonable  to  believe  that  the  “precondi¬ 
tioning”  can  be  minimal  if  the  time  intervals  between  two  repeated 
irradiations  are  two  to  three  times  longer  than  the  time  constants 
during  the  recovery  periods,  such  as  in  the  order  of  1—1.5  h. 

The  interesting  phenomenon,  i.e.  changes  of  oxygenated 
hemoglobin  concentration  are  nearly  two  times  greater  than  those 
of  total  hemoglobin  concentration,  allures  us  to  speculate  that  the 
increase  in  A[Hb02]  comes  partly  from  the  increase  of  A[Hb]totai, 
and  partly  from  the  conversion  of  deoxygenated  hemoglobin  to 
oxygenated  hemoglobin.  The  correlation  of  magnitude  between 
A[Hb02]  and  A[Hb]total  unambiguously  demonstrates  that  the 
increased  total  hemoglobin  concentration,  due  to  blood  vessel 
dilation,  is  mostly  oxygenated  after  the  irradiation  starts.  It 
basically  implies  that  photoirradiation  may  increase  the  tumor 
oxygenation,  which  will  be  beneficial  to  some  nonsurgical  treat¬ 
ments,  such  as  chemotherapy  and  radiotherapy  (37). 

In  summary,  changes  in  breast  tumor  temperature  and  vascular 
oxygenation  have  been  simultaneously  measured  using  a  multi¬ 
channel  thermal  monitor  and  an  NIRS  system,  while  the  tumors 
were  under  photothermal  irradiation.  The  results  have  demonstrat¬ 
ed  that:  (1)  the  compatibility  of  thermal  monitoring  system  and  the 
NIR  system  permits  simultaneous  determinations  of  thermal  and 
vascular  characteristics  of  the  treated  tumors;  (2)  photoirradiation 
did  have  the  ability  to  induce  thermal  effects  inside  the  rat  breast 
tumors  in  a  single  exponential  trend,  thus  leading  to  destruction  of 
the  tumor  cells;  (3)  the  response  of  A[Hb02]  and  A[Hb]totai  to 
photoirradiation  may  be  attributed  to  vessel  dilation  under  both 
thermal  and  photo  effect;  (4)  low-power  laser  irradiations 
predelivered  3  or  4  days  earlier  have  no  effect  on  the  later 
treatments;  the  previous  laser  exposures  within  a  time  interval  of 
10  min  do  influence  the  response  of  tumor  vascular  contents  to 
subsequent  high  irradiation;  (5)  the  fact  that  the  change  in  the 
magnitude  of  A[Hb02]  is  nearly  two  times  greater  than  that  of 
A[Hb]total  suggests  that  photoirradiation  may  enhance  tumor 
vascular  oxygenation.  The  primary  results  from  this  study  clearly 
indicate  that  besides  a  considerable  thermal  effect,  low-power 
density  irradiation  could  also  result  in  a  significant  enhancement  of 
tumor  vascular  oxygenation;  this  enhancement  may  lead  to 
improvement  of  treatment  efficiency. 
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